Search

Microkinetic modeling has drawn increasing attention for quantitatively analyzing catalytic networks in recent decades, in which the speed and stability of the solver play a crucial role. However, for the multi-step complex systems with a wide variation of rate constants, the often encountered stiff problem leads to the low success rate and high computational cost in the numerical solution. Here, we report a new efficient sensitivity-supervised interlock algorithm (SSIA), which enables us to solve the steady state of heterogeneous catalytic systems in the microkinetic modeling with a 100% success rate. In SSIA, we introduce the coverage sensitivity of surface intermediates to monitor the low-precision time-integration of ordinary differential equations, through which a quasi-steady-state is located. Further optimized by the high-precision damped Newton's method, this quasi-steady-state can converge with a low computational cost. Besides, to simulate the large differences (usually by orders of magnitude) among the practical coverages of different intermediates, we propose the initial coverages in SSIA to be generated in exponential space, which allows a larger and more realistic search scope. On examining three representative catalytic models, we demonstrate that SSIA is superior in both speed and robustness compared with its traditional counterparts. This efficient algorithm can be promisingly applied in existing microkinetic solvers to achieve large-scale modeling of stiff catalytic networks.

Heterogeneous catalysis on supported and nonsupported nanoparticles is of fundamental importance in the energy and chemical conversion industries. Rather than laboratory analysis, first-principles calculations give us an atomic-level understanding of the structure and reactivity of nanoparticles and supports, greatly reducing the efforts of screening and design. However, unlike catalysis on low index single crystalline surfaces, nanoparticle catalysis relies on the tandem properties of a support material as well as the metal cluster itself, often with charge transfer processes being of key importance. In this perspective, we examine current state-of-the-art quantum-chemical research for the modeling of reactions that utilize small transition metal clusters on metal oxide supports. This should provide readers with useful insights when dealing with chemical reactions on such systems, before discussing the possibilities and challenges in the field.

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 (N2+3H2→2NH3). 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...

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 (∼200 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 barr...

The metabolism of hydrogen (H2↔2H++2e−) constitutes a central process in the global biological energy cycle. Among all the enzymes that can mediate this process, Fe-only hydrogenases are unique in their particular high reactivity. Recently, some important progresses have been achieved. Following our previous paper [Z.-P. Liu and P. Hu, J. Am. Chem. Soc. 124, 5175 (2002)] that characterizes the individual redox state of the active site of Fe-only hydrogenase, in this work we have determined the feasible reaction pathways and energetics for the H2 metabolism on the active site of Fe-only hydrogenases, using density functional theory. We show that H2 metabolism possesses very low reaction barriers and a proximal base from a nearby protein plays an important role in H2 metabolism.

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 Ea, in contrast to CO oxidation in which the activation of the O(ads) is much greater than CO(ads) activation.

Methane activation is a crucial step in the conversion of methane to valuable oxygenated products. In heterogeneous catalysis, however, methane activation often leads to complete dissociation: If a catalyst can activate the first C–H bond in CH4, it can often break the remaining C–H bonds. In this study, using density functional theory, we illustrate that single C–H bond activation in CH4 is possible. We choose a model system which consists of isolated Pt atoms on a MoO3(010) surface. We find that the Pt atoms on this surface can readily activate the first C–H bond in methane. The reaction barrier of only 0.3 eV obtained in this study is significantly lower than that on a Pt(111) surface. We also find, in contrast to the processes on pure metal surfaces, that the further dehydrogenation of methyl (CH3) is very energetically unfavorable on the MoO3-supported Pt catalyst.

Increasing emissions of carbon dioxide and the depletion of finite fossil fuels have led to many concerns about global warming and energy crises. Consequently, a sustainable and alternative method, photocatalytic CO2 reduction to chemical fuels has received considerable attention. This perspective highlights recent knowledge and the main challenges in CO2 photoreduction primarily from the theoretical field. The fundamental understanding of CO2 adsorption and reaction mechanism at an atomic level is fully addressed, and the relevant effects dominating the process of CO2 photoreduction are also elucidated. Moreover, recent development of photocatalysts including structural modification is presented, which greatly affects the efficiency and selectivity of CO2 conversion. Finally, the possibilities and challenges in this field are discussed. © 2017 Author(s).

Density-functional theory has been used to investigate the chemisorption of S, SH, and H2S as well as the coadsorption of S and H and SH and H on Pt(111). In addition reaction pathways and energy profiles for the conversion of adsorbed S and H into gas-phase H2S have been determined. It has been found that S, SH, and H2S bind preferentially at face-centered-cubic (fcc), bridge, and top sites, respectively. Both the S+H and SH+H reactions have high barriers (∼1 eV) and high exothermicities (∼1 eV). This reveals that adsorbed H2S and SH are highly unstable adsorbates on Pt(111) and that adsorbed S (and H) is the most stable SHX (X=0,1,2) intermediate on Pt(111)

There is some dispute as to whether methanol decomposition occurs by O–H bond scission or C–O bond scission. By carrying out density functional theory calculations, we investigate both scenario of the reaction pathways of methanol decomposition on a Pd(111) surface. It is shown that the O–H bond scission pathway is much more energetically favorable than the C–O bond scission pathway. The high reaction barrier in the latter case is found to be due to the poor bonding abilities of CH3 and OH with the surface at the reaction sites.

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 Oa) 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...

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.

Ruthenium is one of the poorest catalysts for CO oxidation under normal conditions (low or medium O coverage and normal temperature). However, a recent study [Science 285, 1042 (1999)] reveals that, under femtosecond laser irradiation, CO2 can be formed on the Ru surface, and the reaction follows an electron-mediated mechanism. We carried out density functional theory calculations to investigate CO oxidation via an electron-mediated mechanism on Ru(0001). By comparison to the reaction under normal conditions, following features emerge in the electron-mediated mechanism: (i) more reaction channels are open; (ii) the reaction barrier is significantly lowered. The physical origins for these novel features have been analyzed.

Density functional theory (DFT) has been used to determine reaction pathways for several reactions taking place on Pt(111) and Cu(111) surfaces. On Pt(111), the reactions of C+O and C+N were studied, and on Cu(111) we investigated the reaction of C+H. The structures of the transition states accessed in each reaction are similar. An equivalent distance separates the reactants with the first located at a three-fold hollow site and the second close to a bridge site. Previous DFT studies have, in fact, often identified transition states of this type and in every case it is the reactant with the weaker chemisorption energy that is located close to the bridge site. An explanation as to why this is so is provided.

To investigate the softening of CH vibrational frequencies and their implications for dehydrogenation of adsorbed hydrocarbons, an issue of scientific and technological importance, density functional theory calculations have been performed on the chemisorption and dehydrogenation of CH3 on Cu(111) and Pt(111) surfaces. By comparing these results with those of Ni(111) we find that the CH bonds of the adsorbate, when close enough, interact with metal atoms of the surface. It is this interaction and its associated lengthening and weakening of CH bonds that is the physical origin of mode softening. We rule out the possibility of a relationship between the mere presence of mode softening and dehydrogenation. We do show, however, that there is a clear relationship between the extent to which a surface can induce mode softening and the activation energy to dehydrogenation. In addition, periodic trends concerning the extent of mode softening are reproduced.

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.

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 s...

The ground state potential energy surface for CO chemisorption across Pd{110} has been calculated using density functional theory with gradient corrections at monolayer coverage. The most stable site corresponds well with the experimental adsorption heat, and it is found that the strength of binding to sites is in the following order: pseudo-short-bridge>atop>long-bridge>hollow. Pathways and transition states for CO surface diffusion, involving a correlation between translation and orientation, are proposed and discussed.

In heterogeneous catalysis, the two main reaction mechanisms which have been proposed are the Langmuir–Hinshelwood and the Eley–Rideal. For the vast majority of surface catalytic reactions, it has been accepted that the Langmuir–Hinshelwood mechanism is preferred. In this study, we investigate catalytic CO oxidation on Pt(111). It is found that reaction barriers for Langmuir–Hinshelwood mechanisms actually tend to be higher than those for Eley–Rideal ones. An explanation is presented as to why it is still more probable for the reaction to proceed via the Langmuir–Hinshelwood mechanism, despite its higher reaction barrier.

Catalytic ammonia synthesis is believed to proceed via dissociation of N2 and H2 with subsequent stepwise addition reactions from an adsorbed nitrogen atom to NH3. The first step, N2 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 N2 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 N2 dissociation on the same stepped surface, which indicat...

In the catalytic hydrogenation of hydrocarbons, subsurface hydrogen is known experimentally to be much more reactive than surface hydrogen. We use density functional theory to identify low-energy pathways for the hydrogenation of methyl adsorbed on Ni(111) by surface and subsurface hydrogen. The metastability of subsurface hydrogen with respect to chemisorbed hydrogen is mainly responsible for the low activation barrier for subsurface reactions.

A long-standing puzzle in heterogeneous catalysis, CO oxidation on a Ru(0001) surface under high oxygen pressure, was clarified recently by Over et al. [Science 287, 1474 (2000)], who obtained experimental evidence that the active phase is RuO2(110) when the reaction is operated at an elevated temperature. To find microscopic details of the reaction and to understand the high reactivity of RuO2(110), we performed density functional theory calculations for CO oxidation on RuO2(110). The transition state is located and the reaction barrier is determined to be 1.15 eV, being considerably lower than that on Ru(0001). The high reactivity of RuO2(110) for CO oxidation is analyzed.

The crucial roles of the coverage of surface free sites in determining catalytic activity trend are quantitatively addressed with the help of density functional theory and microkinetics. First, by analyzing activity trends of NO oxidation catalyzed by Ru, Rh, Pd, Os, Ir, and Pt surfaces with full kinetic considerations, we identify that the activity trend is in general determined by the competition between the reaction barrier and the coverage of surface free sites. Second, since the dissociation of many important molecules, such as the dissociation of N(2), O(2), and CO, follows the same Bronsted-Evans-Polanyi relationship, the coverage of surface free sites is usually a decisive term that affects the overall activity. Third, an equation is derived for the coverage of surface free sites and it is found that the coverage of surface free sites contains not only all the key thermodynamic parameters but also all the kinetic properties in the catalytic system.

The carbazole moiety is a component of many important pharmaceuticals including anticancer and anti-HIV agents and is commonly utilized in the production of modern polymeric materials with novel photophysical and electronic properties. Simple carbazoles are generally produced via the aromatization of the respective tetrahydrocarbazole (THCZ). In this work, density functional theory calculations are used to model the reaction pathway of tetrahydrocarbazole aromatization over Pd(111). The geometry of each of the intermediate surface species has been determined and how each structure interacts with the metal surface addressed. The reaction energies and barriers of each of the elementary surface reactions have also been calculated, and a detailed analysis of the energetic trends performed. Our calculations have shown that the surface intermediates remain fixed to the surface via the aromatic ring in a manner similar to that of THCZ. Moreover, the aliphatic ring becomes progressively more planer with the dissociation of each subsequent hydrogen atom. Analysis of the reaction energy profile has revealed that the trend in reaction barriers is determined by the two factors: (i) the strength of the dissociating ring-H bond and (ii) the subsequent gain in energy due to the geometric relaxation of the aliphatic ring.

The C-H activation on metal oxides is a fundamental process in chemistry. In this paper, we report a density functional theory study on the process of the C-H activation of CH(4) on Pd(111), Pt(111), Ru(0001), Tc(0001), Cu(111), PdO(001), PdO(110), and PdO(100). A linear relationship between the C-H activation barrier and the chemisorption in the dissociation final state on the metal surfaces is obtained, which is consistent with the work in the literature. However, the relationship is poor on the metal oxide surfaces. Instead, a strong linear correlation between the barrier and the lattice O-H bond strength is found on the oxides. The new linear relationship is analyzed and the physical origin is identified.

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 NH(3) 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.

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.

Hydrogenation is an important process in the Fischer-Tropsch synthesis. In this work, all the elementary steps of the hydrogenation from C to CH4 are studied on both flat and stepped Co(0001) using density functional theory (DFT). We found that (i) CH3 hydrogenation (CH3+H--CH4) is the most difficult one among all the elementary reactions on both surfaces, possessing barriers of around 1.0 eV; (ii) the other elementary reactions have the barriers below 0.9 eV on the flat and stepped surfaces; (iii) CH2 is the least stable species among all the CHx(x=1-3) species on both surfaces; and (iv) surface restructuring may have little effect on the CHx(x=0-3) hydrogenation. The barriers of each elementary step on both flat and stepped surfaces are similar and energy profiles are also similar. The reason as to why CHx hydrogenation is not structure-sensitive is also discussed.

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.

Density functional theory has been used to study the adsorption of hydroxyl at low and high coverages and also to investigate the nature of the intermediate in the H2O formation reaction on Pt(111). At low coverages [1/9 of a monolayer (ML) to 1/3 ML] OH binds preferentially at bridge and top sites with a chemisorption energy of ∼2.25 eV. At high coverages (1/2 ML to 1 ML) H bonding between adjacent hydroxyls causes: (i) an enhancement in OH chemisorption energy by about 15%; (ii) a strong preference for OH adsorption at top sites; and (iii) the formation of OH networks. The activation energy for the diffusion of isolated OH groups along close packed rows of Pt atoms is 0.1 eV. This low barrier coupled with H bonding between neighboring OH groups indicates that hydroxyls are susceptible to island formation at low coverages. Pure OH as well as coadsorbed OH and H can be ruled out as the observed low temperature intermediate in the water formation reaction. Instead we suggest that the intermediate consists ...