Mechanistic Study of 1,3-Butadiene Formation in Acetylene Hydrogenation over the Pd-Based Catalysts Using Density Functional Calculations

Green oil, which leads to the deactivation of the catalysts used for the selective hydrogenation of acetylene, has long been observed but its formation mechanism is not fully understood. In this work, the formation of 1,3-butadiene, known to be the precursor of green oil, on both Pd(111) and Pd(211) surfaces is examined using density functional theory calculations. The pathways containing C-2 + C-2 coupling reactions as well as the corresponding hydrogenation reactions are studied in detail. Three pathways for 1,3-butadiene production, namely coupling plus hydrogenation and further hydrogenation, hydrogenation plus coupling plus hydrogenation, and a two step hydrogenation followed by coupling, are determined. By comparing the effective barriers, we identify the favored pathway on both surfaces. A general understanding toward the deactivation process of the industrial catalysts is also provided. In addition, the effects of the formation of subsurface carbon atoms as well as the Ag alloying on the 1,3-butadiene formation on Pd-based catalysts are also investigated and compared with experimental results.

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.

Significance of β-dehydrogenation in ethanol electro-oxidation on platinum doped with Ru, Rh, Pd, Os and Ir

In the exploration of highly efficient direct ethanol fuel cells (DEFCs), how to promote the CO₂ selectivity is a key issue which remains to be solved. Some advances have been made, for example, using bimetallic electrocatalysts, Rh has been found to be an efficient additive to platinum to obtain high CO₂ selectivity experimentally. In this work, the mechanism of ethanol electrooxidation is investigated using first principles method. It is found that CH₃CHOH* is the key intermediate during ethanol electrooxidation and the activity of β-dehydrogenation is the rate determining factor that affects the completeness of ethanol oxidation. In addition, a series of transition metals (Ru, Rh, Pd, Os and Ir) are alloyed on the top layer of Pt(111) in order to analyze their effects. The elementary steps, α-, β-C-H bond and C-C bond dissociations are calculated on these bimetallic M/Pt(111) surfaces and the formation potential of OH* from water dissociation is also calculated. We find that the active metals increase the activity of β-dehydrogenation but lower the OH* formation potential resulting in the active site being blocked. By considering both β-dehydrogenation and OH* formation, Ru, Os and Ir are identified to be unsuitable for the promotion of CO₂ selectivity and only Rh is able to increase the selectivity of CO₂ in DEFCs.

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. 

Ethanol Steam Reforming on Rh Catalysts: Theoretical and Experimental Understanding

Through combined theoretical and experimental efforts, the reaction mechanism of ethanol steam reforming on Rh catalysts was studied. The results suggest that acetaldehyde (CH3CHO) is an important reaction intermediate in the reaction on nanosized Rh catalyst. Our theoretical work suggests that the H-bond effect significantly modifies the ethanol decomposition pathway. The possible reaction pathway on Rh (211) surface is suggested as CH3CH2OH -> CH3CH2O -> CH3CHO -> CH3CO -> CH3 + CO -> CH2 + CO -> CH + CO -> C + CO, followed by the water gas shift reaction to yield H-2 and CO2. In addition, we found that the water-gas shift reaction, not the ethanol decomposition, is the bottleneck for the overall ethanol steam reforming process. The CO + OH association is considered the key step, with a sizable energy barrier of 1.31 eV. The present work first discusses the mechanisms and the water effect in ethanol steam reforming reactions on Rh catalyst from both theoretical and experimental standpoints, which may shed light on designing improved catalysts.

Evidence To Challenge the Universality of the Horiuti-Polanyi Mechanism for Hydrogenation in Heterogeneous Catalysis:Origin and Trend of the Preference of a Non-Horiuti-Polanyi Mechanism

The Horiuti-Polanyi mechanism has been considered to be universal for explaining the mechanisms of hydrogenation reactions in heterogeneous catalysis for several decades. In this work, we examine this mechanism for the hydrogenation of acrolein, the simplest alpha,beta-unsaturated aldehyde, in gold-based systems as well as some other metals using extensive first-principles calculations. It is found that a non-Horiuti-Polanyi mechanism is favored in some cases. Furthermore, the physical origin and trend of this mechanism are revealed and discussed regarding the geometrical and electronic effects, which will have a significant influence on current understandings on heterogeneous catalytic hydrogenation reactions and the future catalyst design for these reactions.

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.

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.

Molecular electrostatic properties of ions in an ionic liquid

We have analysed the electronic wave functions from an ab initio simulation of the ionic liquid (room temperature molten salt) dimethyl imidazolium chloride ([dmim][Cl] or [C1mim][Cl]) using localized Wannier orbitals. This allows us to assign electron density to individual ions. The probability distributions of the ionic dipole moments for an isolated ion and for ions in solution are compared. The liquid environment is found to polarize the cation by about 0.7 D and to increase the amplitude of the fluctuations in the dipole moments of both cation and anion. The relative changes in nuclear and electronic contributions are shown. The implications for classical force fields are discussed.

Suggested Binding Mechanism of the HIV-gp120 to its CD4 Receptor

The molecular recognition and attachment of the CD4 molecule and the HIV envelope glycoprotein (gp120) might be described as a consecutive three-step molecular recognition process. 1. (a) Long range interaction: electrostatic pre-orientation, 2. (b) short range interaction: electronic attachment followed by a ‘Locking-in’ (via aromatic ring orientation) and 3. (c) internal interaction (induced fit): conformational readjustment of the protein molecules. On the basis of the preliminary investigations (X-ray structures of CD4 and biological studies of CD4 and gp120 point mutants) we described a computational model. This approach consists of empirical calculations as well as ab initio level of quantum chemistry. The conformational analysis of the wild type and mutant CD4 molecules was supported by molecular mechanics and dynamics (Amber force field). The latter analysis involves the application of a novel method, the Amino Acid Conformation Assignment of Proteins (ACAP) software, developed for the notation of secondary protein structures. According to the cardinal role of the electrostatic factors during this interaction, several ab initio investigations were performed for better understanding of the recognition process on submolecular level. Using the above mentioned computational model, we could interpret the basic behaviours and predict some additional features of CD4-gp120 interaction, in spite of the missing gp120 X-ray structure.

Molecular Modeling on Inhibitor Complexes and Active-Site Dynamics of Cytochrome P450 C17, a Target for Prostate Cancer Therapy

A molecular model for the P450 enzyme cytochrome P450 C17 (CYP17) is presented based on sequence alignments of multiple template structures and homology modeling. This enzyme plays a central role in the biosynthesis of testosterone and is emerging as a major target in prostate cancer, with the recently developed inhibitor abiraterone currently in advanced clinical trials. The model is described in detail, together with its validation, by providing structural explanations to available site-directed mutagenesis data. The CYP17 molecule in this model is in the form of a triangular prism, with an edge of similar to 55 angstrom and a thickness of similar to 37 angstrom. It is predominantly helical, comprising 13 alpha helices interspersed by six 3(10) helices and 11 beta-sheets. Multinanosecond molecular dynamics simulations in explicit solvent have been carried out, and principal components analysis has been used to reveal the details of dynamics around the active site. Coarse-grained methods have also been used to verify low-frequency motions, which have been correlated with active-site gating. The work also describes the results of docking synthetic inhibitors, including the drug abiraterone and the natural substrate pregnenolone, in the CYP17 active site together with molecular dynamics simulations on the complexes. (C) 2010 Elsevier Ltd. All rights reserved.

Pinning down the solid-state polymorphism of the ionic liquid [bmim][PF6]

The solid-state polymorphism of the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate, [bmim][PF6], has been investigated via low-temperature and high-pressure crystallisation experiments. The samples have been characterised by single-crystal X-ray diffraction, optical microscopy and Raman spectroscopy. The solid-state phase behaviour of the compound is confirmed and clarified with respect to previous phase diagrams. The structures of the previously reported gamma-form, which essentially exhibits a G'T cation conformation, as well as those of the elusive beta- and alpha-forms, are reported. Crystals of the beta-phase are twinned and the structure is heavily disordered; the cation conformation in this form is predominantly TT, though significant contributions from other less frequently encountered conformers are also observed at low temperature and high pressure. The cation conformation in the alpha-form is GT; the presence of the G'T conformer at 193 K in this phase can be eliminated on cooling to 100 K. Whilst X-ray structural data are overall in good agreement with previous interpretations based on Raman and NMR studies, they also reveal a more subtle interplay of intermolecular interactions, which give rise to a wider range of conformers than previously considered.

Positron scattering and annihilation on noble-gas atoms

Positron scattering and annihilation on noble-gas atoms is studied ab initio using many-body theory methods for positron energies below the positronium formation threshold. We show that in this energy range, the many-body theory yields accurate numerical results and provides a near-complete understanding of the positron–noble-gas atom system. It accounts for positron-atom and electron-positron correlations, including the polarization of the atom by the positron and the nonperturbative effect of virtual positronium formation. These correlations have a large influence on the scattering dynamics and result in a strong enhancement of the annihilation rates compared to the independent-particle mean-field description. Computed elastic scattering cross sections are found to be in good agreement with recent experimental results and Kohn variational and convergent close-coupling calculations. The calculated values of the annihilation rate parameter Zeff (effective number of electrons participating in annihilation) rise steeply along the sequence of noble-gas atoms due to the increasing strength of the correlation effects, and agree well with experimental data.