Kinetic mechanism of methanol decomposition on Ni(111) surface: a theoretical study

The decomposition of methanol on the Ni(111) surface has been studied with the pseudopotential method of density functional theory-generalized gradient approximation (DFT-GGA) and with the repeated slab models. The adsorption energies of possible species and the activation energy barriers of the possible elementary reactions involved are obtained in the present work. The major reaction path on Ni surfaces involves the O-H bond breaking in CH(3)OH and the further decomposition of the resulting methoxy species to CO and H via stepwise hydrogen abstractions from CH(3)O. The abstraction of hydrogen from methoxy itself is the rate-limiting step. We also confirm that the C-O and C-H bond-breaking paths, which lead to the formation of surface methyl and hydroxyl and hydroxymethyl and atom hydrogen, respectively, have higher energy barriers. Therefore, the final products are the adsorbed CO and H atom.     

氢在Ni(511)台阶面的吸附与解离研究

The adsorption and dissociation of hydrogen on stepped surface (511) of nickel are studied with the embedded-atom model (EAM) method. The adsorption energy, the length of the adsorption bond and the adsorption height for a single hydrogen atom are calculated. Three kinds of stable sites are found for hydrogen adsorption. There are the double-fold bridge site B on the step edge, the three-fold hollow site H3′ on the step surface and the four-fold hollow sites H1 and H2 on the terrace surface. Compared with a hydrogen atom adsorbed on low-index (001) surface, there are two other adsorption sites near the step: the two-fold bridge site B on the step edge and the three-fold hollow site H3′ on the step surface. At the same time, the absorbability of the hydrogen atom at the site H1 is intensified. The results show that hydrogen adsorption on Ni (511) is affected by the existence of the step. The active barriers, adsorption energy and corresponding bond length for dissociation of a hydrogen molecule on the stepped surface are presented. The results show that the dissociation is easier at the bottom of the step. It is shown that the steps are the active sites for hydrogen adsorption and dissociation.     

Surface structure and energetics of hydrogen adsorption on the Fe(111) surface

Spin-polarized density functional theory calculations have been performed to characterize the hydrogen adsorption and diffusion on the Fe(111) surface at 2/3-, 1-, and 2-monolayer (ML) coverages. It is found that the most favored adsorption site for atomic hydrogen (H) is the top-shallow bridge site (tsb), followed by the quasi 4-fold site (qff) with the energy difference of about 0.1 eV, while the top site (t) is not competitive. Furthermore, the adsorbed atomic hydrogen (H) has a high mobility, as indicated by the small diffusion barriers. The local density of state (LDOS) analysis reveals that the Fe-H (tsb or qff) bond involves mainly the Fe 4s and 4p and H 1s orbitals with less contribution of the Fe 3d orbital, while the Fe 4s, 4p, and 3d orbitals all participate in the Fe-H (top) bond. In addition, the coverage effects on the adsorption configurations and adsorption energies are addressed.     

氢原子在Pt及Pt系双金属催化表面吸附的密度泛函理论研究

采用密度泛函理论(dFT)考察了Pt(100)、(110)、(111)三种表面氢原子的吸附行为, 计算了覆盖度为0.25 ML时氢原子在Pt 三种表面和M-Pt(111)双金属(M=Al, Fe, Co, Ni, Cu, Pd)上的最稳定吸附位、表面能以及吸附前后金属表面原子层间弛豫情况. 分析了氢原子在不同双金属表面吸附前后的局域态密度变化以及双金属表面d 带中心偏离费米能级的程度并与氢吸附能进行了关联. 计算结果表明, 在Pt(100), Pt(110)和Pt(111)表面, 氢原子的稳定吸附位分别为桥位、短桥位和fcc 穴位. 三种表面中以Pt(111)的表面能最低, 结构最稳定. 氢原子在不同M-Pt(111)双金属表面上的最稳定吸附位均为fcc 穴位, 其中在Ni-Pt 双金属表面的吸附能最低, Co-Pt 次之. 表明氢原子在Ni-Pt 和Co-Pt 双金属表面的吸附最稳定. 通过对氢原子在M-Pt(111)双金属表面吸附前后的局域态密度变化的分析, 验证了氢原子吸附能计算结果的准确性. 掺杂金属Ni、Co、Fe 的3d-Pt(111)双金属表面在吸附氢原子后发生弛豫, 第一层和第二层金属原子均不同程度地向外膨胀. 此外, 3d金属的掺入使得其对应的M-Pt(111)双金属表面d带中心与Pt 相比更靠近费米能级, 吸附氢原子能力增强, 表明3d-Pt系双金属表面有可能比Pt具有更好的脱氢活性.     

A DFT study of H2O dissociation on metal‐precovered Fe (100) surface

The H₂O adsorption and dissociation on the Fe (100) surface with different precovered metals are studied by density functional theory. On both kinds of metal‐precovered surface, H₂O molecules prefer adsorb on hollow sites than bridge and top sites. The impurity energy difference is proportional to the adsorption energy, but the adsorbates are not sensitive to the adsorption orientation and height relative to the surface. The Hirshfeld charge analysis shows that water molecules act as an electron donor while the surface Fe atoms act as an electron acceptor. The rotation and dissociation of H₂O molecule occur on the Co‐ and Mn‐precovered surfaces. Some H₂O molecules are dissociated into OH and H groups. The energy barriers are about 0.5 to 1.0 eV, whose are consistence with the experimental data. H₂O molecules can be dissociated more easily at the top site on Co‐precovered surface 1 than that at bridge site on Mn‐precovered surface 2 because of the lower reaction barrier. The dispersion correction effects on the energies and adsorption configurations on Co‐precovered surface 1 were calculated by OBS + PW91. The dispersion contributions can improve a bit of the bond energy of adsorbates and weaken the hydrogen bond effect between adsorption molecules a little.     

Theoretical study of ethynylbenzene adsorption on Au(111) and implications for a new class of self-assembled monolayer

Density functional calculations of the adsorption of ethynylbenzene on the Au(111) surface show that, after cleavage of the C-H bond, the terminal carbon makes a strong covalent bond to the surface. The bond energy is shown to be about 70 kcal.mol(-1) with the fcc hollow site being most stable and the molecule oriented perpendicular to the surface. Adsorption without elimination of hydrogen is also possible via a hydrogen 1,2 shift to form a vinylidene surface-bound species, or opening of the C-C triple bond and adsorption through the two carbon atoms in a flat conformation. The reaction energy for formation of the surface-bound vinylidene is estimated to be 5 kcal.mol(-1) exothermic relative to the isolated ethynylbenzene and gold substrate.     

苯在Ag(100)表面吸附的周期性密度泛函理论研究

 采用量子化学中的密度泛函理论结合平板周期模型方法,研究了苯在Ag(100)面上的吸附方式和相对稳定性. 通过对不同吸附位置的吸附能和几何构型参数的比较发现,苯在Ag(100)表面的吸附属于较强的化学作用,穴位吸附的稳定性优于桥位,顶位吸附最不稳定. 吸附的苯分子的平衡构型发生扭曲, C-C键有较大程度的伸长; C-H键的键长基本不变,但是偏离苯环平面,并背离Ag(100)表面. 在吸附过程中,电子由苯向表面银原子转移. 本文给出了详细的轨道示意图和电荷布居分析,并且与相关的实验和理论研究结果进行了比较.     

Adsorption of cyclohexene and its dehydrogenation intermediates on nAu/Pt(100)(n＝0,1,2) surfaces:A DFT study

Adsorption of cyclohexene and its dehydrogenation intermediates on the nAu/Pt(100) (n = 0, 1, 2 means clean Pt, one monolayer and two layers of Au covered Pt surfaces, respectively.) has been investigated by self-consistent (GGA-PW91) density functional theory combined with periodic slab model. It is found that on the clean platinum, there are two kinds of favorable adsorption sites, i.e., hollow sites and bridge sites, and the adsorption energy at the hollow site is larger than that at the bridge site. However, on the Au/Pt and 2Au/Pt surfaces, there are three kinds of adsorption sites, and the adsorption energies are alike at both the bridge site and the top site. The magnitude order of the adsorption energies is as follows: clean Pt > Au/Pt > 2Au/Pt. The configurations of cyclohexene molecule have been distorted a little during the geometry optimizations. The lengths of C–M (M = Pt or Au, on the top layer of the slab) bonds are closely related to the corresponding adsorption energies.     

A Density Functional Theory Study on the Effect of Zero‐Point Energy Corrections on the Methanation Profile on Fe(100)

The thermodynamics and kinetics of the surface hydrogenation of adsorbed atomic carbon to methane, following the reaction sequence C+4 H?CH+3 H?CH₂+2 H?CH₃+H?CH₄, are studied on Fe(100) by means of density functional theory. An assessment is made on whether the adsorption energies and overall energy profile are affected when zero‐point energy (ZPE) corrections are included. The C, CH and CH₂ species are most stable at the fourfold hollow site, while CH₃ prefers the twofold bridge site. Atomic hydrogen is adsorbed at both the twofold bridge and fourfold hollow sites. Methane is physisorbed on the surface and shows neither orientation nor site preference. It is easily desorbed to the gas phase once formed. The incorporation of ZPE corrections has a very slight, if any, effect on the adsorption energies and does not alter the trends with regards to the most stable adsorption sites. The successive addition of hydrogen to atomic carbon is endothermic up to the addition of the third hydrogen atom resulting in the methyl species, but exothermic in the final hydrogenation step, which leads to methane. The overall methanation reaction is endothermic when starting from atomic carbon and hydrogen on the surface. Zero‐point energy corrections are rarely provided in the literature. Since they are derived from C? H bonds with characteristic vibrations on the order of 2500–3000 cm^?1, the equivalent ZPE of 1/2 hν is on the order of 0.2–0.3 eV and its effect on adsorption energy can in principle be significant. Particularly in reactions between CH_x and H, the ZPE correction is expected to be significant, as additional C? H bonds are formed. In this instance, the methanation reaction energy of +0.77 eV increased to +1.45 eV with the inclusion of ZPE corrections, that is, less favourable. Therefore, it is crucial to include ZPE corrections when reporting reactions involving hydrogen‐containing species.     

Theoretical Study on the Adsorption and Decomposition of Methanol over the Pt-Mo（111）/C Surface

The density functional theory(DFT) and self-consistent periodic calculation were used to investigate the methanol adsorption on the Pt-Mo(111)/C surface.The adsorption energies,equilibrium geometries and vibration frequencies of CH3OH on nine types of sites on the Pt-Mo(111)/C surface were predicted and the favorite adsorption site for methanol is the top-Pt site.Both sites of valence and conduction bands of doped system have been broadened,which are favorable for electrons to transfer to the cavity.The possible decomposition pathway was investigated with transition state searching and the calculation results indicate that the O-H bond is first broken,and then the methanol decomposes into methoxy.The activation barrier of O-H bond breaking with Pt-Mo catalyst is only 104.8 kJ mol-1,showing that carbon supported Pt-Mo alloys have promoted the decomposition of methanol.Comparing with the adsorption energies of CH3OH on the Pt(111)/C surface and that of CO,the adsorption energies of CO are higher,and Pt(111)/C is liable to be oxidized and loses the activity,which suggests that the catalyst Pt-Mo(111)/C is in favor of decomposing methanol and has better anti-poisoning ability than Pt(111)/C.     

Reaction mechanism of methanol decomposition on Pt‐based model catalysts: A theoretical study

The decomposition mechanisms of methanol on five different Pt surfaces, the flat surface of Pt(111), Pt‐defect, Pt‐step, Pt(110)(1 × 1), and Pt(110)(2 × 1), have been studied with the DFT‐GGA method using the repeated slab model. The adsorption energies under the most stable configuration of the possible species and the activation energy barriers of the possible elementary reactions involved are obtained in this work. Through systematic calculations for the reaction mechanism of methanol decomposition on these surfaces, we found that such a reaction shows the same reaction mechanism on these Pt‐based model catalysts, that is, the final products are all H (H_ads) and CO (CO_ads) via O? H bond breaking in methanol and C? H bond scission in methoxy. These results are in general agreement with the previous experimental observations. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010.     

A DFT Study of Thiophene Adsorption on the Rh（111） Surfaces

Thiophene adsorption on the Rh(111) surfaces has been investigated by density functional theory.The results show that the adsorption at the hollow and bridge sites is the most stable.The molecular plane of the thiophene ring is distorted,the C=C bond is stretched to 1.448  and the C-C bond is shortened to 1.390.The C-H bonds tilt 22～42oaway from the surface.The calculated adsorption geometries are in reasonable agreement with population analysis and density of states.The thiophene molecule obtains 0.74 electrons,reflecting the interaction between the lone pair of sulfur and the d-orbitals of metal.The reaction paths and transition states for desulfurization of the molecule have been investigated.The bridge adsorption structure of thiophene leads to a thiol via an activated reaction with an energetic barrier of 0.30 eV.This second step is slightly difficult,and dissociation into a C4H4 fragment and a sulfur atom is possible,with an energetic barrier of 0.40 eV.     

Density functional study of the interaction of chlorine atom with small neutral and charged silver clusters

Chlorine adsorption on small neutral, anionic, and cationic silver clusters Ag(n) (n=2-7) has been studied using the PW91PW91 density functional method. It was found that the adsorption of chlorine on the lowest-energy bare clusters does not always produce the lowest-energy complexes. In addition, the binding of chlorine can greatly change the geometries of the silver clusters in some cases. Among various possible adsorption sites, bridge site is energetically preferred for the neutral Ag(n) while top site is energetically more preferred for the anionic Ag(n) with n< or =6. For cationic clusters, adsorptions on bridge and face sites have similar binding energies, which are much larger than those on top sites. Natural bond orbital analyses show that irrespective of charge state, electrons always transfer from silver atoms to adsorbate and silver acts like alkali metals in the interaction with chlorine atom. Significant odd-even alternation patterns in the properties of the complexes have been observed: Even-electron clusters often have higher ionization energies, lower electron affinities, and higher dissociation energies than their odd-electron neighbors. It was also found that chlorine atoms bind more strongly with odd-electron bare clusters than with even-electron bare clusters. These patterns reveal that even-electron clusters are more stable than odd-electron clusters.     

Trapping-mediated dissociative chemisorption of cycloalkanes on Ru(001) and Ir(111): influence of ring strain and molecular geometry on the activation of C-C and C-H bonds.

We have measured the initial probabilities of dissociative chemisorption of perhydrido and perdeutero cycloalkane isotopomers on the hexagonally close-packed Ru(001) and Ir(111) single-crystalline surfaces for surface temperatures between 250 and 1100 K. Kinetic parameters (activation barrier and preexponential factor) describing the initial, rate-limiting C-H or C-C bond cleavage reactions were quantified for each cycloalkane isotopomer on each surface. Determination of the dominant initial reaction mechanism as either initial C-C or C-H bond cleavage was judged by the presence or absence of a kinetic isotope effect between the activation barriers for each cycloalkane isotopomer pair, and also by comparison with other relevant alkane activation barriers. On the Ir(111) surface, the dissociative chemisorption of cyclobutane, cyclopentane, and cyclohexane occurs via two different reaction pathways: initial C-C bond cleavage dominates on Ir(111) at high temperature (T > approximately 600 K), while at low temperature (T < approximately 400 K), initial C-H bond cleavage dominates. On the Ru(001) surface, dissociative chemisorption of cyclopentane occurs via initial C-C bond cleavage over the entire temperature range studied, whereas dissociative chemisorption of both cyclohexane and cyclooctane occurs via initial C-H bond cleavage. Comparison of the cycloalkane C-C bond activation barriers measured here with those reported previously in the literature qualitatively suggests that the difference in ring-strain energies between the initial state and the transition state for ring-opening C-C bond cleavage effectively lowers or raises the activation barrier for dissociative chemisorption via C-C bond cleavage, depending on whether the transition state is less or more strained than the initial state. Moreover, steric arguments and metal-carbon bond strength arguments have been evoked to explain the observed trend of decreasing C-H bond activation barrier with decreasing cycloalkane ring size.     

氢原子在Be(0001)表面吸附的密度泛函理论研究

采用第一性原理的密度泛函理论研究单个氢原子和多个氢原子在Be(0001)表面吸附性质.给出了氢吸附Be(0001)薄膜表面的原子结构、吸附能、饱和度、功函数、偶极修正等特性参数.同时也讨论了相关吸附性质与氢原子覆盖度(0.06-1.33ML)的关系.计算结果表明:氢原子的吸附位置与覆盖度之间有强烈的依赖关系,覆盖度低于0.67ML时,氢原子能量上易于占据fcc或hcp的中空位置;覆盖度为0.78ML时,中空位与桥位为氢原子的最佳吸附位;覆盖度在0.89到1.00ML时,桥位是氢原子吸附能量最有利的位置;以上覆盖度中Be(0001)表面最外层铍原子的结构均没有发生明显变化.当覆盖度为1.11-1.33ML,高覆盖度下Be(0001)表面的最外层铍原子部分发生膨胀,近邻氢原子渗入到铍表面次层,氢原子易于占据在hcp和桥位.吸附结构中的氢原子比氢分子中的原子稳定.当覆盖度大1.33ML时,计算结果没有发现相对于氢分子更稳定的吸氢结构.同时从分析偶极修正和氢原子吸附垂直高度随覆盖度的变化关系判断氢覆盖度为1.33ML时,在Be(0001)表面吸附达到饱和.     

Adsorption of atoms on cu surfaces: a density functional theory study

The chemisorption of atoms (H, N, S, O, and C) on Cu surfaces has been systematically studied by the density functional theory generalized gradient approximation method with the slab model. Our calculated results indicate that the orders of the adsorption energy are H < N < S < O < C on Cu(111) and H < N < O < S < C on Cu(110) and Cu(100). Furthermore, the adsorption energies of the given atoms on Cu(100) are larger than those on Cu(111) and Cu(110). The preferred adsorption sites are a 3-fold hollow site on Cu(111) and a 4-fold hollow site on Cu(100), but the preferred adsorption sites on Cu(110) are different for different adatoms. The energy, as well as the geometry, is in good agreement with the experimental and other theoretical data. In addition, this study focuses on the electronic and geometric properties of the metal-atom (M-A) bond to explain the difference in adsorption energies among adatoms. A detailed investigation of the density of states curves explains the nature of the most stable site. Finally, we test the effect of the coverage and find that the surface coverage has no influence on the preferred adsorption sites of the given adatoms on Cu(110) with the exception of hydrogen and oxygen, but has much influence on the value of the adsorption energy.     

Adsorption, mobility, and dimerization of benzaldehyde on Pt(111)

Building on results for the adsorption of benzene on Pt(111), the adsorption of benzaldehyde is investigated using density functional theory. Benzaldehyde is found to chemisorb preferentially with its aromatic ring in the flat-lying bridge geometry that is also preferred for benzene. Across the investigated geometries, adsorption is homogeneously weakened compared to corresponding benzene geometries. This is found to be true for very different adsorption modes, namely, η(6) and η(8) modes, the latter having metal atoms inserted in the carbonyl bond. Reorientation and diffusion of benzaldehyde is found to have low energy barriers. Aggregation of molecules in dimers bound by aryl C-H···O hydrogen bonds is investigated, and specific configurations are found to be up to 0.15 eV more favorable than optimally configured, separated adsorbates. The binding is significantly stronger than what is found for gas phase dimers, suggesting an enhancing effect of the metal interaction.     

Structural requirements and reaction pathways in dimethyl ether combustion catalyzed by supported Pt clusters

The identity and reversibility of the elementary steps required for catalytic combustion of dimethyl ether (DME) on Pt clusters were determined by combining isotopic and kinetic analyses with density functional theory estimates of reaction energies and activation barriers to probe the lowest energy paths. Reaction rates are limited by C-H bond activation in DME molecules adsorbed on surfaces of Pt clusters containing chemisorbed oxygen atoms at near-saturation coverages. Reaction energies and activation barriers for C-H bond activation in DME to form methoxymethyl and hydroxyl surface intermediates show that this step is more favorable than the activation of C-O bonds to form two methoxides, consistent with measured rates and kinetic isotope effects. This kinetic preference is driven by the greater stability of the CH3OCH2* and OH* intermediates relative to chemisorbed methoxides. Experimental activation barriers on Pt clusters agree with density functional theory (DFT)-derived barriers on oxygen-covered Pt(111). Measured DME turnover rates increased with increasing DME pressure, but decreased as the O2 pressure increased, because vacancies (*) on Pt surfaces nearly saturated with chemisorbed oxygen are required for DME chemisorption. DFT calculations show that although these surface vacancies are required, higher oxygen coverages lead to lower C-H activation barriers, because the basicity of oxygen adatoms increases with coverage and they become more effective in hydrogen abstraction from DME. Water inhibits reaction rates via quasi-equilibrated adsorption on vacancy sites, consistent with DFT results indicating that water binds more strongly than DME on vacancies. These conclusions are consistent with the measured kinetic response of combustion rates to DME, O2, and H2O, with H/D kinetic isotope effects, and with the absence of isotopic scrambling in reactants containing isotopic mixtures of 18O2-16O2 or 12CH3O12CH3-13CH3O13CH3. Turnover rates increased with Pt cluster size, because small clusters, with more coordinatively unsaturated surface atoms, bind oxygen atoms more strongly than larger clusters and exhibit lower steady-state vacancy concentrations and a consequently smaller number of adsorbed DME intermediates involved in kinetically relevant steps. These effects of cluster size and metal-oxygen bond energies on reactivity are ubiquitous in oxidation reactions requiring vacancies on surfaces nearly saturated with intermediates derived from O2.     

First principles analysis of the initial oxidation of Si(001) and Si(111) surfaces terminated with H and CH3

Transition state analyses have been carried out within a density functional theory setting to explain and quantify the distinctly different ways in which hydrogen and methyl terminations serve to protect silicon surfaces from the earliest onset of oxidation. We find that oxidation occurs via direct dissociative adsorption, without any energy barrier, on Si(111) and reconstructed Si(001) that have been hydrogen terminated; oxidation initiates with a barrier of only 0.05 eV on unreconstructed Si(001). The commonly measured protection afforded by hydrogen is shown to derive from a coverage-dependent dissociation rate combined with barriers to the hopping of adsorbed oxygen atoms. Methyl termination, in contrast, offers an additional level of protection because oxygen must first undergo interactions with these ligands in a three-step process with significant energy barriers: adsorption of O(2) into a C-H bond to form a C-O-O-H intermediate; decomposition of C-O-O-H into C-O-H and C=O intermediates; and, finally, hopping of oxygen atoms from ligands to the substrate.     

甲醇在FeS2(100)完整表面的吸附和分解

采用基于第一性原理的密度泛函理论结合周期平板模型方法, 研究了甲醇分子在FeS₂(100)完整表面的吸附与解离. 通过比较不同吸附位置的吸附能和构型参数发现: 表面Fe位为有利吸附位, 甲醇分子通过氧原子吸附在表面Fe位, 吸附后甲醇分子中的C―O键和O―H键都有伸长, 振动频率发生红移; 甲醇分子易于解离成甲氧基CH3O和H, 表面Fe位仍然是二者有利吸附位. 通过计算得出甲醇在FeS₂(100)表面解离吸附的可能机理: 甲醇分子首先发生O―H键的断裂, 生成甲氧基中间体, 继而甲氧基C―H键断裂, 得到最后产物HCHO和H₂.