Oxygen reduction reaction mechanism on nitrogen-doped graphene:A density functional theory study

Nitrogen-doped graphene (N-graphene) was reported to exhibit a good activity experimentally as an electrocatalyst of oxygen reduction reaction (ORR) on the cathode of fuel cells under the condition of electropotential of similar to 0.04 V (vs. NNE) and pH of 14. This material is promising to replace or partially replace the conventionally used Pt. In order to understand the experimental results. ORR catalyzed by N-graphene is studied using density functional theory (DFT) calculations under experimental conditions taking the solvent, surface adsorbates, and coverages into consideration. Two mechanisms, i.e., dissociative and associative mechanisms, over different N-doping configurations are investigated. The results show that N-graphene surface is covered by O with 1/6 monolayer, which is used for reactions in this work. The transition state of each elementary step was identified using four different approaches, which give rise to a similar chemistry. A full energy profile including all the reaction barriers shows that the associative mechanism is more energetically favored than the dissociative one and the removal of O species from the surface is the rate-determining step. (C) 2011 Elsevier Inc. All rights reserved.

Acrolein hydrogenation on Pt(211) and Au(211) surfaces:a density functional theory study

Partial hydrogenation of acrolein, the simplest alpha, beta-unsaturated aldehyde, is not only a model system to understand the selectivity in heterogeneous catalysis, but also technologically an important reaction. In this work, the reaction on Pt(211) and Au(211) surfaces is thoroughly investigated using density functional theory calculations. The formation routes of three partial hydrogenation products, namely propenol, propanal and enol, on both metals are studied. It is found that the pathway to produce enol is kinetically favoured on Pt while on Au the route of forming propenol is preferred. Our calculations also show that the propanal formation follows an indirect pathway on Pt(211). An energy decomposition method to analyze the barrier is utilized to understand the selectivities at Pt(211) and Au(211), which reveals that the interaction energies between the reactants involved in the transition states play a key role in determining the selectivity difference.

Density Functional Theory Study of Iron and Cobalt Carbides for Fischer-Tropsch Synthesis

Carbides are important phases in heterogeneous catalysis. However, the understanding of carbide phases is inadequate: Fe and Co are the two commercial catalysts for Fischer-Tropsch (FT) synthesis, and experimental work showed that Fe carbide is the active phase in FT synthesis, whereas the appearance of Co carbide is considered as a possible deactivation cause, TO understand very different catalytic roles of carbides, all the key elementary steps in FT synthesis, that is, CO dissociation, C(1) hydrogenation, and C(1)+C(1) coupling, are extensively investigated on both carbide surfaces using first principles calculations. In particular, the most important issues in FT synthesis, the activity and methane selectivity, on the carbide surfaces are quantitatively determined and analyzed. They are also discussed together with metallic Fe and Co surfaces. It is found that (i) Fe carbide is more active than metallic Fe and has similar methane selectivity to Fe, being consistent with the experiments; and (ii) Co carbide is less active than Co and has higher methane selectivity, providing evidence on the molecular level to support the suggestion that the formation of Co carbide is a cause of relatively high methane selectivity and deactivation on Co catalysts.

An Energy Descriptor To Quantify Methane Selectivity in Fischer-Tropsch Synthesis:A Density Functional Theory Study

Selectivity is a fundamental issue in heterogeneous catalysis. In this study, the CH(4) selectivity in Fischer-Tropsch synthesis is chosen to be investigated: CH4 selectivity on Rh, Co, Ru, Fe, and Re surfaces is computed by first-principles methods. In conjunction with kinetic analyses, we are able to derive the effective barrier difference between methane formation and chain growth (Delta E(eff)) to quantify the CH(4) selectivity. By using this energy descriptor, the ranking of methane selectivity predicted from density functional theory (DFT) calculations is consistent with experimental work. Moreover, a linear correlation between Delta E(eff) and the chemisorption energy of C + 4H (Delta H) is found. This fundamental finding possesses the following significance: (i) it shows that the selectivity, which appears to have kinetic characteristics, is largely determined by thermodynamic properties; and (ii) it suggests that an increase of the binding strength of C + 4H will suppress methane selectivity.

A density functional theory study of CO oxidation on Ru(0001) at low coverage

We have performed ab initio density functional theory calculations with the generalized gradient approximation to investigate CO oxidation on Ru(0001). Several reaction pathways and transition states are identified. A much higher reaction barrier compared to that on Pt(111) is determined, confirming that the Ru is very inactive for CO oxidation under UHV conditions. The origin of the reaction barrier was analyzed. It is found that in the transition state the chemisorbed O atom sits in an unfavorable bonding site and a significant competition for bonding with the same substrate atoms occurs between the CO and the chemisorbed O, resulting in the high barrier. Ab initio molecular dynamics calculations show that the activation of the chemisorbed O atom from the initial hcp hollow site (the most stable site) to the bridge site is the crucial step for the reaction. The CO oxidation on Ru(0001) via the Eley-Rideal mechanism has also been investigated. A comparison with previous theoretical work has been made. (C) 2000 American Institute of Physics. [S0021-9606(00)31223-5].

A density functional theory study of CH2 and H adsorption on Ni(111)

Ab initio total energy calculations within the density functional theory framework have been used to study the adsorption of CH2 and H as well as the coadsorption of CH2 and H on Ni(111). H binds strongly at threefold hollow sites with calculated adsorption energies of 2.60 and 2.54 eV at the face-centered-cubic (fcc) and hexagonal-close-packed (hcp) hollow sites, respectively. Adsorption energies and H-Ni distances are found to agree well with both experimental and theoretical results. CH2 adsorbs strongly at all high symmetry sites with calculated adsorption energies of 3.26, 3.22, 3.14 and 2.36 eV at the fcc, hcp, bridge and top sites, respectively. Optimized structures are reported at all sites, and, in the most stable hollow sites there is considerable internal reorganization of the CH2 fragment. The CH2 molecule is tilted, the hydrogens are inequivalent and the C-H bonds are lengthened relative to the gas phase. In the CH2-H coadsorption systems the adsorbates have a tendency to move toward bridge sites. The bonding of all adsorbates to the surface is analyzed in detail. (C) 2000 American Institute of Physics. [S0021-9606(00)71213-X].

Density functional theory study of the interaction between CO and on a Pt surface:CO/Pt(111), O/Pt(111), and CO/O/Pt(111)

Ab initio total energy calculations within the Density Functional Theory framework were carried out for Pt(111), Pt(111)-p(2x2)-CO, Pt(111)-p(2x2)-O, and Pt(111)-p(2x2)-(CO+O) to provide an insight into the interaction between CO and O on metal surfaces, an important issue in CO oxidation, and also in promotion and poisoning effects of catalysis. The geometrical structures of these systems were optimized with respect to the total energy, the results of which agree with existing experimental values very well. It is found that (i) the local structures of Pt(111)-p(2x2)-(CO+O), such as the bond lengths of C-O, C-Pt, and O-Pt (chemisorbed O atom with Pt), are almost the same as that in Pt(111)-p(2x2)-CO and Pt(111)-p(2x2)-O, respectively, (ii) the total valence charge density distributions in Pt(111)-p(2x2)-(CO+O) are very similar to that in Pt(111)-p(2x2)-CO, except in the region of the chemisorbed oxygen atom, and also nearly identical to that in Pt(111)-p(2x2)-O, apart from in the region of the chemisorbed CO, and (iii) the chemisorption energy of CO on a precovered Pt(111)-p(2x2)-O and the chemisorption energy of O on a precovered Pt(111)-p(2x2)CO are almost equal to that in Pt(111)-p(2x2)-CO and Pt(111)-p(2x2)-O, respectively. These results indicate that the interaction between CO and chemisorbed oxygen on a metal surface is mainly shore range in nature. The discussions of Pt-CO and Pt-O bonding and the interaction between CO and the chemisorbed oxygen atom on Pt(111) are augmented by local densities of states and real space distributions of quantum states.

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.

Density functional and Monte Carlo studies of sulfur. II. Equilibrium polymerization of the liquid phase

The equilibrium polymerization of sulfur is investigated by Monte Carlo simulations. The potential energy model is based on density functional results for the cohesive energy, structural, and vibrational properties as well as reactivity of sulfur rings and chains [Part I, J. Chem. Phys. 118, 9257 (2003)]. Liquid samples of 2048 atoms are simulated at temperatures 450less than or equal toTless than or equal to850 K and P=0 starting from monodisperse S-8 molecular compositions. Thermally activated bond breaking processes lead to an equilibrium population of unsaturated atoms that can change the local pattern of covalent bonds and allow the system to approach equilibrium. The concentration of unsaturated atoms and the kinetics of bond interchanges is determined by the energy DeltaE(b) required to break a covalent bond. Equilibrium with respect to the bond distribution is achieved for 15less than or equal toDeltaE(b)less than or equal to21 kcal/mol over a wide temperature range (Tgreater than or equal to450 K), within which polymerization occurs readily, with entropy from the bond distribution overcompensating the increase in enthalpy. There is a maximum in the polymerized fraction at temperature T-max that depends on DeltaE(b). This fraction decreases at higher temperature because broken bonds and short chains proliferate and, for Tless than or equal toT(max), because entropy is less important than enthalpy. The molecular size distribution is described well by a Zimm-Schulz function, plus an isolated peak for S-8. Large molecules are almost exclusively open chains. Rings tend to have fewer than 24 atoms, and only S-8 is present in significant concentrations at all T. The T dependence of the density and the dependence of polymerization fraction and degree on DeltaE(b) give estimates of the polymerization temperature T-f=450+/-20 K. (C) 2003 American Institute of Physics.

Hydrogen bonding and collective proton modes in clusters and periodic layers of squaric acid: A density functional study

Hydrogen bonding in clusters and extended layers of squaric acid molecules has been investigated by density functional computations. Equilibrium geometries, harmonic vibrational frequencies, and energy barriers for proton transfer along hydrogen bonds have been determined using the Car-Parrinello method. The results provide crucial parameters for a first principles modeling of the potential energy surface, and highlight the role of collective modes in the low-energy proton dynamics. The importance of quantum effects in condensed squaric acid systems has been investigated, and shown to be negligible for the lowest-energy collective proton modes. This information provides a quantitative basis for improved atomistic models of the order-disorder and displacive transitions undergone by squaric acid crystals as a function of temperature and pressure. (C) 2001 American Institute of Physics.