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

The dissociative adsorption of N-2 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.

Insight into the adsorption competition and the relationship between dissociation and association reactions in ammonia synthesis

Ammonia synthesis on three metal surfaces (Zr, Ru, and Pd) is investigated using density functional theory calculations. In addition to N-2 dissociation, all the transition states of the hydrogenation reactions from N to NH3 are located and the reaction energy profiles at both low and high surface coverages are compared and analyzed. The following are found: (i) Surface coverage effect on dissociation reactions is more significant than that on association reactions. (ii) The difference between N and H chemisorption energies, the so-called chemisorption energy gap which is a measure of adsorption competition, is vital to the reactivity of the catalysts. (iii) The hydrogenation barriers can considerably affect the overall rate of ammonia synthesis. A simple model to describe the relationship between dissociation and association reactions is proposed. (c) 2007 American Institute of Physics.

An insight into alkali promotion: a density functional theory study of CO dissociation on K/Rh(111)

The important role of alkali additives in heterogeneous catalysis is, to a large extent, related to the high promotion effect they have on many fundamental reactions. The wide application of alkali additives in industry does not, however, reflect a thorough understanding of the mechanism of their promotional abilities. To investigate the physical origin of the alkali promotion effect, we have studied CO dissociation on clean Rh(111) and K-covered Rh(111) surfaces using density functional theory. By varying the position of potassium atoms relative to a dissociating CO, we have mapped out the importance of different K effects on the CO dissociation reactions. The K-induced changes in the reaction pathways and reaction barriers have been determined; in particular, a large reduction of the CO dissociation barrier has been identified. A thorough analysis of this promotion effect allows us to rationalize both the electronic and the geometrical factors that govern alkali promotion effect: (i) The extent of barrier reductions depends strongly on how close K is to the dissociating CO. (ii) Direct K-O bonding that is in a very short range plays a crucial role in reducing the barrier. (iii) K can have a rather long-range effect on the TS structure, which could reduce slightly the barriers.

General trends in CO dissociation on transition metal surfaces

Dissociative adsorption is one of the most important reactions in catalysis. In this communication we propose a model aiming to generalize the important factors that affect dissociation reactions. Specifically, for a dissociation reaction, say AB -->A + B, the model connects the dissociation barrier with the association barrier, the chemisorption energies of A and B at the final state and the bonding energy of AB in the gas phase. To apply this model, we have calculated CO dissociation on Ru(0001), Rh(111), Pd(111) (4d transition metals), Os(0001), Ir(111), and Pt(111) (5d transition metals) using density function theory (DFT). All the barriers are determined. We find that the DFT results can be rationalized within the model. The model can also be used to explain many experimental observations. (C) 2001 American Institute of Physics.

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.

The dissociation of molecularly adsorbed CO and CN over the 4d transition metals: A universal relationship between the reaction barriers and the reaction enthalpies

Density functional theory calculations are used to study the stability of molecularly adsorbed CO and CN over transition metal surfaces. The minimum energy reaction pathways, corresponding reaction barriers (E-a), and reaction enthalpies (Delta H) for the dissociation of CO and CN to atomic products over the 4d transition metals from Zr to Pd have been determined. CO is found to be the more stable adsorbate on the right hand side of the period (from Tc onwards), whereas CN is the more stable surface species on the early metals (Zr, Nb and Mo). A single linear relationship is found to exist that correlates the barriers of both reactions with their respective reaction enthalpies. (c) 2006 Elsevier B.V. All rights reserved.

Reactivity of the 4d transition metals toward N hydrogenation and NH dissociation: A DFT-based HSAB analysis

The density functional theory (DFT) based hard-soft acid-base (HSAB) reactivity indices, including the electrophilicity index, have been successfully applied to many areas of molecular chemistry. In this work we test the applicability of such an approach to fundamental surface chemistry. We have considered, as prototypical surface reactions, both the hydrogenation of atomic nitrogen and the dissociative adsorption of the NH molecular radical. By use of a DFT methodology, the minimum energy reaction pathways, and corresponding reaction barriers, of the above reactions over Zr(001), Nb(110), Mo(110), Tc(001), Ru(001), Rh(111), and Pd(111) have been determined. By consideration of the chemical potential and chemical hardness of the surface metal atoms, and the principle of electronegativity equalization, it is found that the charge transferred to the NH radical during the process of dissociative adsorption correlates very well with that determined by Mulliken population analysis. Furthermore, it is found that the stability of the NH/surface transition state complex relates directly to this charge transfer and that the trend in transition state stability predicted by a HSAB; treatment correlates very strongly with that determined by DFT calculations. With regards to N hydrogenation, we find that during the course of the reaction, H loses cohesion to the surface, as it must migrate from a 3-fold hollow site to either a bridge or top site, to react with N. Partial density of states (PDOS) and Mulliken population analysis reveal that this loss of bonding is accompanied by charge transfer from H to the surface metal atoms. Moreover, by simple modeling, we show that the reaction barriers are directly proportional to this mandatory charge transfer. Indeed, it is found that the reaction barriers correlate very well with the electrophilicity index of the metal atoms.

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.

Stepwise addition reactions in ammonia synthesis: A first principles study

Catalytic ammonia synthesis is believed to proceed via dissociation of N-2 and H-2 with subsequent stepwise addition reactions from an adsorbed nitrogen atom to NH3. The first step, N-2 dissociation, has been thoroughly studied. However, little is known about the microscopic details of the stepwise addition reactions. To shed light on these stepwise addition reactions, density functional theory calculations with the generalized gradient approximation are employed to investigate NHx (x=1,3) formation on Ru(0001). Transition states and reaction barriers are determined in each elementary step. It is found that the reaction barriers for 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. In addition, one of the stepwise addition reactions on a stepped surface is also considered. The reaction barrier is found to be much higher than that of N-2 dissociation on the same stepped surface, which indicates the importance of stepwise addition reactions in ammonia synthesis. (C) 2001 American Institute of Physics.

Universal tight binding model for chemical reactions in solution and at surfaces. II. Water

A revised water model intended for use in condensed phase simulations in the framework of the self consistent polarizable ion tight binding theory is constructed. The model is applied to water monomer, dimer, hexamers, ice, and liquid, where it demonstrates good agreement with theoretical results obtained by more accurate methods, such as DFT and CCSD(T), and with experiment. In particular, the temperature dependence of the self diffusion coefficient in liquid water predicted by the model, closely reproduces experimental curves in the temperature interval between 230 K and 350 K. In addition, and in contrast to standard DFT, the model properly orders the relative densities of liquid water and ice. A notable, but inevitable, shortcoming of the model is underestimation of the static dielectric constant by a factor of two. We demonstrate that the description of inter and intramolecular forces embodied in the tight binding approximation in quantum mechanics leads to a number of valuable insights which can be missing from ab initio quantum chemistry and classical force fields. These include a discussion of the origin of the enhanced molecular electric dipole moment in the condensed phases, and a detailed explanation for the increase of coordination number in liquid water as a function of temperature and compared with ice-leading to insights into the anomalous expansion on freezing. The theory holds out the prospect of an understanding of the currently unexplained density maximum of water near the freezing point. 

Kinetics study on thermal dissociation of levoglucosan during cellulose pyrolysis

(Chemical Equation Presented) The mechanisms and kinetics studies of the levoglucosan (LG) primary decomposition during cellulose pyrolysis have been carried out theoretically in this paper. Three decomposition mechanisms (C-O bond scission, C-C bond scission, and LG dehydration) including nine pathways and 16 elementary reactions were studied at the B3LYP/6-31 + G(D,P) level based on quantum mechanics. The variational transi-tion- state rate constants for every elementary reaction and every pathway were calculated within 298-1550 K. The first-order Arrhenius expressions for these 16 elementary reactions and nine pathways were suggested. It was concluded that computational method using transition state theory (TST) without tunneling correction gives good description for LG decomposition by comparing with the experimental result. With the temperature range of 667-1327 K, one dehydration pathway, with one water molecule composed of a hydrogen atom from C3 and a hydroxyl group from C2, is a preferred LG decomposition pathway by fitting well with the experimental results. The calculated Arrhenius plot of C-O bond scission mechanism is better agreed with the experimental Arrhenius plot than that of C-C bond scission. This C-O bond scission mechanism starts with breaking of C1-O5 and C6-O1 bonds with formation of CO molecule (C1-O1) simultaneously. C-C bond scission mechanism is the highest energetic barrier pathway for LG decomposition. © 2013 Elsevier Ltd. All rights reserved.

Catalytic Water Formation on Platinum: A First-Principles Study

The study of catalytic behavior begins with one seemingly simple process, namely the hydrogenation of O to H2O on platinum. Despite the apparent simplicity its mechanism has been much debated. We have used density functional theory with,gradient corrections to examine microscopic reaction pathways for several elementary steps implicated in this fundamental catalytic process. We find that H2O formation from chemisorbed O and H atoms is a highly activated process. The largest barrier along this route, with a value of similar to1 eV, is the addition of the first H to O to produce OH. Once formed, however, OH groups are easily hydrogenated to H2O with a barrier of similar to0.2 eV. Disproportionation reactions with 1:1 and 2:1 stoichiometries of H2O and O have been examined as alternative routes for OH formation. Both stoichiometries of reaction produce OH groups with barriers that are much lower than that associated with the O + H reaction. H2O, therefore, acts as an autocatalyst in the overall H O formation process. Disproportionation with a 2:1 stoichiometry is thermodynamically and kinetically favored over disproportionation with a l:I stoichiometry. This highlights an additional (promotional) role of the second H2O molecule in this process. In support of our previous suggestion that the key intermediate in the low-temperature H2O formation reaction is a mixed OH and H2O overlayer we find that then is a very large barrier for the dissociation of the second H2O molecule in the 2:1 disproportionation process. We suggest that the proposed intermediate is then hydrogenated to H2O through a very facile proton transfer mechanism.