Chemisorption of (CHx and C2Hy) hydrocarbons on Pt(111) clusters and surfaces from DFT studies

We used the B3LYP flavor of density functional theory (DFT) to study the chemisorption of all CH(x) and C(2)H(y) intermediates on the Pt(111) surface. The surface was modeled with the 35 atom Pt(14.13.8) cluster, which was found to be reliable for describing all adsorption sites. We find that these hydrocarbons all bind covalently (sigma-bonds) to the surface, in agreement with the studies by Kua and Goddard on small Pt clusters. In nearly every case the structure of the adsorbed hydrocarbon achieves a saturated configuration in which each C is almost tetrahedral with the missing H atoms replaced by covalent bonds to the surface Pt atoms. Thus, (Pt(3))CH prefers a mu(3) hollow site (fcc), (Pt(2))CH(2) prefers a mu(2) bridge site, and PtCH(3) prefers mu(1) on-top sites. Vinyl leads to (Pt(2))CH-CH(2)(Pt), which prefers a mu(3) hollow site (fcc). The only exceptions to this model are ethynyl (CCH), which binds as (Pt(2))C=CH(Pt), retaining a CC pi-bond while binding at a mu(3) hollow site (fcc), and HCCH, which binds as (Pt)HC=CH(Pt), retaining a pi bond that coordinates to a third atom of a mu(3) hollow site (fcc) to form an off center structure. These structures are in good agreement with available experimental data. For all species we calculated heats of formation (DeltaH(f)) to be used for considering various reaction pathways on Pt(111). For conditions of low coverage, the most strongly bound CH(x) species is methylidyne (CH, BE = 146.61 kcal/mol), and ethylidyne (CCH(3), BE = 134.83 kcal/mol) among the C(2)H(y) molecules. We find that the net bond energy is nearly proportional to the number of C-Pt bonds (48.80 kcal/mol per bond) with the average bond energy decreasing slightly with the number of C ligands.     

An interaction model for OH + H2O-mixed and pure H2O overlayers adsorbed on Pt(111)

A model potential for the adsorbate-adsorbate interaction among OH and H2O molecules adsorbed on a Pt(111) surface has been developed solely based on first-principle calculations. By combining this directional-dependent model potential for the lateral interaction with a lattice model of Ising type, large length scale structure calculations can be made. The strength of different hydrogen bonds can be analyzed in detail from this model potential. It is found that the hydrogen bond between OH and H2O molecules is stronger than that between two H2O molecules (0.4 eV per pair as compared to 0.2 eV per pair, respectively). Via the computed chemical potential for water in mixed OH + H2O overlayers the water uptake as a function of oxygen precoverage on Pt(111) has been determined. The results compare very well with recent experiments.     

Adsorption and decomposition of cyclohexanone (C6H10O) on Pt(111) and the (2 × 2) and (√3 × √3)r30°-Sn/Pt(111) surface alloys

Adsorption and decomposition of cyclohexanone (C(6)H(10)O) on Pt(111) and on two ordered Pt-Sn surface alloys, (2 × 2)-Sn/Pt(111) and (√3 × √3)R30°-Sn/Pt(111), formed by vapor deposition of Sn on the Pt(111) single crystal surface were studied with TPD, HREELS, AES, LEED, and DFT calculations with vibrational analyses. Saturation coverage of C(6)H(10)O was found to be 0.25 ML, independent of the Sn surface concentration. The Pt(111) surface was reactive toward cyclohexanone, with the adsorption in the monolayer being about 70% irreversible. C(6)H(10)O decomposed to yield CO, H(2)O, H(2), and CH(4). Some C-O bond breaking occurred, yielding H(2)O and leaving some carbon on the surface after TPD. HREELS data showed that cyclohexanone decomposition in the monolayer began by 200 K. Intermediates from cyclohexanone decomposition were also relatively unstable on Pt(111), since coadsorbed CO and H were formed below 250 K. Surface Sn allowed for some cyclohexanone to adsorb reversibly. C(6)H(10)O dissociated on the (2 × 2) surface to form CO and H(2)O at low coverages, and methane and H(2) in smaller amounts than on Pt(111). Adsorption of cyclohexanone on (√3 × √3)R30°-Sn/Pt(111) at 90 K was mostly reversible. DFT calculations suggest that C(6)H(10)O adsorbs on Pt(111) in two configurations: by bonding weakly through oxygen to an atop Pt site and more strongly through simultaneously oxygen and carbon of the carbonyl to a bridged Pt-Pt site. In contrast, on alloy surfaces, C(6)H(10)O bonds preferentially to Sn. The presence of Sn, furthermore, is predicted to make the formation of the strongly bound C(6)H(10)O species bonding through O and C, which is a likely decomposition precursor, thermodynamically unfavorable. Alloying with Sn, thus, is shown to moderate adsorptive and reactive activity of Pt(111).     

Comparative study of water dissociation on Rh(111) and Ni(111) studied with first principles calculations

The dissociation and formation of water on the Rh(111) and Ni(111) surfaces have been studied using density functional theory with generalized gradient approximation and ultrasoft pseudopotentials. Calculations have been performed on 2x2 surface unit cells, corresponding to coverages of 0.25 ML, with spot checks on 3x3 surface unit cells (0.11 ML). On both surfaces, the authors find that water adsorbs flat on top of a surface atom, with binding energies of 0.35 and 0.25 eV, respectively, on Rh(111) and Ni(111), and is free to rotate in the surface plane. Barriers of 0.92 and 0.89 eV have to be overcome to dissociate the molecule into OH and H on the Rh(111) and Ni(111) surfaces, respectively. Further barriers of 1.03 and 0.97 eV need to be overcome to dissociate OH into O and H. The barriers for the formation of the OH molecule from isolated adsorbed O and H are found to be 1.1 and 1.3 eV, and the barriers for the formation of the water molecule from isolated adsorbed OH and H are 0.82 and 1.05 eV on the two surfaces. These barriers are found to vary very little as coverage is changed from 0.25 to 0.11 ML. The authors have also studied the dissociation of OH in the presence of coadsorbed H or O. The presence of a coadsorbed H atom only weakly affects the energy barriers, but the effect of O is significant, changing the dissociation barrier from 1.03 to 1.37 and 1.15 eV at 0.25 or 0.11 ML coverage on the Rh(111) surface. Finally, the authors have studied the dissociation of water in the presence of one O atom on Rh(111), at 0.11 ML coverage, and the authors find a barrier of 0.56 eV to dissociate the molecule into OH+OH.     

Thermodynamics and kinetics of oxygen-induced segregation of 3d metals in Pt-3d-Pt(111) and Pt-3d-Pt(100) bimetallic structures

The stability of subsurface 3d transition metals (3d represents Ni, Co, Fe, Mn, Cr, V, and Ti) in Pt(111) and Pt(100) was examined in vacuum and with 0.5 ML atomic oxygen by a combined experimental and density functional theory (DFT) approach. DFT was used to predict the trends in the binding energy of oxygen and in the stability of 3d metals to remain in the subsurface layer. DFT calculations predicted that for both (111) and (100) crystal planes the subsurface Pt-3d-Pt configurations were thermodynamically preferred in vacuum and that the surface 3d-Pt-Pt configurations were preferred with the adsorption of 0.5 ML atomic oxygen. Experimentally, the DFT predictions were verified by using Auger electron spectroscopy to monitor the segregation of Ni and Co in Pt-3d-Pt structures on polycrystalline Pt foil, composed of mainly (111) and (100) facets. The activation barrier for the oxygen-induced segregation of Ni was found to be 17+/-1 kcal/mol attributed to the Pt(111) areas and 27+/-1 kcal/mol attributed to the Pt(100) areas of the Pt foil. For Pt-Co-Pt, the activation barrier was found to be 10+/-1 kcal/mol and was attributed to the Pt(111) areas of the Pt foil. The Bronsted-Evans-Polanyi relationship was utilized to predict the activation barriers for segregation of the other Pt-3d-Pt(111) and Pt-3d-Pt(100) systems. These results are further discussed in connection to the activity and stability for cathode bimetallic electrocatalysts for proton exchange membrane fuel cells.     

Quantum chemical study of adsorption and dissociation of H2S on the gallium-rich GaAs (001)-4 x 2 surface

H(2)S adsorption and dissociation on the gallium-rich GaAs(001)-4 x 2 surface is investigated using hybrid density functional theory. Starting from chemisorbed H(2)S on the GaAs(001)-4 x 2 surface, two possible reaction routes have been proposed. We find that H(2)S adsorbs molecularly onto GaAs(001)-4 x 2 via the formation of a dative bond, and this process is exothermic with adsorption energy of 6.6 kcal/mol. For the first reaction route, one of the H atoms from the chemisorbed H(2)S is transferred to a second-layer As atom and the dissociated SH is inserted into the Ga-As bond with an activation barrier of 8.2 kcal/mol, which is found to be 29.3 kcal/mol more stable than the reactants. For the second case, the dissociated species may insert themselves into the Ga-Ga dimer resulting in the Ga-H-Ga and Ga-HS-Ga bridge-bonded states, which are found to be 29.8 and 22.2 kcal/mol more stable than the reactants, respectively. However, the calculations also show that the activation barrier (16.1 kcal/mol) for chemisorbed H(2)S dissociation through the second route is higher than the transfer of one H atom into a second-layer As atom. As a result, we conclude that sulfur insertion into the Ga-As bond is more kinetically favorable.     

Mechanism for degradation of Nafion in PEM fuel cells from quantum mechanics calculations

We report results of quantum mechanics (QM) mechanistic studies of Nafion membrane degradation in a polymer electrolyte membrane (PEM) fuel cell. Experiments suggest that Nafion degradation is caused by generation of trace radical species (such as OH(●), H(●)) only when in the presence of H(2), O(2), and Pt. We use density functional theory (DFT) to construct the potential energy surfaces for various plausible reactions involving intermediates that might be formed when Nafion is exposed to H(2) (or H(+)) and O(2) in the presence of the Pt catalyst. We find a barrier of 0.53 eV for OH radical formation from HOOH chemisorbed on Pt(111) and of 0.76 eV from chemisorbed OOH(ad), suggesting that OH might be present during the ORR, particularly when the fuel cell is turned on and off. Based on the QM, we propose two chemical mechanisms for OH radical attack on the Nafion polymer: (1) OH attack on the S-C bond to form H(2)SO(4) plus a carbon radical (barrier: 0.96 eV) followed by decomposition of the carbon radical to form an epoxide (barrier: 1.40 eV). (2) OH attack on H(2) crossover gas to form hydrogen radical (barrier: 0.04 eV), which subsequently attacks a C-F bond to form HF plus carbon radicals (barrier as low as 1.00 eV). This carbon radical can then decompose to form a ketone plus a carbon radical with a barrier of 0.86 eV. The products (HF, OCF(2), SCF(2)) of these proposed mechanisms have all been observed by F NMR in the fuel cell exit gases along with the decrease in pH expected from our mechanism.     

CO在Pt低指数面上吸附行为的理论研究

用密度泛函理论(DFT)中的广义梯度近似(GGA)方法对CO-Pt低指数面吸附体系进行了结构优化, 并对吸附体系的吸附热、C—O键和C—Pt键的键长、布居数分析、电子态密度进行了研究. 计算结果表明, 在0.25 ML(monolayer)的覆盖率下, CO最容易在Pt(100)晶面的桥位、Pt(110)晶面的短桥位、Pt(111)晶面的hcp三重位吸附, 吸附热分别达到了2.11、2.37、1.96 eV; CO在吸附成键过程中伴有电子在CO分子和Pt之间的转移. 吸附后, C—O键被削弱, 键长变长, 金属内部的作用亦被削弱, 其表层Pt 原子的布居数明显降低; 态密度分析表明, CO在吸附过程中, 其4σ、1π、5σ、2π轨道均参与成键.     

Sulfate-enhanced catalytic destruction of 1,1,1-trichlorethane over Pt(111)

The catalytic destruction of 1,1,1-trichloroethane (TCA) over model sulfated Pt(111) surfaces has been investigated by fast X-ray photoelectron spectroscopy and mass spectrometry. TCA adsorbs molecularly over SO4 precovered Pt(111) at 100 K, with a saturation coverage of 0.4 monolayer (ML) comparable to that on the bare surface. Surface crowding perturbs both TCA and SO4 species within the mixed adlayer, evidenced by strong, coverage-dependent C 1s and Cl and S 2p core-level shifts. TCA undergoes complete dechlorination above 170 K, accompanied by C-C bond cleavage to form surface CH3, CO, and Cl moieties. These in turn react between 170 and 350 K to evolve gaseous CO2, C2H6, and H2O. Subsequent CH3 dehydrogenation and combustion occurs between 350 and 450 K, again liberating CO2 and water. Combustion is accompanied by SO4 reduction, with the coincident evolution of gas phase SO2 and CO2 suggesting the formation of a CO-SOx surface complex. Reactively formed HCl desorbs in a single state at 400 K. Only trace (<0.06 ML) residual atomic carbon and chlorine remain on the surface by 500 K.     

Hydrogen fluoride adsorption and reaction on the α-Al2O3(0001) surface: a density functional theory study

The adsorption and reaction behaviors of HF on the α-Al(2)O(3)(0001) surface are systematically investigated using density functional theory method. By increasing the number of HF molecules in a p(2 × 1) α-Al(2)O(3)(0001) slab, we find that HF is chemically dissociated at low coverage; while both physical and dissociative adsorption occurs at a 3/2 monolayer (ML) coverage. At the same coverage (1.0 ML), diverse configurations of the dissociated HF are obtained in the p(2 × 1) model; while only one is observed in the p(1 × 1) slab due to its smaller surface area compared with the former one. Preliminary fluorination reaction study suggests that the total energy of two dissociated HF in the p(2 × 1) slab increases by 1.00 and 0.72 eV for the formation and desorption of water intermediate, respectively. The coadsorption behaviors of HF and H(2)O indicate that the pre-adsorbed water is unfavorable for the fluorination of Al(2)O(3), which is well consistent with the experimental results. The calculated density of states show that the peak of σ(H-F) disappears, while the peaks of σ(H-O) and σ(Al-F) are observed at -8.4 and -5 to -3 eV for the dissociated HF. Charge density difference analysis indicates that the dissociated F atom attracts electrons, while no obvious changes on electrons are observed for the surface Al atoms.     

Energetics of cyclohexene adsorption and reaction on Pt(111) by low-temperature microcalorimetry

The heat of adsorption and sticking probability of cyclohexene on Pt(111) were measured as a function of coverage using single-crystal adsorption calorimetry in the temperature range from 100 to 300 K. At 100 K, cyclohexene adsorbs as intact di-sigma bonded cyclohexene on Pt(111), and the heat of adsorption is well described by a second-order polynomial (130 - 47 theta - 1250 theta(2)) kJ/mol, yielding a standard enthalpy of formation of di-sigma bonded cyclohexene on Pt(111) at low coverages of -135 kJ/mol and a C-Pt sigma bond strength of 205 kJ/mol. At 281 K, cyclohexene dehydrogenates upon adsorption, forming adsorbed 2-cyclohexenyl (c-C6H(9,a)) and adsorbed hydrogen, and the heat of adsorption is well described by another second-order polynomial (174 - 700 theta + 761 theta(2)) kJ/mol. This yields a standard enthalpy of formation of adsorbed 2-cyclohexenyl on Pt(111) at a low coverage of -143 kJ/mol. At coverages below 0.10 ML, the sticking probability of cyclohexene on Pt(111) is close to unity (>0.95), independent of temperature.     

Heat of adsorption of naphthalene on Pt(111) measured by adsorption calorimetry

The heat of adsorption of naphthalene on Pt(111) at 300 K was measured with single-crystal adsorption calorimetry. The heat of adsorption on the ideal, defect-free surface is estimated to be (300 - 34 - 199(2)) kJ/mol. From this, a C-Pt bond energy for aromatic hydrocarbons on Pt(111) of approximately 30 kJ/mol is estimated, consistent with earlier results for benzene on Pt(111). There is higher heat of adsorption at very low coverage, attributed to step sites where the adsorption heat is >/=330 kJ/mol. Saturation coverage, = 1 ML, corresponds to 1.55 x 10(14) molecules/cm(2). Sticking probability measurements of naphthalene on Pt(111) give a high initial value of 1.0 and a Kisliuk-type coverage dependence that implies precursor-mediated sticking. The ratio of the hopping rate to the desorption rate of this precursor is approximately 51. Naphthalene adsorbs transiently on top of chemisorbed naphthalene molecules with a heat of adsorption of 83-87 kJ/mol.     

Theoretical calculation of the dehydrogenation of ethanol on a Rh/CeO2(111) surface

We applied periodic density-functional theory (DFT) to investigate the dehydrogenation of ethanol on a Rh/CeO2 (111) surface. Ethanol is calculated to have the greatest energy of adsorption when the oxygen atom of the molecule is adsorbed onto a Ce atom in the surface, relative to other surface atoms (Rh or O). Before forming a six-membered ring of an oxametallacyclic compound (Rh-CH2CH2O-Ce(a)), two hydrogen atoms from ethanol are first eliminated; the barriers for dissociation of the O-H and the beta-carbon (CH2-H) hydrogens are calculated to be 12.00 and 28.57 kcal/mol, respectively. The dehydrogenated H atom has the greatest adsorption energy (E(ads) = 101.59 kcal/mol) when it is adsorbed onto an oxygen atom of the surface. The dehydrogenation continues with the loss of two hydrogens from the alpha-carbon, forming an intermediate species Rh-CH2CO-Ce(a), for which the successive barriers are 34.26 and 40.84 kcal/mol. Scission of the C-C bond occurs at this stage with a dissociation barrier Ea = 49.54 kcal/mol, to form Rh-CH(2(a)) + 4H(a) + CO(g). At high temperatures, these adsorbates desorb to yield the final products CH(4(g)), H(2(g)), and CO(g).     

CaO对Pt/C催化硬脂酸脱羧反应性能的影响

以硬脂酸(C₁₇H₃₅COOH)脱羧为探针反应,研究了亚临界水中CaO对Pt/C催化脱羧反应性能的影响.实验表明,CaO可以促进脱羧反应,对十七烷选择性没有影响.CaO/硬脂酸物质的量比为0.5时,在330 ℃反应1 h,硬脂酸转化率由未添加CaO时的46.06%提高到66.99%.硬脂酸催化脱羧的最佳反应温度为350 ℃,高于这一温度时,烷基链上碳碳键断裂的副反应增加,导致十七烷选择性降低.亚临界水中硬脂酸脱羧反应符合一级动力学,建立的动力学方程可以较好地预测不同反应条件下十七烷的产率.根据实验结果推测,氧化钙与硬脂酸反应生成硬脂酸钙,催化剂表面的吸附态羧酸盐增加,从而提高了脱羧反应的速率.     

1-Phenyl-1-propyne on Cu(111): TOFMS TPD, XPS, UPS, and 2PPE Studies

The bonding properties of 1-phenyl-1-propyne (PP, C6H5CCCH3) on Cu(111) at 100 K have been studied using temperature-programmed desorption (TPD), and X-ray, ultraviolet, and two-photon photoemission spectroscopies (XPS, UPS, and 2PPE). In TPD, there is no evidence for dissociation. Multilayer desorption occurs at 187 K, and monolayer desorption occurs at 320 (83.5 kJ/mol) and 390 K (102.4 kJ/mol), with the latter dominating. Based on the calibrated C(1s) XPS, the saturation monolayer coverage is one PP per four surface Cu atoms. The broad and asymmetric C(1s) intensity profile of the monolayer can be resolved into three symmetric components, with peaks at 283.6, 284.5, and 285.2 eV and intensities of 2:6:1, respectively. These are attributed, respectively, to acetylenic carbons bound to Cu, phenyl, and methyl carbons. The monolayer valence band ultraviolet photoemission spectrum profile contains four resonances attributable to PP perturbed by interactions with the Cu(111) substrate. With the exception of the highest occupied molecular orbital (HOMO) that is shifted by 0.4 eV, these are uniformly shifted by 1 eV further from the Fermi level for the multilayer. Calculated electron density plots of the occupied orbitals coupled with UPS profiles suggest a spectator role for the phenyl group and bonding to Cu via the acetylenic carbons. The adsorption of 1.0 monolayer (ML) of PP on Cu(111) lowers the work function by 0.85 eV. Using 2PPE, two unoccupied orbitals were identified at 1.0 (U1*-LUMO) and 0.6 eV (U2*-image state) below the vacuum level. A chemisorption model consistent with these spectroscopic results and the major chemisorption peak in TPD involve di-sigma-bonding of the acetylenic carbons to a pair of second-nearest neighbor surface Cu atoms (cross-bridge).     

Undissociated versus dissociated structures for water clusters and ammonia-water clusters: (H2O)n and NH3(H2O)n-1 (n = 5, 8, 9, 21). Theoretical study

To understand the autoionization of pure water and the solvation of ammonia in water, we investigated the undissociated and dissociated (ion-pair) structures of (H2O) n and NH3(H2O)n-1 (n = 5, 8, 9, 21) using density functional theory (DFT) and second order Moller-Plesset perturbation theory (MP2). The stability, thermodynamic properties, and infrared spectra were also studied. The dissociated (ion-pair) form of the clusters tends to favor the solvent-separated ion-pair of H3O+/NH4+ and OH-. As for the NH3(H2O)20 cluster, the undissociated structure has the internal conformation, in contrast to the surface conformation for the (H2O)21 cluster, whereas the dissociated structure of NH3(H2O)20 has the surface conformation. As the cluster size of (H2O)n/NH3(H2O)n-1 increases, the difference in standard free energy between undissociated and dissociated (ion-pair) clusters is asymptotically well corroborated with the experimental free energy change at infinite dilution of H3O+/NH4+ and OH-. The predicted NH and OH stretching frequencies of the undissociated and dissociated (ion-pair) clusters are discussed.     

Theoretical study of CCl(4) adsorption and hydrogenation on a Pt (111) surface

The adsorption and hydrogenation of carbon tetrachloride (CCl(4)) on a Pt (111) surface have been investigated using density functional theory (DFT). We have performed calculations on the adsorption energies and structures of CCl(4) on four different adsorption sites of a Pt (111) surface using the full adsorbate geometry optimization method. The results show that the adsorption energy of all of the potential sites is less than -17 kcal/mol, which indicates that CCl(4) is physiosorbed on a Pt (111) surface through van der Waals interactions. The dissociation and hydrogenation pathways were investigated by a transition state search. For the Pt(15), Pt(19), and Pt(25) cluster surfaces, the activation energies of dissociation obtained in this work are 15.69, 16.94, and 16.77 kcal/mol, respectively. The hydrogenation of CCl(3). was studied at the on-top site of the Pt(15) cluster, and the calculated activation energy is 5.06 kcal/mol. The small activation energies indicate that the Pt (111) surface has high catalytic activity for the CCl(4) hydrogenation reaction. In addition, the Hirshfeld population analysis reveals that the charge transfer from the Pt (111) surface to the adsorbates occurs in both the dissociation and hydrogenation pathways.     

Reduction of the (001) surface of gamma-V2O5 compared to alpha-V2O5

The defect-free gamma-V(2)O(5)(001) surface and ordered structures of oxygen vacancies have been studied for a wide range of defect concentrations, Theta ((1)/(6) monolayer (ML) < or = Theta < or = 1 ML), combining density functional theory and statistical thermodynamics. The gamma polymorph of V(2)O(5) is characterized by two structurally different vanadium sites, V(A) and V(B). The V(A) sites having a weaker bond to an adjacent crystal layer are easier to reduce. Up to (1)/(2) ML, the V(A) defect structures with defects aligned along the [010] direction are increasingly more stable as in alpha-V(2)O(5)(001). At higher defect concentrations, the different coordination of the V(B) vanadium atoms at the gamma-V(2)O(5) surface causes an increase in the vacancy formation energy of approximately 0.8 eV/atom at Theta = 1.0 compared to Theta = (1)/(2). For alpha-V(2)O(5), this increase amounts to 0.2 eV/atom only. Under conditions (low oxygen partial pressures and high temperatures) at which the alpha-V(2)O(5)(001) surface would be fully reduced, the gamma-V(2)O(5)(001) surface is only partially reduced. The presence of surface vanadyl oxygen groups at V(B) sites may change the surface reactivity compared to that of alpha-V(2)O(5)(001).     

Coverage dependence and hydroperoxyl-mediated pathway of catalytic water formation on Pt (111) surface

Hydrogen oxidation on Pt (111) surface is modeled by density functional theory (DFT). Previous DFT calculations showed too large O2 dissociation barriers, but we find them highly coverage dependent: when the coverage is low, dissociation barriers close to experimental values (approximately 0.3 eV) are obtained. For the whole reaction, a new pathway involving hydroperoxyl (OOH) intermediate is found, with the highest reaction barrier of only approximately 0.4 eV. This may explain the experimental observation of catalytic water formation on Pt (111) surface above the H2O desorption temperature of 170 K, despite that the direct reaction between chemisorbed O and H atoms is a highly activated process with barrier approximately 1 eV as previous calculations showed.