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.

First-principles study of oxygenates on Co surfaces in Fischer-Tropsch synthesis

Extensive density function theory calculations are performed to study the mechanism of the formation of aldehyde and alcohol on Co surfaces in Fischer-Tropsch synthesis, a challenging issue in heterogeneous catalysis. Three possible pathways for the production of formaldehyde and methanol on flat and stepped Co(0001) surfaces are investigated: (i) CO + 4H -> CHO + 3H -> CH2O + 2H -> CH3O + H -> CH3OH; (ii) CO + 4H -> COH + 3H -> CHOH + 2H -> CH2OH + H -> CH3OH; and (iii) the coupling reactions of CH2 + O -> CH2O and CH3 + OH -> CH3OH. It is found that these pathways are generally favored at step sites, and the preferred mechanism is pathway (i) via CHO. Furthermore, the three traditional chain growth mechanisms in Fischer-Tropsch synthesis are semi quantitatively compared and discussed. Our results suggest that the two mechanisms involving oxygenate intermediates (the CO-insertion and hydroxycarbene mechanisms) are less important than the carbene mechanism in the production of long chain hydrocarbons. However, the CO-insertion mechanism may be responsible for the production of long-chain oxygenates.

General trends in CO dissociation on transition metal surfaces

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. (C) 2001 American Institute of Physics.

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.

Resolution of an ancient surface science anomaly: Work function change induced by N adsorption on W{100}

For many decades it has been assumed that an adsorbate centered above a metal surface and with a net negative charge should increase the work function of the surface. However, despite their electronegativity, N adatoms on W{100} cause a significant work function decrease. Here we present a resolution of this anomaly. Using density functional theory, we demonstrate that while the N atom carries a negative charge, of overriding importance is a reduction in the surface overspill electron density into the vacuum, when that charge is engaged in bonding to the adatom. This novel interpretation is fundamentally important in the general understanding of work function changes induced by atomic adsorbates.

Large-scale density functional calculations of the surface properties of the Wigner crystal

The surface properties of the jellium model have been investigated by large supercell computations in the density functional theory-local spin-density (DFT-LSD) approach for planar slabs with up to 1000 electrons. A wide interval of densities has been explored, extending into the stability range of the Wigner crystal. Most computations have been carried out on nominally paramagnetic samples with an equal number of spin-up and spin-down electrons. The results show that within DFT-LSD spontaneous spin polarization and charge localization start nearly simultaneously at the surface for r(s) similar to 20, then, with decreasing density, they progress toward the center of the slab. Electrons are fully localized and spin polarized at r(s) = 30. At this density the charge distribution is the superposition of disjoint charge blobs, each corresponding to one electron. The distribution of blobs displays both regularities and disorder, the first being represented by well-defined planes and simple in-plane geometries, and the latter by a variety of surface defects. The surface energy, surface dipole, electric polarisability, and magnetization pattern have been determined as a function of density. All these quantities display characteristic anomalies at the density of the localization transition. The analysis of the low-frequency electric conductivity shows that in the fluid paramagnetic regime the in-plane current preferentially flows in the central region of the slab and the two spin channels are equally conducting. In the charge localized, spin-polarized regime, conductivity is primarily a surface effect, and an apparent asymmetry is observed in the two spin currents.

Theory of the near K-edge structure in electron energy loss spectroscopy

Arguments are given that lead to a formalism for calculating near K-edge structure in electron energy loss spectroscopy (EELS). This is essentially a one electron picture, while many body effects may be introduced at different levels, such as the local density approximation to density functional theory or the GW approximation to the electron self-energy. Calculations are made within the all electron LMTO scheme in crystals with complex atomic and electronic structures, and these are compared with experiment. (c) 2004 Elsevier B.V. All rights reserved.

Escape times for rigid Brownian rotators in a bistable potential from the time evolution of the Green function and the characteristic time of the probability evolution

The greatest relaxation time for an assembly of three- dimensional rigid rotators in an axially symmetric bistable potential is obtained exactly in terms of continued fractions as a sum of the zero frequency decay functions (averages of the Legendre polynomials) of the system. This is accomplished by studying the entire time evolution of the Green function (transition probability) by expanding the time dependent distribution as a Fourier series and proceeding to the zero frequency limit of the Laplace transform of that distribution. The procedure is entirely analogous to the calculation of the characteristic time of the probability evolution (the integral of the configuration space probability density function with respect to the position co-ordinate) for a particle undergoing translational diffusion in a potential; a concept originally used by Malakhov and Pankratov (Physica A 229 (1996) 109). This procedure allowed them to obtain exact solutions of the Kramers one-dimensional translational escape rate problem for piecewise parabolic potentials. The solution was accomplished by posing the problem in terms of the appropriate Sturm-Liouville equation which could be solved in terms of the parabolic cylinder functions. The method (as applied to rotational problems and posed in terms of recurrence relations for the decay functions, i.e., the Brinkman approach c.f. Blomberg, Physica A 86 (1977) 49, as opposed to the Sturm-Liouville one) demonstrates clearly that the greatest relaxation time unlike the integral relaxation time which is governed by a single decay function (albeit coupled to all the others in non-linear fashion via the underlying recurrence relation) is governed by a sum of decay functions. The method is easily generalized to multidimensional state spaces by matrix continued fraction methods allowing one to treat non-axially symmetric potentials, where the distribution function is governed by two state variables. (C) 2001 Elsevier Science B.V. All rights reserved.