General rules for predicting where a catalytic reaction should occur on metal surfaces: A density functional theory study of C-H and C-O bond breaking/making on flat, stepped, and kinked metal surfaces

To predict where a catalytic reaction should occur is a fundamental issue scientifically. Technologically, it is also important because it can facilitate the catalyst's design. However, to date, the understanding of this issue is rather limited. In this work, two types of reactions, CH4 CH3 + H and CO C + 0 on two transition metal surfaces, were chosen as model systems aiming to address in general where a catalytic reaction should occur. The dissociations of CH4 - CH3 + H and CO --> C + O and their reverse reactions on flat, stepped, and kinked Rh and Pd surfaces were studied in detail. We find the following: First, for the CH4 Ch(3) + H reaction, the dissociation barrier is reduced by similar to0.3 eV on steps and kinks as compared to that on flat surfaces. On the other hand, there is essentially no difference in barrier for the association reaction of CH3 + H on the flat surfaces and the defects. Second, for the CO C + 0 reaction, the dissociation barrier decreases dramatically (more than 0.8 eV on Rh and Pd) on steps and kinks as compared to that on flat surfaces. In contrast to the CH3 + H reaction, the C + 0 association reaction also preferentially occurs on steps and kinks. We also present a detailed analysis of the reaction barriers in which each barrier is decomposed quantitatively into a local electronic effect and a geometrical effect. Our DFT calculations show that surface defects such as steps and kinks can largely facilitate bond breaking, while whether the surface defects could promote bond formation depends on the individual reaction as well as the particular metal. The physical origin of these trends is identified and discussed. On the basis of our results, we arrive at some simple rules with respect to where a reaction should occur: (i) defects such as steps are always favored for dissociation reactions as compared to flat surfaces; and (ii) the reaction site of the association reactions is largely related to the magnitude of the bonding competition effect, which is determined by the reactant and metal valency. Reactions with high valency reactants are more likely to occur on defects (more structure-sensitive), as compared to reactions with low valency reactants. Moreover, the reactions on late transition metals are more likely to proceed on defects than those on the early transition metals.

A Density Functional Theory Study on the Active Center of Fe-Only Hydrogenase: Characterization and Electronic Structure of the Redox States

We have carried out extensive density functional theory (DFT) calculations for possible redox states of the active center in Fe-only hydrogenases. The active center is modeled by [(H(CH(3))S)(CO)(CN(-))Fe(p)(mu-DTN)(mu-CO)Fe(d)(CO)(CN(-))(L)](z) (z is the net charge in the complex; Fe(p)= the proximal Fe, Fe(d) = the distal Fe, DTN = (-SCH(2)NHCH(2)S-), L is the ligand that bonds with the Fed at the trans position to the bridging CO). Structures of possible redox states are optimized, and CO stretching frequencies are calculated. By a detailed comparison of all the calculated structures and the vibrational frequencies with the available experimental data, we find that (i) the fully oxidized, inactive state is an Fe(II)-Fe(II) state with a hydroxyl (OH(-)) group bonded at the Fe(d), (ii) the oxidized, active state is an Fe(II)-Fe(l) complex which is consistent with the assignment of Cao and Hall (J. Am. Chem. Soc. 2001, 123, 3734), and (iii) the fully reduced state is a mixture with the major component being a protonated Fe(l)-Fe(l) complex and the other component being its self-arranged form, Fe(II)-Fe(II) hydride, Our calculations also show that the exogenous CO can strongly bond with the Fe(II)-Fe(l) species, but cannot bond with the Fe(l)-Fe(l) complex. This result is consistent with experiments that CO tends to inhibit the oxidized, active state, but not the fully reduced state. The electronic structures of all the redox states have been analyzed. It is found that a frontier orbital which is a mixing state between the e(g) of Fe and the 2pi of the bridging CO plays a key role concerning the reactivity of Fe-only hydrogenases: (1) it is unoccupied in the fully oxidized, inactive state, half-occupied in the oxidized, active state, and fully occupied in the fully reduced state; (ii) the e(g)-2pi orbital is a bonding state, and this is the key reason for stability of the low oxidation states, such as Fe(l)-Fe(l) complexes; and (iii) in the e(g)-2pi orbital more charge accumulates between the bridging CO and the Fe(d) than between the bridging CO and the Fe(p), and the occupation increase in this orbital will enhance the bonding between the bridging CO and the Fe(d), leading to the bridging-CO shift toward the Fe(d).

Reactions of polycarbonate with cyclohexene oxide and phosphites: A density functional study

Epoxides and phosphites are often used as additives to stabilize the properties of polymers, including bisphenol A polycarbonate (BPA-PC). We describe density functional (DF) calculations of the reactions of cyclohexene oxide (CHO, cyclohexane epoxide) and phosphites with chain segments of BPA-PC, with the aim of identifying possible reaction paths and energy barriers. The reactions of CHO with the OH-terminated PC chains and with the carbonate group are exothermic, although there is an energy barrier in each case of more than 10 kcal/mol. A comparison of results for different CHO isomers demonstrates the importance of steric effects. The reactions between the same groups of the PC chain and the phosphites 2-[2,4-bis(tert-butyl)phenoxy]-5,5-dimethyl-1,3,2-dioxaphosphorinane] (BPDD) and trimethyl phosphite (TMP), and their phosphonate isomers are characterized by large energy barriers.

Understanding the dehydrogenation mechanism of tetrahydrocarbazole over palladium using a combined experimental and density functional theory approach

The mechanism of the dehydrogenation of tetrahydrocarbazole to carbazole over palladium has been examined for the first time. By use of a combination of deuterium exchange experiments and density functional theory calculations, a detailed reaction profile for the aromatization of tetrahydrocarbazole has been identified and validated by experiment. As with many dehydrogenation reactions, the initial hydrogen abstraction is found to have the highest reaction barrier. Tetrahydrocarbazole has four hydrogens which can, in principle, be cleaved initially; however, the theory and experiment show that the reaction is dominated by the cleavage of the carbon hydrogens at the carbon atoms in positions 1 and 4. The two pathways originating from these two C-H bond cleavage processes are found to have similar reaction energy profiles and both contribute to the overall reaction.

Au on (111) and (110) surfaces of CeO2: A density-functional theory study

The atomic structures of gold supported on (111) and (110) surfaces of CeO2 have been studied using density-functional theory calculations. A single Au atom is placed on three adsorption sites on the surfaces; the stoichiometric surfaces, an oxygen vacancy and a Ce-vacancy. It is found that (i) the Au adsorption energies are in the following order: E-ad(Ce-vacancy) > E-ad(O-vacancy) > E-ad(stoichiometric surface); and (ii) the Au atom adsorption on the Ce-vacancy activates O atoms nearby. One 0 atom is less stable than that in O-2 in the gas phase and another O atom is much easier to remove compared to that of the stoichiometric surfaces. These results suggest that the Au adsorption on Ce-vacancies not only creates an O-vacancy but also activates an O atom nearby. This provides a piece of direct evidence that Au adsorption on a Ce-vacancy may be responsible for some unique catalytic properties of Au/CeO2. (C) 2008 Elsevier B.V. All rights reserved.

Insight into the solvent effect: A density functional theory study of cisplatin hydrolysis

The solvent effect on reactions in solutions is crucial for many systems. In this study, the reaction barrier with respect to the number of solvent molecules included in the system is systematically studied using density function theory calculations. Our results show that the barriers rapidly converge with respect to the number of solvent molecules. The solvent effect is investigated by calculating cisplatin hydrolysis in several types of solvents. The results are analyzed and a linear relationship between the reaction barrier and the interaction strength of solvent-reactants is found. Insight into the general solvent effect is obtained. (c) 2006 American Institute of Physics.

A density functional theory study of the a-olefin selectivity in Fischer-Tropsch synthesis

We studied the alpha-olefin selectivity in Fischer-Tropsch (FT) synthesis using density functional theory (131717) calculations. We calculated the relevant elementary steps from C-2 to C-6 species. Our results showed that the barriers of hydrogenation and dehydrogenation reactions were constant with different chain lengths, and the chemisorption energies of alpha-olefins from DFT calculations also were very similar, except for C-2 species. A simple expression of the paraffin/olefin ratio was obtained based on a kinetic model. Combining the expression of the paraffin/olefin ratio and our calculation results, experimental findings are satisfactorily explained. We found that the physical origin of the chain length dependence of paraffin/olefin ratio is the chain length dependence of both the van der Waals interaction between adsorbed alpha-olefins and metal surfaces and the entropy difference between adsorbed and gaseous alpha-olefins, and that the greater chemisorption energy of ethylene is the main reason for the abnormal ethane/ethylene ratio. (c) 2008 Elsevier Inc. All rights reserved.

A density functional theory study of sulfur poisoning

Density functional theory calculations have been used to investigate the chemisorption of H, S, SH, and H2S as well as the hydrogenation reactions S+H and SH+H on a Rh surface with steps, Rh(211), aiming to explain sulfur poisoning effect. In the S hydrogenation from S to H2S, the transition state of the first step S+H-->SH is reached when the S moves to the step-bridge and H is on the off-top site. In the second step, SH+H-->H2S, the transition state is reached when SH moves to the top site and H is close to another top site nearby. Our results show that it is difficult to hydrogenate S and they poison defects such as steps. In order to address why S is poisoning, hydrogenation of C, N, and O on Rh(211) has also been calculated and has been found that the reverse and forward reactions possess similar barriers in contrast to the S hydrogenation. The physical origin of these differences has been analyzed and discussed. (C) 2005 American Institute of Physics.

General insight into CO oxidation: A density functional theory study of the reaction mechanism on platinum oxides

CO oxidation on PtO2(110) has been studied using density functional theory calculations. Four possible reaction mechanisms were investigated and the most feasible one is the following: (i) the O at the bridge site of PtO2(110) reacts with CO on the coordinatively unsaturated site (CUS) with a negligible barrier; (ii) O-2 adsorbs on the bridge site and then interacts with CO on the CUS to form an OO-CO complex; (iii) the bond of O-OCO breaks to produce CO2 with a small barrier (0.01 eV). The CO oxidation mechanisms on metals and metal oxides are rationalized by a simple model: The O-surface bonding determines the reactivity on surfaces; it also determines whether the atomic or molecular mechanism is preferred. The reactivity on metal oxides is further found to be related to the 3rd ionization energy of the metal atom.

CO oxidation and NO reduction on metal surfaces: density functional theory investigations

This article reviews the accumulated theoretical results, in particular density functional theory calculations, on two catalytic processes, CO oxidation and NO reduction on metal surfaces. Owing to their importance in automotive emission control, these two reactions have generated a lot of interest in the last 20 years. Here the pathways and energetics of the involved elementary reactions under different catalytic conditions are described in detail and the understanding of the reactions is generalized. It is concluded that density functional theory calculations can be applied to catalysis to elucidate mechanisms of complex surface reactions and to understand the electronic structure of chemical processes in general. The achieved molecular knowledge of chemical reactions is certainly beneficial to new catalyst design.

Catalytic role of gold in gold-based catalysts: A density functional theory study on the CO oxidation on gold

Gold-based catalysts have been of intense interests in recent years, being regarded as a new generation of catalysts due to their unusually high catalytic performance. For example, CO oxidation on Au/TiO2 has been found to occur at a temperature as low as 200 K. Despite extensive studies in the field, the microscopic mechanism of CO oxidation on Au-based catalysts remains controversial. Aiming to provide insight into the catalytic roles of Au, we have performed extensive density functional theory calculations for the elementary steps in CO oxidation on Au surfaces. O atom adsorption, CO adsorption, O-2 dissociation, and CO oxidation on a series of Au surfaces, including flat surfaces, defects and small clusters, have been investigated in detail. Many transition states involved are located, and the lowest energy pathways are determined. We find the following: (i) the most stable site for O atom on Au is the bridge site of step edge, not a kink site; (ii) O-2 dissociation on Au (O-2-->20(ad)) is hindered by high barriers with the lowest barrier being 0.93 eV on a step edge; (iii) CO can react with atomic O with a substantially lower barrier, 0.25 eV, on Au steps where CO can adsorb; (iv) CO can react with molecular O-2 on Au steps with a low barrier of 0.46 eV, which features an unsymmetrical four-center intermediate state (O-O-CO); and (v) O-2 can adsorb on the interface of Au/TiO2 with a reasonable chemisorption energy. On the basis of our calculations, we suggest that (i) O-2 dissociation on Au surfaces including particles cannot occur at low temperatures; (ii) CO oxidation on Au/inactive-materials occurs on Au steps via a two-step mechanism: CO+O-2-->CO2+O, and CO+O-->CO2; and (iii) CO oxidation on Au/active-materials also follows the two-step mechanism with reactions occurring at the interface.

A density functional theory study of stepwise addition reactions in ammonia synthesis on Ru(0001)

To shed light on stepwise addition reactions in ammonia synthesis, density functional theory calculations are carried out to investigate NHx (x = 1-3) formation on Ru(0001). The reactions on a flat surface are first examined. Transition states and reaction barriers are determined. It is found that the reaction barriers for these stepwise addition reactions are rather high. For example, the barrier for NH hydrogenation is calculated to be 1.28 eV, which is comparable with that of N-2 dissociation. One of the stepwise addition reactions, NH + H --> NH2, on a stepped surface is also considered. Interestingly, the reaction barrier is found to be significantly lower than that on the flat surface, but is considerably higher than that of N-2 dissociation on the same stepped surface. In addition, the coverage effect on the reaction energetics is also addressed. (C) 2001 Published by Elsevier Science B.V.

A density functional theory study of carbon monoxide oxidation on the Cu3Pt(111) alloy surface: Comparison with the reactions on Pt(111) and Cu(111)

Alloying metals is often used as an effective way to enhance the reactivity of surfaces. Aiming to shed light on the effect of alloying on reaction mechanisms, we carry out a comparative study of CO oxidation on Cu3Pt(111), Pt(111), and Cu(111) by means of density functional theory calculations. Alloying effects on the bonding sites and bonding energies of adsorbates, and the reaction pathways are investigated. It is shown that CO preferentially adsorbs on an atop site of Pt and O preferentially adsorbs on a fcc hollow site of three Cu atoms on Cu3Pt(111). It is also found that the adsorption energies of CO (or O-a) decreases on Pt (or Cu) on the alloy surface with respect to those on pure metals. More importantly, having identified the transition states for CO oxidation on those three surfaces, we found an interesting trend for the reaction barrier on the three surfaces. Similar to the adsorption energies, the reaction barrier on Cu3Pt possesses an intermediate value of those on pure Pt and Cu metals. The physical origin of these results has been analyzed in detail. (C) 2001 American Institute of Physics.

A density functional theory study of the surface relaxation and reactivity of Cu2O(100)

Density functional theory has been used to investigate the surface relaxation of Cu2O(100) and the adsorption of NO. The calculations indicate the formation of surface copper dimers on relaxation coupled with a large contraction of the spacing between the first and second layers. Local density of states for atoms in the top three layers shows that the third layer copper atoms have the greatest change in bonding character. Adsorption energies have been calculated for the N-down and O-down adsorption of NO on the Cu2O(100) surface. These indicate that N-down adsorption is favoured and that in this case NO-lattice oxygen interactions dominate the adsorbate structure. (C) 2000 Elsevier Science B.V. All rights reserved.

Reaction Mechanisms of Crotonaldehyde Hydrogenation on Pt(111): Density Functional Theory and Microkinetic Modeling

The microkinetics based on density function theory (DFT) calculations is utilized to investigate the reaction mechanism of crotonaldehyde hydrogenation on Pt(111) in the free energy landscape. The dominant reaction channel of each hydrogenation product is identified. Each of them begins with the first surface hydrogenation of the carbonyl oxygen of crotonaldehyde on the surface. A new mechanism, 1,4-addition mechanism generating enols (butenol), which readily tautomerize to saturated aldehydes (butanal), is identified as a primary mechanism to yield saturated aldehydes instead of the 3,4-addition via direct hydrogenation of the ethylenic bond. The calculation results also show that the full hydrogenation product, butylalcohol, mainly stems from the deep hydrogenation of surface open-shell dihydrogenation intermediates. It is found that the apparent barriers of the dominant pathways to yield three final products are similar on P(111), which makes it difficult to achieve a high selectivity to the desired crotyl alcohol (COL).