Pathways for methanol steam reforming involving adsorbed formaldehyde and hydroxyl intermediates on Cu(111): density functional theory studies

Plane-wave density functional theory calculations have been carried out to explore possible pathways in methanol steam reforming (MSR) on Cu(111). We focus on reactions involving the adsorbed formaldehyde intermediate (CH(2)O) produced by methanol decomposition and the surface hydroxyl (OH) species generated by dissociative adsorption of H(2)O. Several possible pathways leading to the H(2) + CO(2) products have been identified. The two most likely pathways involve the formate (CHOO), rather than the carboxyl (COOH), intermediate, and they possess barriers lower than that of the rate-limiting step of MSR, namely the dehydrogenation of adsorbed methoxyl (CH(3)O) species.

A dimensionless reaction coordinate for quantifying the lateness of transition states

The Hammond-Leffler postulate asserts that transition states of exothermic reactions are reactant-like (early), whereas transition states of endothermic reactions are product-like (late). Related postulates have been proposed to describe the sensitivity of activation barriers for reactions occurring on catalytic surfaces to the catalyst structure. To evaluate the validity of these postulates for different chemical reactions, a general method for classifying transition states as either early or late is needed. One can envision a dimensionless reaction coordinate that changes continuously and monotonically from 0 to 1 along a minimum energy reaction pathway. The value of the dimensionless reaction coordinate for the transition state (W(TS)) classifies transition states as (a) early when W(TS) < 0.5, (b) late when W(TS) > 0.5, and (c) equidistant between reactants and products when W(TS) = 0.5. In this article, we derive such a dimensionless reaction coordinate and illustrate its usefulness for several different chemical reactions.

On the mechanism of low-temperature water gas shift reaction on copper

Periodic, self-consistent density functional theory (DFT-GGA) calculations are used to investigate the water gas shift reaction (WGSR) mechanism on Cu(111). The thermochemistry and activation energy barriers for all the elementary steps of the commonly accepted redox mechanism, involving complete water activation to atomic oxygen, are presented. Through our calculations, we identify carboxyl, a new reactive intermediate, which plays a central role in WGSR on Cu(111). The thermochemistry and activation energy barriers of the elementary steps of a new reaction path, involving carboxyl, are studied. A detailed DFT-based microkinetic model of experimental reaction rates, accounting for both the previous and the new WGSR mechanism show that, under relevant experimental conditions, (1) the carboxyl-mediated route is the dominant path, and (2) the initial hydrogen abstraction from water is the rate-limiting step. Formate is a stable "spectator" species, formed predominantly through CO2 hydrogenation. In addition, the microkinetic model allows for predictions of (i) surface coverage of intermediates, (ii) WGSR apparent activation energy, and (iii) reaction orders with respect to CO, H2O, CO2, and H2.

Total Oxidation of Methanol on Cu(110): A Density Functional Theory Study

The partial and total oxidation of methanol on clean and oxygen-precovered Cu(110) has been studied by periodic density functional theory calculations within the generalized gradient approximation. Reaction paths including the geometry and the energetics of several reaction intermediates and the activation barriers between them have been determined, thus creating a complete scheme for methanol oxidation on copper. The calculations demonstrate that the specific structure of oxygen on copper plays an important role in both the partial and the total oxidation of methanol. For lower oxygen concentrations on the surface, the partial oxidation of methanol to formaldehyde is promoted by the presence of oxygen on the surface through the removal of hydrogen in the form of water, which prevents the recombinative desorption of methanol. At larger oxygen concentrations, the presence of isolated oxygen atoms reduces the C-H bond breaking barrier of adsorbed methoxy considerably, thus accelerating the formation of formaldehyde. Furthermore, oxygen also promotes the formation of dioxymethylene from formaldehyde, which then easily decays to formate. Formate is the most stable reaction intermediate in the total oxidation. Thus the formate decomposition represents the rate-limiting step in the total oxidation of methanol on copper.

Adsorption of atoms on cu surfaces: a density functional theory study

The chemisorption of atoms (H, N, S, O, and C) on Cu surfaces has been systematically studied by the density functional theory generalized gradient approximation method with the slab model. Our calculated results indicate that the orders of the adsorption energy are H < N < S < O < C on Cu(111) and H < N < O < S < C on Cu(110) and Cu(100). Furthermore, the adsorption energies of the given atoms on Cu(100) are larger than those on Cu(111) and Cu(110). The preferred adsorption sites are a 3-fold hollow site on Cu(111) and a 4-fold hollow site on Cu(100), but the preferred adsorption sites on Cu(110) are different for different adatoms. The energy, as well as the geometry, is in good agreement with the experimental and other theoretical data. In addition, this study focuses on the electronic and geometric properties of the metal-atom (M-A) bond to explain the difference in adsorption energies among adatoms. A detailed investigation of the density of states curves explains the nature of the most stable site. Finally, we test the effect of the coverage and find that the surface coverage has no influence on the preferred adsorption sites of the given adatoms on Cu(110) with the exception of hydrogen and oxygen, but has much influence on the value of the adsorption energy.

Effect of pre-covered oxygen on the dehydrogenation reactions over copper surface: a density functional theory study

Dehydrogenation of five species including CH3OH, CH3O, H2COO, NH3, and H2O over clean and oxygen-modified copper surfaces has been investigated by the first-principle density functional calculations within the generalized gradient approximation. The reaction enthalpies and the activation energies have been calculated for 10 elementary steps corresponding to the direct and oxygen-assisted cleavage of X-H bonds (X = O, N, C). The DFT-GGA results showed that the pre-adsorbed oxygen always facilitates the dehydrogenation reaction by decreasing the reaction enthalpies and the activation energies. The obtained results are in general agreement with experimental observations.

Ab Initio Molecular Dynamics Study of Methanol Adsorption on Copper Clusters

The preferential structures of small copper clusters Cun (n=2-9) and the adsorption of methanol molecules on these clusters are examined with first principles, molecular dynamics simulations. The results show that the copper clusters undergo systematic changes in bond length and bond order associated with altering their preferential structures from one-dimensional structures, to two-dimensional and three-dimensional structures. The results also indicate that low coordination number sites on the copper clusters are both the most favorable for methanol adsorption and have the greatest localization of electronic charge. The simulations predict that charge transfer between the neutral copper clusters and the incident methanol molecules is a key process by which adsorption is stabilized. Importantly, the changes in the dimensionality of the copper clusters do not significantly influence methanol adsorption.

Comparative Theoretical Study of Formaldehyde Decomposition on PdZn, Cu, and Pd Surfaces

Methanol steam reforming, catalyzed by Pd/ZnO (PdZn alloy), is a potential source of hydrogen for on-board fuel cells. CO has been reported to be a minor side product of methanol decomposition that occurs in parallel to methanol steam reforming on PdZn catalysts. However, fuel cells currently used in vehicles are very sensitive to CO poisoning. To contribute to the understanding of pertinent reaction mechanisms, we employed density functional slab model calculations to study the decomposition of formaldehyde, a key intermediate in methanol decomposition and steam reforming reactions, on planar surfaces of Pd, Cu, and PdZn as well as on a stepped surface of PdZn. The calculated activation energies indicate that dehydrogenation of formaldehyde is favorable on Pd(111), but unfavorable on Cu(111) and PdZn(111). On the stepped PdZn(221) surface, the dehydrogenation process was calculated to be more competitive to formaldehyde desorption than on PdZn(111). Thus, we ascribe the experimentally observed small amount of CO, formed during steam reforming of methanol on the Pd/ZnO catalyst, to occur at metallic Pd species of the catalyst or at defect sites of PdZn alloy.

A Systematic Density Functional Theory Study of the C−N Bond Cleavage of Methylamine on Metals

The C-N bond cleavage for the relative large molecule of methylamine on Cu(1111), Ag(111), Au(111), Ni(111), Rh(111), Pd(111), Pt(111), and Mo(100) has been systemically studied using the DFT-GGA method; the reaction energy changes and the activation energies were obtained. The calculated results show that the activation energy of C-N bond cleavage decreases as the metal element goes up and to the left across the periodic table, which is in general agreement with the experimental observation. Moreover, it was found that the steric effect should be considered for the metals with high activity and small radius such as Ni, which is much different from the case for the small molecule decomposition in which the steric effect may be ignored. The linear relationships between the activation energies and electronic properties (d-band center) are presented. It is expected that such a rule can be used to predict the reactivity of metal for other dissociative adsorption systems.

Transition State for β-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.

Density Functional Theory of Water−Gas Shift Reaction on Molybdenum Carbide

The density functional theory (DFT) of the water-gas shift (WGS) reaction over molybdenum carbide was studied with the aim of understanding the dissociation of H(2)O, the OH group, and CO to determine on what sections of molybdenum carbide CO(2) and H(2) formed and whether they played a role in the reaction. The energy diagram of each elementary step, the reaction of the hydrogen and oxygen atoms with CO, and the transition state for this elementary step were also studied. The IR spectra of the CO adsorption was experimentally analyzed for the identification of several candidates of the CO adsorption modes. The adsorptions of the threefold Mo site (a) with and (b) without the underlying C atom of the second layer have the second and highest adsorption energies of -281.59 and -321.00 kJ/mol, respectively. The IR data showed that the bands at 1626 cm(-1) from the IR experiments are (a) the nearest adsorption of the threefold Mo site with the underlying C atom at the calculated/corrected band of 1621 cm(-1). The calculated/corrected threefold adsorption (b) had the highest adsorption energy but exhibited an IR band at 1147 cm(-1) which was not observed in the experimental data. The C-O bond length increased to 1.49 from 1.36 after the H(2)O adsorption (b), suggesting the dissociation of C-O after the H(2)O coadsorption. The WGS reaction on the beta-Mo(2)C(001) slab carbide was calculated and took place as follows: H(2)O was dissociated into OH and H on the Mo(2)C surface and the OH subsequently dissociated into H and O atoms. CO approached the O atom to form CO(2).

Kinetic Mechanism of Methanol Decomposition on Ni(111) Surface: A Theoretical Study

The decomposition of methanol on the Ni(111) surface has been studied with the pseudopotential method of density functional theory-generalized gradient approximation (DFT-GGA) and with the repeated slab models. The adsorption energies of possible species and the activation energy barriers of the possible elementary reactions involved are obtained in the present work. The major reaction path on Ni surfaces involves the O-H bond breaking in CH(3)OH and the further decomposition of the resulting methoxy species to CO and H via stepwise hydrogen abstractions from CH(3)O. The abstraction of hydrogen from methoxy itself is the rate-limiting step. We also confirm that the C-O and C-H bond-breaking paths, which lead to the formation of surface methyl and hydroxyl and hydroxymethyl and atom hydrogen, respectively, have higher energy barriers. Therefore, the final products are the adsorbed CO and H atom.

Theoretical Study of CO and NO Chemisorption on RhCu(111) Surfaces

A systematic density functional theory study using periodic models is presented concerning the chemisorption of CO and NO on various sites of RhCu(111) surfaces. The properties of the adsorbed molecules on various mono- and bimetallic sites of these alloy surfaces have been obtained and compared to those corresponding to the pure Rh(111) and Cu(111) surfaces. It is shown that that the interaction of small probe molecules such as CO or NO on RhCu alloys is essentially dominated by the atomic nature of the surface active site with little influence of the rest of the metallic system. Moreover, it is suggested that it is possible to control the adsorption site of these molecules by appropriate choice of the surface composition.

Characterization of methoxy adsorption on some transition metals: A first principles density functional theory study

Based on the gradient-density functional theory, calculation results of methoxy adsorption on Au(111), Ag(111), Cu(111), Pt(111), Pd(111), Ni(111), Rh(111), and Fe(100) surfaces are presented, and a consistent picture for some key physical properties determining the reactivity of metals appears. These eight metals belong to two groups: either with filled d electrons (group IB) or with unfilled but more than half filled d electrons (group VIII). The calculated adsorption energies are quite in agreement with the experimental data as well as the previous theoretical calculation results. Importantly, using the analysis of B. Hammer and J. K. Norskov, Nature (London) 376, 232 (1995) and in Chemisorption and Reactivity on Supported Clusters and Thin Films, edited by R. M. Lambert and G. Pacchioni (Kluwer Academic, Dordrecht, 1997), pp. 285-351, the binding energies have selectively been linearly correlated to the d-band center and to the size of the metal d-band orbital overlapping with the adsorbate (coupling matrix element) for these two groups of metals. And by analyzing the nature of the adsorption bonding, the possible reason of this difference is suggested.

CH3O decomposition on PdZn(111), Pd(111), and Cu(111). A theoretical study

Methanol steam re-forming, catalyzed by Pd/ZnO, is a potential hydrogen source for fuel cells, in particular in pollution-free vehicles. To contribute to the understanding of pertinent reaction mechanisms, density functional slab model studies on two competing decomposition pathways of adsorbed methoxide (CH(3)O) have been carried out, namely, dehydrogenation to formaldehyde and C-O bond breaking to methyl. For the (111) surfaces of Pd, Cu, and 1:1 Pd-Zn alloy, adsorption complexes of various reactants, intermediates, transition states, and products relevant for the decomposition processes were computationally characterized. On the surface of Pd-Zn alloy, H and all studied C-bound species were found to prefer sites with a majority of Pd atoms, whereas O-bound congeners tend to be located on sites with a majority of Zn atoms. Compared to Pd(111), the adsorption energy of O-bound species was calculated to be larger on PdZn(111), whereas C-bound moieties were less strongly adsorbed. C-H scission of CH(3)O on various substrates under study was demonstrated to proceed easier than C-O bond breaking. The energy barrier for the dehydrogenation of CH(3)O on PdZn(111) (113 kJ mol(-)(1)) and Cu(111) (112 kJ mol(-)(1)) is about 4 times as high as that on Pd(111), due to the fact that CH(3)O interacts more weakly with Pd than with PdZn and Cu surfaces. Calculated results showed that the decomposition of methoxide to formaldehyde is thermodynamically favored on Pd(111), but it is an endothermic process on PdZn(111) and Cu(111) surfaces.