Examining the redox and formate mechanisms for water-gas shift reaction on Au/CeO2 using density functional theory

We perform DFT calculations to investigate the redox and formate mechanisms of water-gas-shift (WGS) reaction on Au/CeO2 catalysts. In the redox mechanism, we analyze all the key elementary steps and find that the OH cleavage is the key step. Three possible pathways of OH cleavage are calculated: (1) OHad '' + *'--> H-ad' + O-ad"; (2) H-ad' + OHad '' --> H-2(g) + O-ad '' + *'; and (3) OHad" + OHad '' --> 2O(ad '') + H-2(g) (*': the free adsorption sites on the oxides; ad': adsorption on the metal; ad": adsorption on the oxide, respectively). In the formate mechanism, we identify all the possible pathways for the formation and decomposition of surface formates in the WGS reaction. It is found that there is a shortcoming in the redox and formate mechanisms which is related to surface oxygen reproduction. Four possible pathways for producing surface oxygen are studied, and all the barriers of the four pathways are more than 1 eV. Our results suggest that the processes to reproduce surface oxygen in the reaction circle are not kinetically easy. (C) 2008 Elsevier B.V. All rights reserved.

A density functional theory study of CO and atomic oxygen chemisorption on Pt(111)

Ab initio total energy calculations within a density functional theory framework have been performed for CO and atomic oxygen chemisorbed on the Pt(111) surface. Optimised geometries and chemisorption energies for CO and O on four high-symmetry sites, namely the top, bridge, fee hollow and hcp hollow sites, are presented, the coverage in all cases being 0.25 ML. The differences in CO adsorption energies between these sites are found to be small, suggesting that the potential energy surface for CO diffusion across Pt(111) is relatively flat. The 5 sigma and 2 pi molecular orbitals of CO are found to contribute to bonding with the metal. Some mixing of the 4 sigma and 1 pi molecular orbitals with metal states is also observed. For atomic oxygen, the most stable adsorption site is found to be the fee hollow site, followed in decreasing order of stability by the hcp hollow and bridge sites, with the top site being the least stable. The differences in chemisorption energies between sites for oxygen are larger than in the case of CO, suggesting a higher barrier to diffusion for atomic oxygen. The co-adsorption of CO and O has also been investigated. Calculated chemisorption energies for CO on an O/fcc-precovered surface show that of the available chemisorption sites, the top site at the oxygen atom's next-nearest neighbour surface metal atom is the most stable, with the other four sites calculated bring at least 0.29 eV less stable. The trend of CO site stability in the coadsorption system is explained in terms of a 'bonding competition' model. (C) 2000 Elsevier Science B.V. All rights reserved.

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 of CH2 and H adsorption on Ni(111)

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 sites. The bonding of all adsorbates to the surface is analyzed in detail. (C) 2000 American Institute of Physics. [S0021-9606(00)71213-X].

Methyl chemisorption on Ni(111) and C-H-M multicentre bonding:a density functional theory study

Density functional theory has been used to study the adsorption of CH3 on Ni(111). CH3 is found to adsorb strongly at all four high symmetry sites of the Ni(111) surface. Calculated adsorption energies of CH3 on the different sites are in the following order: hcp approximate to fcc>bridge>top. The bonding and structures of CH3 on the different sites are analysed in detail. An important factor, namely three-centre bonding between carbon, hydrogen and nickel which contributes to the 'soft' C-H vibrational frequency of CH3 on Ni(111), and may determine the preferred chemisorption site, is stressed. (C) 1999 Elsevier Science B.V. All rights reserved.

Density functional theory study of the interaction between CO and on a Pt surface:CO/Pt(111), O/Pt(111), and CO/O/Pt(111)

Ab initio total energy calculations within the Density Functional Theory framework were carried out for Pt(111), Pt(111)-p(2x2)-CO, Pt(111)-p(2x2)-O, and Pt(111)-p(2x2)-(CO+O) to provide an insight into the interaction between CO and O on metal surfaces, an important issue in CO oxidation, and also in promotion and poisoning effects of catalysis. The geometrical structures of these systems were optimized with respect to the total energy, the results of which agree with existing experimental values very well. It is found that (i) the local structures of Pt(111)-p(2x2)-(CO+O), such as the bond lengths of C-O, C-Pt, and O-Pt (chemisorbed O atom with Pt), are almost the same as that in Pt(111)-p(2x2)-CO and Pt(111)-p(2x2)-O, respectively, (ii) the total valence charge density distributions in Pt(111)-p(2x2)-(CO+O) are very similar to that in Pt(111)-p(2x2)-CO, except in the region of the chemisorbed oxygen atom, and also nearly identical to that in Pt(111)-p(2x2)-O, apart from in the region of the chemisorbed CO, and (iii) the chemisorption energy of CO on a precovered Pt(111)-p(2x2)-O and the chemisorption energy of O on a precovered Pt(111)-p(2x2)CO are almost equal to that in Pt(111)-p(2x2)-CO and Pt(111)-p(2x2)-O, respectively. These results indicate that the interaction between CO and chemisorbed oxygen on a metal surface is mainly shore range in nature. The discussions of Pt-CO and Pt-O bonding and the interaction between CO and the chemisorbed oxygen atom on Pt(111) are augmented by local densities of states and real space distributions of quantum states.

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.

CO oxidation on Pt(111):An ab initio density functional theory study

CO oxidation on Pt(111) is studied with ab initio density functional theory. The low energy pathway and transition state for the reaction are identified. The key event is the breaking of an O-metal bond prior to the formation of a chemisorbed CO2 molecule. The pathway can be rationalized in terms of competition of the O and C atoms for bonding with the underlying surface, and the predominant energetic barrier is the strength of the O-metal bond.

GRADIENT CORRECTIONS IN DENSITY-FUNCTIONAL THEORY CALCULATIONS FOR SURFACES - CO ON PD(110)

Ab initio total energy calculations have been performed for CO chemisorption on Pd(110). Local density approximation (LDA) calculations yield chemisorption energies which are significantly higher than experimental values but inclusion of the generalised gradient approximation (GGA) gives better agreement. In general, sites with higher coordination of the adsorbate to surface atoms lead to a larger degree of overbinding with LDA, and give larger corrections with GGA. The reason is discussed using a first-order perturbation approximation. It is concluded that this may be a general failure of LDA for chemisorption energy calculations. This conclusion may be extended to many surface calculations, such as potential energy surfaces for diffusion.

Elucidating the mechanism and active site of the cyclohexanol dehydrogenation on copper-based catalysts: A density functional theory study

Abstract The dehydrogenation of cyclohexanol to cyclohexanone is very important in the manufacture of nylon. Copper-based catalysts are the most popular catalysts for this reaction, and on these catalysts the reaction mechanism and active site are in debate. In order to elucidate the mechanism and active site of the cyclohexanol dehydrogenation on copper-based catalysts, density functional theory with dispersion corrections were performed on up to six facets of copper in two different oxidation states: monovalent copper and metallic copper. By calculating the surface energies of these facets, Cu(111) and Cu2O(111) were found to be the most stable facets for metallic copper and for monovalent copper, respectively. On these two facets, all the possible elementary steps in the dehydrogenation pathway of cyclohexanol were calculated, including the adsorption, dehydrogenation, hydrogen coupling and desorption. Two different reaction pathways for dehydrogenation were considered on both surfaces. It was revealed that the dehydrogenation mechanisms are different on these two surfaces: on Cu(111) the hydrogen belonging to the hydroxyl is removed first, then the hydrogen belonging to the carbon is subtracted, while on Cu2O(111) the hydrogen belonging to the carbon is removed followed by the subtraction of the hydrogen in the hydroxyl group. Furthermore, by comparing the energy profiles of these two surfaces, Cu2O(111) was found to be more active for cyclohexanol dehydrogenation than Cu(111). In addition, we found that the coordinatively unsaturated copper sites on Cu2O(111) are the reaction sites for all the steps. Therefore, the coordinatively unsaturated copper site on Cu2O(111) is likely to be the active site for cyclohexanol dehydrogenation on the copper-based catalysts.

Some Attempts in the Rational Design of Heterogeneous Catalysts Using Density Functional Theory Calculations

Heterogeneous catalysis is of great importance both industrially and academically. Rational design of heterogeneous catalysts is highly desirable, and the computational screening and design method is one of the most promising approaches for rational design of heterogeneous catalysts. Herein, we review some attempts towards the rational catalyst design using density functional theory from our group. Some general relationships and theories on the activity and selectivity are covered, such as the Brønsted–Evans–Polanyi relation, volcano curves/surfaces, chemical potentials, optimal adsorption energy window and energy descriptor of selectivity. Furthermore, the relations of these relationships and theories to the rational design are discussed, and some examples of computational screening and design method are given.

Understanding catalytic reactions over zeolites: A density functional theory study of selective catalytic reduction of NOX by NH3 over Cu-SAPO-34

Metal exchanged CHA-type (SAPO-34 and SSZ-13) zeolites are promising catalysts for selective catalytic reduction (SCR) of NOx by NH3. However, the understanding of the process at the molecular level is still limited, which hinders the identification of its mechanism and the design of more efficient zeolite catalysts. In this work, modelling the reaction over Cu-SAPO-34, a periodic density functional theory (DFT) study of NH3-SCR was performed using hybrid functional with the consideration of van der Waals (vdW) interactions. A mechanism with a low N–N coupling barrier is proposed to account for the activation of NO. The redox cycle of Cu2+ and Cu+, which is crucial for the SCR process, is identified with detailed analyses. Besides, the decomposition of NH2NO is shown to readily occur on the Brønsted acid site by a hydrogen push-pull mechanism, confirming the collective efforts of Brønsted acid and Lewis acid (Cu2+) sites. The special electronic and structural properties of Cu-SAPO-34 are demonstrated to play an essential role the reaction, which may have a general implication on the understanding of zeolite catalysis.

Dielectrics in a time-dependent electric field: a real-time approach based on density-polarization functional theory

In the presence of a (time-dependent) macroscopic electric field the electron dynamics of dielectrics cannot be described by the time-dependent density only. We present a real-time formalism that has the density and the macroscopic polarization P as key quantities. We show that a simple local function of P already captures long-range correlation in linear and nonlinear optical response functions. Specifically, after detailing the numerical implementation, we examine the optical absorption, the second- and third-harmonic generation of bulk Si, GaAs, AlAs and CdTe at different level of approximation. We highlight links with ultranonlocal exchange-correlation functional approximations proposed within linear response time-dependent density functional theory framework.

Electronic Structure of an XUV Photogenerated Solid-Density Aluminum Plasma

By use of high intensity XUV radiation from the FLASH free-electron laser at DESY, we have created highly excited exotic states of matter in solid-density aluminum samples. The XUV intensity is sufficiently high to excite an inner-shell electron from a large fraction of the atoms in the focal region. We show that soft-x-ray emission spectroscopy measurements reveal the electronic temperature and density of this highly excited system immediately after the excitation pulse, with detailed calculations of the electronic structure, based on finite-temperature density functional theory, in good agreement with the experimental results.

On the required shape corrections to the local density and generalized gradient approximations to the Kohn-Sham potentials for molecular response calculations of (hyper)polarizabilities and excitation energies

It is well known that shape corrections have to be applied to the local-density (LDA) and generalized gradient (GGA) approximations to the Kohn-Sham exchange-correlation potential in order to obtain reliable response properties in time dependent density functional theory calculations. Here we demonstrate that it is an oversimplified view that these shape corrections concern primarily the asymptotic part of the potential, and that they affect only Rydberg type transitions. The performance is assessed of two shape-corrected Kohn-Sham potentials, the gradient-regulated asymptotic connection procedure applied to the Becke-Perdew potential (BP-GRAC) and the statistical averaging of (model) orbital potentials (SAOP), versus LDA and GGA potentials, in molecular response calculations of the static average polarizability alpha, the Cauchy coefficient S-4, and the static average hyperpolarizability beta. The nature of the distortions of the LDA/GGA potentials is highlighted and it is shown that they introduce many spurious excited states at too low energy which may mix with valence excited states, resulting in wrong excited state compositions. They also lead to wrong oscillator strengths and thus to a wrong spectral structure of properties like the polarizability. LDA, Becke-Lee-Yang-Parr (BLYP), and Becke-Perdew (BP) characteristically underestimate contributions to alpha and S-4 from bound Rydberg-type states and overestimate those from the continuum. Cancellation of the errors in these contributions occasionally produces fortuitously good results. The distortions of the LDA, BLYP, and BP spectra are related to the deficiencies of the LDA/GGA potentials in both the bulk and outer molecular regions. In contrast, both SAOP and BP-GRAC potentials produce high quality polarizabilities for 21 molecules and also reliable Cauchy moments and hyperpolarizabilities for the selected molecules. The analysis for the N-2 molecule shows, that both SAOP and BP-GRAC yield reliable energies omega(i) and oscillator strengths f(i) of individual excitations, so that they reproduce well the spectral structure of alpha and S-4.(C) 2002 American Institute of Physics.