On the HSAB based estimate of charge transfer between adsorbates and metal surfaces

The applicability of the HSAB based electron charge transfer parameter, ΔN, is analyzed for molecular and atomic adsorbates on metal surfaces by means of explicit DFT calculations. For molecular adsorbates ΔN gives reasonable trends of charge transfer if work function is used for electronegativity of metal surface. For this reason, calculated work functions of low Miller index surfaces for 11 different metals are reported. As for reactive atomic adsorbates, e.g., N, O, and Cl, the charge transfer is proportional to the adatom valence times the electronegativity difference between the metal surface and the adatom, where the electronegativity of metal is represented by a linear combination of atomic Mulliken electronegativity and the work function of metal surface. It is further shown that the adatom-metal bond strength is linearly proportional to the metal-to-adatom charge transfer thus making the ΔN parameter a useful indicator to anticipate the corresponding adsorption energy trends.     

Water Adsorption and Dissociation on Ceria‐Supported Single‐Atom Catalysts: A First‐Principles DFT+U Investigation

Single‐atom catalysts have attracted wide attention owing to their extremely high atom efficiency and activities. In this paper, we applied density functional theory with the inclusion of the on‐site Coulomb interaction (DFT+U) to investigate water adsorption and dissociation on clean CeO₂(111) surfaces and single transition metal atoms (STMAs) adsorbed on the CeO₂(111) surface. It is found that the most stable water configuration is molecular adsorption on the clean CeO₂(111) surface and dissociative adsorption on STMA/CeO₂(111) surfaces, respectively. In addition, our results indicate that the more the electrons that transfer from STMA to the ceria substrate, the stronger the binding energies between the STMA and ceria surfaces. A linear relationship is identified between the water dissociation barriers and the d band centers of STMA, known as the generalized Brønsted–Evans–Polanyi principle. By combining the oxygen spillovers, single‐atom dispersion stabilities, and water dissociation barriers, Zn, Cr, and V are identified as potential candidates for the future design of ceria‐supported single‐atom catalysts for reactions in which the dissociation of water plays an important role, such as the water–gas shift reaction.     

Principle and mechanism of direct porphyrin metalation: joint experimental and theoretical investigation

The direct metalation of tetraphenylporphyrin with bare metal atoms (Co and Zn) was studied with X-ray photoelectron spectroscopy, scanning tunneling microscopy, and temperature-programmed reaction measurements on ordered monolayer films of the molecules adsorbed on a Ag(111) surface. The mechanism of this novel type of surface reaction was investigated using density functional theory (DFT) calculations for the related gas-phase reactions of the unsubstituted porphyrin with the metals Fe, Co, Ni, Cu, and Zn. The reaction starts with the formation of an initial complex, in which the metal atom is coordinated by the intact unreduced porphyrin. This complex resembles the sitting-atop complex proposed for porphyrin metalation with metal ions in solution. In two subsequent steps, the pyrrolic hydrogen atoms are transferred to the metal atom, forming H2, which is eventually released. The activation barriers of the H-transfer steps vary for the different metal atoms. DFT calculations suggest that metalations with Fe, Co, and Ni show two-state reactivity, while those with Cu and Zn proceed on a single potential energy surface. For metalation with Zn, we calculated a barrier of the first hydrogen transfer step of 32.6 kcal mol(-1), in good agreement with the overall experimental activation energy of 31 kcal mol(-1).     

Interaction between benzenedithiolate and gold: classical force field for chemical bonding

We have constructed a group of classical potentials based on ab initio density-functional theory (DFT) calculations to describe the chemical bonding between benzenedithiolate (BDT) molecule and gold atoms, including bond stretching, bond angle bending, and dihedral angle torsion involved at the interface between the molecule and gold clusters. Three DFT functionals, local-density approximation (LDA), PBE0, and X3LYP, have been implemented to calculate single point energies (SPE) for a large number of molecular configurations of BDT-1, 2 Au complexes. The three DFT methods yield similar bonding curves. The variations of atomic charges from Mulliken population analysis within the molecule/metal complex versus different molecular configurations have been investigated in detail. We found that, except for bonded atoms in BDT-1, 2 Au complexes, the Mulliken partial charges of other atoms in BDT are quite stable, which significantly reduces the uncertainty in partial charge selections in classical molecular simulations. Molecular-dynamics (MD) simulations are performed to investigate the structure of BDT self-assembled monolayer (SAM) and the adsorption geometry of S adatoms on Au (111) surface. We found that the bond-stretching potential is the most dominant part in chemical bonding. Whereas the local bonding geometry of BDT molecular configuration may depend on the DFT functional used, the global packing structure of BDT SAM is quite independent of DFT functional, even though the uncertainty of some force-field parameters for chemical bonding can be as large as approximately 100%. This indicates that the intermolecular interactions play a dominant role in determining the BDT SAMs global packing structure.     

Adsorption and Decarbonylation Reaction of Furfural on the Pd(111) and M/Pd(111)(M=Ni,Cu and Ru) Surfaces

By performing with density functional theory(DFT) method, the detailed adsorption process and the catalytic decarbonylation mechanisms of furfural over Pd(111) and M/Pd(111)(M = Ni, Cu, Ru) surfaces toward furan were clarified. The results of atomic size factor, formation energy and d-band center showed that Ru/Pd(111) surface was the most stable and active. The adsorption energies of furfural on the different surfaces followed the order Ru/Pd(111) Cu/Pd(111) Pd(111) Ni/Pd(111). After analyzing Mulliken atomic charge population and the deformation density, we can find that on Ru/Pd(111) surface, the number of charge transfer was the most and the interaction was the strongest. Therefore, its adsorption energy was the highest. Furthermore, the furfural decarbonylation pathway is more kinetically feasible on bimetallic surface, and the reaction is the most likely to occur on Ru/Pd(111).     

Initial hydrogenations of pyridine on MoP(001): a density functional study

The initial hydrogenations of pyridine on MoP(001) with various hydrogen species are studied using self-consistent periodic density functional theory (DFT). The possible surface hydrogen species are examined by studying interaction of H(2) and H(2)S with the surface, and the results suggest that the rational hydrogen source for pyridine hydrogenations should be surface hydrogen atoms, followed by adsorbed H(2)S and SH. On MoP(001), pyridine has two types of adsorption modes, i.e., side-on and end-on; and the most stable η(5)(N,C(α),C(β),C(β),C(α)) configuration of the side-on mode facilitates the hydrogenation of pyridine. The optimal hydrogenation path of pyridine with surface hydrogen atoms in the Langmuir-Hinshelwood mechanism is the formation of 3-monohydropyridine, followed by producing 3,5-dihydropyridine, in which the two-step hydrogenations take place on the C(β) atoms. When adsorbed H(2)S is considered as the source of hydrogen, slightly higher hydrogenation barriers are always involved, while the energy barriers for hydrogenations involving adsorbed SH are much lower. However, the hydrogenation of pyridine should be suppressed by the adsorption of H(2)S, and the promotion effect of adsorbed SH is limited.     

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.     

Importance of electronegativity differences and surface structure in molecular dissociation reactions at transition metal surfaces

The dissociative adsorption of N2 has been studied at both monatomic steps and flat regions on the surfaces of the 4d transition metals from Zr to Pd. Using density functional theory (DFT) calculations, we have determined and analyzed the trends in both straight reactivity and structure sensitivity across the periodic table. With regards to reactivity, we find that the trend in activation energy (Ea) is determined mainly by a charge transfer from the surface metal atoms to the N atoms during transition state formation, namely, the degree of ionicity of the N-surface bond at the transition state. Indeed, we find that the strength of the metal-N bond at the transition state (and therefore the trend in Ea) can be predicted by the difference in Mulliken electronegativity between the metal and N. Structure sensitivity is analyzed in terms of geometric and electronic effects. We find that the lowering of Ea due to steps is more pronounced on the right-hand side of the periodic table. It is found that for the early transition metals the geometric and electronic effects work in opposition when going from terrace to step active site. In the case of the late 4d metals, however, these effects work in combination, producing a more marked reduction in Ea.     

Transition state for beta-hydride elimination in alkyl groups on Pt(111)

The transition state for beta-hydride elimination in alkyl groups on the Pt(111) surface has been probed by studying the effects of fluorine substitution on the barriers to beta-hydride elimination, DeltaE++(betaH). Four different fluoroalkyl groups have been formed on the Pt(111) surface by dissociative adsorption of four fluoroalkyl iodides: RCH(2)CH(2)-I (R = CF(3), CF(3)CH(2), and CF(3)CF(2)) and (CF(3))(2)CHCH(2)-I. In the absence of preadsorbed hydrogen, fluoroalkyl groups on the Pt(111) surface dehydrogenate via beta-hydride elimination to form unsaturated fluorocarbons and deposit hydrogen atoms onto the surface. Those hydrogen atoms then hydrogenate the remaining fluoroalkyl groups to produce fluoroalkanes that desorb rapidly from the surface. The kinetics of hydrogenation and fluoroalkane desorption are rate limited by the beta-hydride elimination step and thus serve as measures of the kinetics of beta-hydride elimination. The field effects of the fluorinated substituents increase the barriers to beta-hydride elimination with a reaction constant of rho(F) = 19 +/- 2 kJ/mol. The interpretation of this effect is that the beta-carbon atom in the transition state is cationic, [RC(delta+...)H]++, with respect to the reactant. The field effect of the fluorinated substituent energetically destabilizes the electron deficient beta-carbon atom in the transition state. This is consistent with observations made on the Cu(111) surface; however, the substituent effect is significantly smaller on the Pt(111) surface. On the Pt(111) surface, the transition state for beta-hydride elimination is less polarized with respect to the initial state alkyl group than on the Cu(111) surface.     

Local HSAB principle in the conjugate addition of p‐substituted thiophenols to cyclohexenone

In the formation of new bonds through the Michel‐type addition of thiolates to a α,β‐conjugate system, it is observed that the reaction takes place faster if the thio donor compound contains an electron‐rich group or the enone acceptor is highly electron deficient. In contrast, the hard–soft acid–base (HSAB) principle predicts that this reaction is favored when a soft–soft interaction between the reactants takes place. Taking into account that softness is related to a barely electronegative atom of high polarizability, we consider it of interest to analyze the effect of charge transfer of a p‐substituent on the softness of sulfur in thiophenols, as well as its impact in the conjugate addition to 2‐cyclohexen‐1‐one. Experiments‐in‐competition, net charge of X‐groups at the p‐position of the aromatic ring, the global and local softness at sulfur, and the electrophilicity, obtained by density functional theory (DFT), led to the observation that the reaction is faster for electron attractor thiophenols. The softness at sulfur increases by delocalization of charge through the aromatic ring. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

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.     

Pt (111)面和Pt14团簇上肉桂醛吸附及选择性加氢机理研究

利用密度泛函理论研究了Pt（111）面及Pt₁₄团簇对肉桂醛（CAL）的吸附作用和不完全加氢的反应机理。分析吸附能结果表明,肉桂醛分子以C=O与C=C键协同吸附在Pt（111）面上的六角密积（Hcp）位最稳定,以C=C键吸附在Pt₁₄团簇上最稳定,且在Pt₁₄团簇上的吸附作用较Pt（111）面更强。由过渡态搜索并计算得到的反应能垒及反应热可知,肉桂醛在Pt（111）面和Pt₁₄团簇上均较容易对C=O键加氢得到肉桂醇（COL）。其中,优先加氢O原子为最佳反应路径,即Pt无论是平板还是团簇对肉桂醛加氢均有较好的选择性。同时发现,肉桂醛分子在Pt（111）面的加氢反应能垒较Pt₁₄团簇上更低,即Pt的催化活性及对肉桂醛加氢产物选择性与其结构密切相关,其中,Pt（111）面对生成肉桂醇更加有利。     

The importance of hydrogen's potential-energy surface and the strength of the forming R-H bond in surface hydrogenation reactions

An understanding of surface hydrogenation reactivity is a prevailing issue in chemistry and vital to the rational design of future catalysts. In this density-functional theory study, we address hydrogenation reactivity by examining the reaction pathways for N+H-->NH and NH+H-->NH(2) over the close-packed surfaces of the 4d transition metals from Zr-Pd. It is found that the minimum-energy reaction pathway is dictated by the ease with which H can relocate between hollow-site and top-site adsorption geometries. A transition state where H is close to a top site reduces the instability associated with bond sharing of metal atoms by H and N (NH) (bonding competition). However, if the energy difference between hollow-site and top-site adsorption energies (DeltaE(H)) is large this type of transition state is unfavorable. Thus we have determined that hydrogenation reactivity is primarily controlled by the potential-energy surface of H on the metal, which is approximated by DeltaE(H), and that the strength of N (NH) chemisorption energy is of less importance. DeltaE(H) has also enabled us to make predictions regarding the structure sensitivity of these reactions. Furthermore, we have found that the degree of bonding competition at the transition state is responsible for the trend in reaction barriers (E(a)) across the transition series. When this effect is quantified a very good linear correlation is found with E(a). In addition, we find that when considering a particular type of reaction pathway, a good linear correlation is found between the destabilizing effects of bonding competition at the transition state and the strength of the forming N-H (HN-H) bond.     

Nature of Cl bonding on the Au(111) surface: evidence of a mainly covalent interaction

We report theoretical evidence from first-principles density functional theory (DFT) calculations that the bonding between Cl and the Au(111) surface is primarily covalent in character, which is in contrast to the generally held view that the bonding of halogens to metal surfaces is ionic. We observe the transfer of charge density into the region between interacting Au and Cl atoms, which would not be expected in the case of Cl- anion formation (symmetric charge accumulation on Cl). Importantly, we also find a clear directionality of dz2 orbitals of the Au atoms pointing to the adsorbed Cl and the mixing of electronic states between the gold surface and the adsorbed chlorine.     

Charge Polarization at a Au–TiC Interface and the Generation of Highly Active and Selective Catalysts for the Low‐Temperature Water–Gas Shift Reaction

Au atoms in contact with TiC(001) undergo significant charge polarization. Strong metal–support interactions make Au/TiC(001) an excellent catalyst for the low‐temperature water–gas shift (WGS), with turnover frequencies orders of magnitude larger than those observed for conventional metal/oxide catalysts. DFT calculations indicate that the WGS reaction follows an associative mechanism with HOCO as a key intermediate.     

Experimental and theoretical investigation of the stability of Pt-3d-Pt(111) bimetallic surfaces under oxygen environment

The stability of the Pt-3d-Pt(111) (3d = Ti, V, Cr, Mn, Fe, Co, or Ni) bimetallic surface structures in the presence of adsorbed oxygen has been investigated by means of density functional theory (DFT). The dissociative binding energies of oxygen on Pt-3d-Pt(111) (i.e., subsurface 3d monolayer) and 3d-Pt-Pt(111) (i.e., surface 3d monolayer) were calculated. All of the Pt-3d-Pt(111) surfaces were found to have weaker oxygen binding energies than pure Pt(111) whereas all of the 3d-Pt-Pt(111) surfaces were found to have stronger oxygen binding energies than pure Pt(111). The total heat of reaction was calculated for the segregation for 3d metal atoms from Pt-3d-Pt(111) to 3d-Pt-Pt(111) when exposed to a half monolayer of oxygen. All of the Pt-3d-Pt(111) subsurface structures were predicted to be thermodynamically unstable with adsorbed oxygen. In addition, the segregation of subsurface Ni and Co to the surfaces of Pt-Ni-Pt(111) and Pt-Co-Pt(111) was investigated experimentally using Auger electron spectroscopy (AES) and high-resolution electron energy loss spectroscopy (HREELS). AES and HREELS confirmed the trend predicted by DFT modeling and showed that both the Pt-Ni-Pt(111) and Pt-Co-Pt(111) surface structures were unstable in the presence of adsorbed oxygen. The activation barrier of the segregation of surbsurface Ni and Co atoms was determined to be 15 +/- 2 and 7 +/- 1 kcal/mol, respectively. These results are further discussed for their implication in the design and selection of cathode bimetallic electrocatalysts for the oxygen reduction reaction (ORR) in polymer electrode membrane (PEM) fuel cells.     

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.     

Theoretical analysis of the conversion mechanism of acetylene to ethylidyne on Pt(111)

The conversion of acetylene to ethylidyne on Pt(111) has been comprehensively investigated using self-consistent periodic density functional theory. Geometries and energies for all of the intermediates involved as well as the conversion mechanism were analyzed. On Pt(111), the carbon atoms in the majority of stable C(2)H(x) (x = 1-4) intermediates prefer saturated sp(3) configurations with the missing H atoms substituted by the adjacent metal atoms. The most favorable conversion pathway for acetylene to ethylidyne is via a three-step reaction mechanism, acetylene → vinyl → vinylidene → ethylidyne. The first step, acetylene → vinyl, depends on the availability of surface H atoms: without preadsorbed H the reaction occurs via the initial disproportionation of acetylene, which resulted in adsorbed vinyl; with an abundance of preadsorbed H, acetylene could transform to vinyl via both the disproportionation and hydrogenation reactions. Conversions through initial dehydrogenation of acetylene and isomerizations of acetylene and vinyl are unfavorable due to high energy barriers along the relevant pathways. The conversion rate involving vinylidene as an intermediate is at least 100 times larger than that involving ethylidene.     

A DFT Study for NH3 Adsorption on the GaN （0001） Surface

1 INTRODUCTION Recently, the reaction of NH3 with III-V com- pounds[1～5] has attracted much attention, especially for the wurtzite GaN, since NH3 is the predominant raw stuff for growing crystalline GaN by both of the most important growth techniques[6～9], i.e., organo- metallic chemical vapor deposition (OMCVD) and molecular-beam epitaxy (MBE). Experimentally, Shekhar et al.[2] reported the chemisorption and reaction of hydrogen and ammonia on the single- crystalline GaN (0001…     

MnO2上氧气第一个电子转移步骤的从头计算研究

采用裸露簇和嵌入簇模型, 对β-MnO₂ (001), (110), (111)三个晶面以及O₂在(110)晶面的单址吸附模式(Pauling和Griffths模式), 进行从头计算. 从β-MnO₂ (001), (110), (111)三个晶面的电子结构差异以及O₂在(110)晶面吸附的吸附能、几何结构、集居数以及净电荷数分析得到: (001), (110), (111)三个晶面中(110)晶面的催化活性最高, 其活性顺序为(110)＞(111)＞(001). 氧气在(110)晶面的吸附, Pauling和Griffths两种吸附模式均存在, 属于化学吸附中的离子吸附. 氧气与MnO₂固体间发生了单电子转移, 氧气得到电子被还原成O₂^-, 转移电子属于整个体系, 具有离域性.