Influence of surface structures, subsurface carbon and hydrogen, and surface alloying on the activity and selectivity of acetylene hydrogenation on Pd surfaces: A density functional theory study

The selective hydrogenation of acetylene to ethylene on several Pd surfaces (Pd(111), Pd(100), Pd(211), and Pd(211)-defect) and Pd surfaces with subsurface species (carbon and hydrogen) as well as a number of Pd-based alloys (Pd-M/Pd(111) and Pd-M/Pd(211) (M = Cu, Ag and Au)) are investigated using density functional theory calculations to understand both the acetylene hydrogenation activity and the selectivity of ethylene formation. All the hydrogenation barriers are calculated, and the reaction rates on these surfaces are obtained using a two-step model. Pd(211) is found to have the highest activity for acetylene hydrogenation while Pd(100) gives rise to the lowest activity. In addition, more open surfaces result in over-hydrogenation to form ethane, while the close-packed surface (Pd(111)) is the most selective. However, we also find that the presence of subsurface carbon and hydrogen significantly changes the reactivity and selectivity of acetylene toward hydrogenation on Pd surfaces. On forming surface alloys of Pd with Cu, Ag and Au, the selectivity for ethylene is also found to be changed. A new energy decomposition method is used to quantitatively analyze the factors in determining the changes in selectivity. These surface modifiers are found to block low coordination unselective sites, leading to a decreased ethane production. (C) 2013 The Authors. Published by Elsevier Inc. All rights reserved.

Selective Hydrogenation of Acetylene over Pd-Boron Catalysts: A Density Functional Theory Study

Boron-modified Pd catalysts have shown excellent performance for the selective hydrogenation of alkynes experimentally. In the current work, we investigated the hydrogenation of acetylene on boron-modified Pd(111) and Pd(211) surfaces, utilizing density functional theory calculations. The activity of acetylene hydrogenation has been studied by estimating the effective barrier of the whole process. The selectivity of ethylene formation is investigated from a comparison between the desorption and the hydrogenation of ethylene as well as comparison between the ethylene and the 1,3-butadiene formation. Formation of subsurface carbon and hydrogen on both boron-modified Pd(111) and Pd(211) surfaces has also been evaluated, since these have been reported to affect both the activity and the selectivity of acetylene hydrogenation to produce ethylene on Pd surfaces. Our results provide some important insights into the Pd B catalysts for selective hydrogenation of acetylene and also for more complex hydrogenation systems, such as stereoselective hydrogenation of longer chain alkynes and selective hydrogenation of vegetable oil.

Density functional theory study of mixed-phase TiO2: heterostructures and electronic properties

In this work, density functional theory calculations have been performed to study the geometric, electronic, and energetic properties of two-phase TiO2 composites built by joining two single-phase TiO2 slabs, aiming at verifying possible improvement of the photo-activities of the composites through phase separation of excitons. We find that such desired electronic properties can be determined by several factors. When both the HOMO and LUMO levels of one of the two single-phase TiO2 slabs are higher than the corresponding ones of the other, the composite may have native electronic structures with phase-separated HOMO-LUMO states, especially when the two slabs exhibit highly matched surface lattices. For those pairs of TiO2 slabs with the HOMO and LUMO levels of one phase being within the range of those of the other, though the energetically favored composite give HOMO-LUMO states within one phase, one may still be able to separate them and move the HOMO state to the interface region by destabilizing the interactions between the two slabs.

Activity and coke formation of nickel and nickel carbide in dry reforming: A deactivation scheme from density functional theory

Dry reforming is a promising reaction to utilise the greenhouse gases CO2 and CH4. Nickel-based catalysts are the most popular catalysts for the reaction, and the coke formation on the catalysts is the main obstacle to the commercialisation of dry reforming. In this study, the whole reaction network of dry reformation on both flat and stepped nickel catalysts (Ni(111) and Ni(211)) as well as nickel carbide (flat: Ni3C(001); stepped: Ni3C(111)) is investigated using density functional theory calculations. The overall reaction energy profiles in the free energy landscape are obtained, and kinetic analyses are utilised to evaluate the activity of the four surfaces. By careful examination of our results, we find the following regarding the activity: (i) flat surfaces are more active than stepped surfaces for the dry reforming and (ii) metallic nickel catalysts are more active than those of nickel carbide, and therefore, the phase transformation from nickel to nickel carbide will reduce the activity. With respect to the coke formation, the following is found: (i) the coke formation probability can be measured by the rate ratio of CH oxidation pathway to C oxidation pathway (r(CH)/r(C)) and the barrier of CO dissociation, (ii) on Ni(111), the coke is unlikely to form, and (iii) the coke formations on the stepped surfaces of both nickel and nickel carbide can readily occur. A deactivation scheme, using which experimental results can be rationalised, is proposed. 

A density functional theory study of hydrogen dissociation and diffusion at the perimeter sites of Au/TiO2

Reactivity of supported gold catalysts is a hot topic in catalysis for many years. This communication reports an investigation on the dissociation of molecular hydrogen at the perimeter sites of Au/TiO2 and the spillover of hydrogen atoms from the gold to the support using density functional theory calculations. It is found that the heterolytic dissociation is favoured in comparison with homolytic dissociation of molecular hydrogen at the perimeter sites. However, the surface oxygen of the rutile TiO2(110) surface at these sites can be readily passivated by the formed OH, suggesting that further dissociation of molecular hydrogen may occur at pure gold sites.

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.

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.

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 on the interaction between chemisorbed CO and S on Rh(111)

Density functional theory calculations are carried out for Rh(111)-p(2 x 2)-CO, Rh(111)-p(2 x 2)-S, Rh(111)-p(2 x 2)-(S + CO), Rh(111)-p(3 x 3)-CO, Rh(111)-p(3 x 3)-S and Rh(111)-p(3 x 3)-(S + CO), aiming to shed some light on the S poisoning effect. Geometrical structures of these systems are optimized and chemisorption energies are determined. The presence of S does not significantly influence the geometrical structure and chemisorption energy of CO and vice versa, which strongly suggests that the interaction between CO and S on the Rh(111) surface is mainly short-range in nature. The long range electronic effect for the dramatic attenuation of the CO methanation activity by sulfur is likely to be incorrect. It is suggested that an ensemble effect may be dominant in the catalytic deactivation. (C) 1999 Elsevier Science B.V. All rights reserved.

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.

A Density Functional Theory Study on the Active Center of Fe-Only Hydrogenase: Characterization and Electronic Structure of the Redox States

We have carried out extensive density functional theory (DFT) calculations for possible redox states of the active center in Fe-only hydrogenases. The active center is modeled by [(H(CH(3))S)(CO)(CN(-))Fe(p)(mu-DTN)(mu-CO)Fe(d)(CO)(CN(-))(L)](z) (z is the net charge in the complex; Fe(p)= the proximal Fe, Fe(d) = the distal Fe, DTN = (-SCH(2)NHCH(2)S-), L is the ligand that bonds with the Fed at the trans position to the bridging CO). Structures of possible redox states are optimized, and CO stretching frequencies are calculated. By a detailed comparison of all the calculated structures and the vibrational frequencies with the available experimental data, we find that (i) the fully oxidized, inactive state is an Fe(II)-Fe(II) state with a hydroxyl (OH(-)) group bonded at the Fe(d), (ii) the oxidized, active state is an Fe(II)-Fe(l) complex which is consistent with the assignment of Cao and Hall (J. Am. Chem. Soc. 2001, 123, 3734), and (iii) the fully reduced state is a mixture with the major component being a protonated Fe(l)-Fe(l) complex and the other component being its self-arranged form, Fe(II)-Fe(II) hydride, Our calculations also show that the exogenous CO can strongly bond with the Fe(II)-Fe(l) species, but cannot bond with the Fe(l)-Fe(l) complex. This result is consistent with experiments that CO tends to inhibit the oxidized, active state, but not the fully reduced state. The electronic structures of all the redox states have been analyzed. It is found that a frontier orbital which is a mixing state between the e(g) of Fe and the 2pi of the bridging CO plays a key role concerning the reactivity of Fe-only hydrogenases: (1) it is unoccupied in the fully oxidized, inactive state, half-occupied in the oxidized, active state, and fully occupied in the fully reduced state; (ii) the e(g)-2pi orbital is a bonding state, and this is the key reason for stability of the low oxidation states, such as Fe(l)-Fe(l) complexes; and (iii) in the e(g)-2pi orbital more charge accumulates between the bridging CO and the Fe(d) than between the bridging CO and the Fe(p), and the occupation increase in this orbital will enhance the bonding between the bridging CO and the Fe(d), leading to the bridging-CO shift toward the Fe(d).

Insight into the solvent effect: A density functional theory study of cisplatin hydrolysis

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. (c) 2006 American Institute of Physics.

A density functional theory study of the a-olefin selectivity in Fischer-Tropsch synthesis

We studied the alpha-olefin selectivity in Fischer-Tropsch (FT) synthesis using density functional theory (131717) calculations. We calculated the relevant elementary steps from C-2 to C-6 species. Our results showed that the barriers of hydrogenation and dehydrogenation reactions were constant with different chain lengths, and the chemisorption energies of alpha-olefins from DFT calculations also were very similar, except for C-2 species. A simple expression of the paraffin/olefin ratio was obtained based on a kinetic model. Combining the expression of the paraffin/olefin ratio and our calculation results, experimental findings are satisfactorily explained. We found that the physical origin of the chain length dependence of paraffin/olefin ratio is the chain length dependence of both the van der Waals interaction between adsorbed alpha-olefins and metal surfaces and the entropy difference between adsorbed and gaseous alpha-olefins, and that the greater chemisorption energy of ethylene is the main reason for the abnormal ethane/ethylene ratio. (c) 2008 Elsevier Inc. All rights reserved.

A density functional theory study of the surface relaxation and reactivity of Cu2O(100)

Density functional theory has been used to investigate the surface relaxation of Cu2O(100) and the adsorption of NO. The calculations indicate the formation of surface copper dimers on relaxation coupled with a large contraction of the spacing between the first and second layers. Local density of states for atoms in the top three layers shows that the third layer copper atoms have the greatest change in bonding character. Adsorption energies have been calculated for the N-down and O-down adsorption of NO on the Cu2O(100) surface. These indicate that N-down adsorption is favoured and that in this case NO-lattice oxygen interactions dominate the adsorbate structure. (C) 2000 Elsevier Science B.V. All rights reserved.

Reaction Mechanisms of Crotonaldehyde Hydrogenation on Pt(111): Density Functional Theory and Microkinetic Modeling

The microkinetics based on density function theory (DFT) calculations is utilized to investigate the reaction mechanism of crotonaldehyde hydrogenation on Pt(111) in the free energy landscape. The dominant reaction channel of each hydrogenation product is identified. Each of them begins with the first surface hydrogenation of the carbonyl oxygen of crotonaldehyde on the surface. A new mechanism, 1,4-addition mechanism generating enols (butenol), which readily tautomerize to saturated aldehydes (butanal), is identified as a primary mechanism to yield saturated aldehydes instead of the 3,4-addition via direct hydrogenation of the ethylenic bond. The calculation results also show that the full hydrogenation product, butylalcohol, mainly stems from the deep hydrogenation of surface open-shell dihydrogenation intermediates. It is found that the apparent barriers of the dominant pathways to yield three final products are similar on P(111), which makes it difficult to achieve a high selectivity to the desired crotyl alcohol (COL).