Insight into association reactions on metal surfaces: Density-functional theory studies of hydrogenation reactions on Rh(111)

Hydrogenation reaction, as one of the simplest association reactions on surfaces, is of great importance both scientifically and technologically. They are essential steps in many industrial processes in heterogeneous catalysis, such as ammonia synthesis (N-2+3H(2)-->2NH(3)). Many issues in hydrogenation reactions remain largely elusive. In this work, the NHx (x=0,1,2) hydrogenation reactions (N+H-->NH, NH+H-->NH2 and NH2+H-->NH3) on Rh(111) are used as a model system to study the hydrogenation reactions on metal surfaces in general using density-functional theory. In addition, C and O hydrogenation (C+H-->CH and O+H-->OH) and several oxygenation reactions, i.e., C+O, N+O, O+O reactions, are also calculated in order to provide a further understanding of the barrier of association reactions. The reaction pathways and the barriers of all these reactions are determined and reported. For the C, N, NH, and O hydrogenation reactions, it is found that there is a linear relationship between the barrier and the valency of R (R=C, N, NH, and O). Detailed analyses are carried out to rationalize the barriers of the reactions, which shows that: (i) The interaction energy between two reactants in the transition state plays an important role in determining the trend in the barriers; (ii) there are two major components in the interaction energy: The bonding competition and the direct Pauli repulsion; and (iii) the Pauli repulsion effect is responsible for the linear valency-barrier trend in the C, N, NH, and O hydrogenation reactions. For the NH2+H reaction, which is different from other hydrogenation reactions studied, the energy cost of the NH2 activation from the IS to the TS is the main part of the barrier. The potential energy surface of the NH2 on metal surfaces is thus crucial to the barrier of NH2+H reaction. Three important factors that can affect the barrier of association reactions are generalized: (i) The bonding competition effect; (ii) the local charge densities of the reactants along the reaction direction; and (iii) the potential energy surface of the reactants on the surface. The lowest energy pathway for a surface association reaction should correspond to the one with the best compromise of these three factors. (C) 2003 American Institute of Physics.

The catalytic role of water in CO oxidation

Water, one of the most popular species in our planet, can play a catalytic role in many reactions, including reactions in heterogeneous catalysis. In a recent experimental work, Bergeld, Kasemo, and Chakarov demonstrated that water is able to promote CO oxidation under low temperatures (similar to200 K). In this study, we choose CO oxidation on Pt(111) in the presence of water as a model system to address the catalytic role of water for surface reactions in general using density functional theory. Many elementary steps possibly involved in the CO oxidation on Pt(111) at low temperatures have been investigated. We find the following. First, in the presence of water, the CO oxidation barrier is reduced to 0.33 eV (without water the barrier is 0.80 eV). This barrier reduction is mainly due to the H-bonding between the H in the H2O and the O at the transition state (TS), which stabilizes the TS. Second, CO can readily react with OH with a barrier of 0.44 eV, while COOH dissociation to produce CO2 is not easy (the barrier is 1.02 eV). Third, in the H2O+OH mixed phase, CO can be easily converted into CO2. It occurs through two steps: CO reacts with OH, forming COOH; and COOH transfers the H to a nearby H2O and, at the same time, an H in the H2O transfers to a OH, leading to CO2 formation. The reaction barrier of this process is 0.60 eV under CO coverage of 1/6 ML and 0.33 eV under CO coverage of 1/3 ML. The mechanism of CO oxidation at low temperatures is discussed. On the basis of our calculations, we propose that the water promotion effect can in general be divided into two classes: (i) By H-bonding between the H of H2O and an electron negative species such as the O in the reaction of CO+O+H2O-->CO2+H2O, H2O can stabilize the TS of the reaction and hence reduce the barrier. (ii) H2O first dissociates into H and OH and then OH or H participates directly in the reaction to induce new reaction mechanism with more favorable routes, in which OH or H can act as an intermediate. (C) 2003 American Institute of Physics.

N2O and NO2 formation on Pt(111): A density functional theory study

Catalytic formation of N2O and NO2 were studied employing density functional theory with generalized gradient approximations, in order to investigate the microscopic reaction pathways of these catalytic processes on a Pt(111) surface. Transition states and reaction barriers for the addition of chemisorbed N or chemisorbed O to NO(ads) producing N2O and NO2, respectively, were calculated. The N2O transition state involves bond formation across the hcp hollow site with an associated reaction barrier of 1.78 eV. NO2 formation favors a fcc hollow site transition state with a barrier of 1.52 eV. The mechanisms for both reactions are compared to CO oxidation on the same surface. The activation of the chemisorbed NO and the chemisorbed N or O from the energetically stable initial state to the transition state are both significant contributors to the overall reaction barrier E-a, in contrast to CO oxidation in which the activation of the O-(ads) is much greater than CO(ads) activation. (C) 2002 American Institute of Physics.

Methane transformation to carbon and hydrogen on Pd(100): Pathways and energetics from density functional theory calculations

Density functional theory with gradient corrections has been employed to study the reaction pathways and the reaction energetics for the transformations of CH4 to C and H on a Pd(100) surface. On examination of transition state structures identified in each elementary reaction, a clear relationship between the valencies of the CHx fragments and the locations of the transition states emerges. The higher the valency of the CHx fragment, the higher the coordination number of the CHx with the surface atoms. The calculated reaction energetics are in good agreement with the experiments. In addition, calculation results are also used to illustrate an interesting issue concerning the CH3 stability on Pd surfaces. (C) 2002 American Institute of Physics.

Metal-Ion Catalysis In the Tetrahymena Ribozyme Reaction

LL catalytic RNAs (ribozymes) require or are stimulated by divalent metal ions, but it has been difficult to separate the contribution of these metal ions to formation of the RNA tertiary structure1 from a more direct role in catalysis. The Tetrahymena ribozyme catalyses cleavage of exogenous RNA2,3 or DNA4,5 substrates with an absolute requirement for Mg2+ or Mn2+ (ref. 6). A DNA substrate, in which the bridging 3' oxygen atom at the cleavage site is replaced by sulphur, is cleaved by the ribozyme about 1,000 times more slowly than the corresponding unmodified DNA substrate when Mg2+ is present as the only divalent metal ion. But addition of Mn2+ or Zn2+ to the reaction relieves this negative effect, with the 3' S–P bond being cleaved nearly as fast as the 3' O–P bond. Considering that Mn2+ and Zn2+ coordinate sulphur more strongly than Mg2+ does7,8, these results indicate that the metal ion contributes directly to catalysis by coordination to the 3' oxygen atom in the transition state, presumably stabilizing the developing negative charge on the leaving group. We conclude that the Tetrahymena ribozyme is a metalloenzyme, with mechanistic similarities to several protein enzymes9–12.

The Plasmodium falciparum Malaria M1 Alanyl Aminopeptidase (PfA-M1): Insights of Catalytic Mechanism and Function from MD Simulations

Malaria caused by several species of Plasmodium is major parasitic disease of humans, causing 1-3 million deaths worldwide annually. The widespread resistance of the human parasite to current drug therapies is of major concern making the identification of new drug targets urgent. While the parasite grows and multiplies inside the host erythrocyte it degrades the host cell hemoglobin and utilizes the released amino acids to synthesize its own proteins. The P. falciparum malarial M1 alanyl-aminopeptidase (PfA-M1) is an enzyme involved in the terminal stages of hemoglobin digestion and the generation of an amino acid pool within the parasite. The enzyme has been validated as a potential drug target since inhibitors of the enzyme block parasite growth in vitro and in vivo. In order to gain further understanding of this enzyme, molecular dynamics simulations using data from a recent crystal structure of PfA-M1 were performed. The results elucidate the pentahedral coordination of the catalytic Zn in these metallo-proteases and provide new insights into the roles of this cation and important active site residues in ligand binding and in the hydrolysis of the peptide bond. Based on the data, we propose a two-step catalytic mechanism, in which the conformation of the active site is altered between the Michaelis complex and the transition state. In addition, the simulations identify global changes in the protein in which conformational transitions in the catalytic domain are transmitted at the opening of the N-terminal 8 angstrom-long channel and at the opening of the 30 angstrom-long C-terminal internal chamber that facilitates entry of peptides to the active site and exit of released amino acids. The possible implications of these global changes with regard to enzyme function are discussed.

Improved description of chemical barriers with generalized gradient approximations (GGAs) and meta-GGAs

The performance of exchange and correlation (xc) functionals of the generalized gradient approximation (GGA) type and of the meta-GGA type in the calculation of chemical reactions is related to topological features of the electron density which, in turn, are connected to the orbital structure of chemical bonds within the Kohn-Sham (KS) theory. Seventeen GGA and meta-GGA xc functionals are assessed for 15 hydrogen abstraction reactions and 3 symmetrical S(N)2 reactions. Systems that are problematic for standard GGAs characteristically have enhanced values of the dimensionless gradient argument s(sigma)(2) with local maxima in the bonding region. The origin of this topological feature is the occupation of valence KS orbitals with an antibonding or essentially nonbonding character. The local enhancement of s(sigma)(2) yields too negative exchange-correlation energies with standard GGAs for the transition state of the S(N)2 reaction, which leads to the reduced calculated reaction barriers. The unwarranted localization of the effective xc hole of the standard GGAs, i.e., the nondynamical correlation that is built into them but is spurious in this case, wields its effect by their s(sigma)(2) dependence. Barriers are improved for xc functionals with the exchange functional OPTX as x component, which has a modified dependence on s(sigma)(2). Standard GGAs also underestimate the barriers for the hydrogen abstraction reactions. In this case the barriers are improved by correlation functionals, such as the Laplacian-dependent (LAP3) functional, which has a modified dependence on the Coulomb correlation of the opposite- and like-spin electrons. The best overall performance is established for the combination OLAP3 of OPTX and LAP3.

Probing myo-inositol 1-phosphate synthase with multisubstrate adducts

The synthesis of a series of carbohydrate-nucleotide hybrids, designed to be multisubstrate adducts mimicking myo-inositol 1-phosphate synthase first oxidative transition state, is reported. Their ability to inhibit the synthase has been assessed and results have been rationalised computationally to estimate their likely binding mode.

GHMP Kinases - Structures, Mechanisms and Potential for Therapeutically Relevant Inhibition

The GHMP kinases are a structurally related family of small molecule kinases named after four of its members - galactokinase, homoserine kinase, mevalonate kinase and phosphomevalonate kinase. The group also includes the enzymes N-acetylgalactosamine kinase, arabinose kinase, mevalonate 5-diphosphate decarboxylase, archeal shikimate kinase and 4-(cytidine 5'-diphospho)-2-c-methyl-D-erythritol kinase. In addition the group includes two members not known to be catalytically active, the Caenorhabditis elegans sex-fate determining protein XOL-1 and the Saccharomyces cerevisiae transcriptional activator Gal3p. Two catalytic mechanisms have been proposed for GHMP kinases. The structure of mevalonate kinase suggests that an aspartate residue acts as an active site base, removing a proton from the substrate to facilitate attack on the ? phosphate of MgATP. In contrast, in homoserine kinase there is no potential catalytic base and it is proposed that catalysis is driven by transition state stabilisation. Potential chemotherapeutic interventions against GHMP kinases fall into three main categories: inhibition of galactokinase to assist suffers of galactosemia, inhibition of mevalonate kinase or mevalonate 5-diphosphate decarboxylase to reduce flux through the cholesterol biosynthesis pathway and inhibition of bacterial GHMP kinases for novel anti-microbial therapies. These are in the early stages of development, but the accumulation of structural and mechanistic data will assist future progress.

Linear Optics Simulation of Quantum Non-Markovian Dynamics

The simulation of open quantum dynamics has recently allowed the direct investigation of the features of system-environment interaction and of their consequences on the evolution of a quantum system. Such interaction threatens the quantum properties of the system, spoiling them and causing the phenomenon of decoherence. Sometimes however a coherent exchange of information takes place between system and environment, memory effects arise and the dynamics of the system becomes non-Markovian. Here we report the experimental realisation of a non-Markovian process where system and environment are coupled through a simulated transverse Ising model. By engineering the evolution in a photonic quantum simulator, we demonstrate the role played by system-environment correlations in the emergence of memory effects.

Kinetics of levoglucosan and formaldehyde formation during cellulose pyrolysis process

The mechanisms and kinetics studies of the formation of levoglucosan and formaldehyde from anhydroglucose radical have been carried out theoretically in this paper. The geometries and frequencies of all the stationary points are calculated at the B3LYP/6-31+G(D,P) level based on quantum mechanics, Six elementary reactions are found, and three global reactions are involved. The variational transition-state rate constants for the elementary reactions are calculated within 450-1500 K. The global rate constants for every pathway are evaluated from the sum of the individual elementary reaction rate constants. The first-order Arrhenius expressions for these six elementary reactions and the three pathways are suggested. By comparing with the experimental data, computational methods without tunneling correction give good description for Path1 (the formation of levoglucosan); while methods with tunneling correction (zero-curvature tunneling and small-curvature tunneling correction) give good results for Path2 (the first possibility for the formation of formaldehyde), all the test methods give similar results for Path3 (the second possibility for the formation of formaldehyde), all the modeling results for Path3 are in good agreement with the experimental data, verifying that it is the most possible way for the formation of formaldehyde during cellulose pyrolysis. © 2012 Elsevier Ltd. All rights reserved.

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.

Aza-annulation of enaminones with itaconic anhydride:kinetic preference for exocyclic enamide products

Enaminones react with itaconic anhydride in methylene chloride at room temperature to give exocyclic enamides sis the major products. These can be readily equilibrated to the thermodynamically more stable endocyclic enamides. In substrates where the exocyclic isomer could not be formed only intractable materials were produced from this reaction. An intermediate in this two step process was detected and identified by proton and C-13 NMR spectroscopy. In two cases chiral enaminones were employed and the relative stereochemistry at the new chiral centre in the product was established by a crystal structure of compound 27.

Density Functional Theory Studies of Ethanol Decomposition on Rh(211)

Density functional theory calculations were carried out to examine the mechanism of ethanol decomposition on the Rh(211) surface. We found that there are two possible decomposition pathways: (1) CH(3)CH(2)OH -> CH(3)CHOH -> CH(3)COH -> CH(3)CO -> CH(3) + CO -> CH(2) + CO -> CH + CO -> C + CO and (2) CH(3)CH(2)OH -> CH(3)CHOH -> CH(3)COH -> CH(2)COH -> CHCOH -> CHCO -> CH + CO -> C + CO. Both pathways have a common intermediate of CH(3)COH, and the key step is the formation of CH(3)CHOH species. According to our calculations, the mechanism of ethanol decomposition on Rh(211) is totally different from that on Rh(111): the reaction proceeds via CH(3)COH rather than an oxametallacycle species (-CH(2)CH(2)O- for Rh( 111)), which implies that the decomposition process is structure sensitive. Further analyses on electronic structures revealed that the preference of the initial C(alpha)-H path is mainly due to the significant reduction of d-electron energy in the presence of the transition state (TS) complex, which may stabilize the TS-surface system. The present work first provides a clear picture for ethanol decomposition on stepped Rh(211), which is an important first step to completely understand the more complicated reactions, like ethanol steam reforming and electrooxidation.

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.