Structural manifestation of the delocalization error of density functional approximations: C(4N+2) rings and C(20) bowl, cage, and ring isomers.

The ground state structure of C(4N+2) rings is believed to exhibit a geometric transition from angle alternation (N < or = 2) to bond alternation (N > 2). All previous density functional theory (DFT) studies on these molecules have failed to reproduce this behavior by predicting either that the transition occurs at too large a ring size, or that the transition leads to a higher symmetry cumulene. Employing the recently proposed perspective of delocalization error within DFT we rationalize this failure of common density functional approximations (DFAs) and present calculations with the rCAM-B3LYP exchange-correlation functional that show an angle-to-bond-alternation transition between C(10) and C(14). The behavior exemplified here manifests itself more generally as the well known tendency of DFAs to bias toward delocalized electron distributions as favored by Huckel aromaticity, of which the C(4N+2) rings provide a quintessential example. Additional examples are the relative energies of the C(20) bowl, cage, and ring isomers; we show that the results from functionals with minimal delocalization error are in good agreement with CCSD(T) results, in contrast to other commonly used DFAs. An unbiased DFT treatment of electron delocalization is a key for reliable prediction of relative stability and hence the structures of complex molecules where many structure stabilization mechanisms exist.

Reformulating time-dependent density functional theory with non-orthogonal localized molecular orbitals.

Time-dependent density functional theory (TDDFT) has broad application in the study of electronic response, excitation and transport. To extend such application to large and complex systems, we develop a reformulation of TDDFT equations in terms of non-orthogonal localized molecular orbitals (NOLMOs). NOLMO is the most localized representation of electronic degrees of freedom and has been used in ground state calculations. In atomic orbital (AO) representation, the sparsity of NOLMO is transferred to the coefficient matrix of molecular orbitals (MOs). Its novel use in TDDFT here leads to a very simple form of time propagation equations which can be solved with linear-scaling effort. We have tested the method for several long-chain saturated and conjugated molecular systems within the self-consistent charge density-functional tight-binding method (SCC-DFTB) and demonstrated its accuracy. This opens up pathways for TDDFT applications to large bio- and nano-systems.

Experimental and density functional theory study of the vibrational properties of 2-mercaptobenzimidazole in interaction with gold

The vibrational properties of the 2-mercaptobenzimidazole (MBI) molecule in interaction with gold were examined by a combined approach of FTIR measurements and density functional theory (DFT). A complete assignment of the 42 normal modes of MBI has been performed on the basis of DFT calculations at the B3PW91 level in complement to the Raman and FTIR spectra. Calculations demonstrated that, on the deprotonated MBI molecule, the negative charge is localized on the sulfur atom, favoring the formation of a gold-sulfur bond upon reaction of MBI with gold. This was confirmed by the very good agreement between the calculated spectrum and the experimental spectra of different gold-MBI compounds, indicating that the vibrational properties of adsorbed MBI are chiefly determined by the coordination through the sulfur atom. © 2006 American Chemical Society.

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

Reactions of polycarbonate with cyclohexene oxide and phosphites: A density functional study

Epoxides and phosphites are often used as additives to stabilize the properties of polymers, including bisphenol A polycarbonate (BPA-PC). We describe density functional (DF) calculations of the reactions of cyclohexene oxide (CHO, cyclohexane epoxide) and phosphites with chain segments of BPA-PC, with the aim of identifying possible reaction paths and energy barriers. The reactions of CHO with the OH-terminated PC chains and with the carbonate group are exothermic, although there is an energy barrier in each case of more than 10 kcal/mol. A comparison of results for different CHO isomers demonstrates the importance of steric effects. The reactions between the same groups of the PC chain and the phosphites 2-[2,4-bis(tert-butyl)phenoxy]-5,5-dimethyl-1,3,2-dioxaphosphorinane] (BPDD) and trimethyl phosphite (TMP), and their phosphonate isomers are characterized by large energy barriers.

Understanding the dehydrogenation mechanism of tetrahydrocarbazole over palladium using a combined experimental and density functional theory approach

The mechanism of the dehydrogenation of tetrahydrocarbazole to carbazole over palladium has been examined for the first time. By use of a combination of deuterium exchange experiments and density functional theory calculations, a detailed reaction profile for the aromatization of tetrahydrocarbazole has been identified and validated by experiment. As with many dehydrogenation reactions, the initial hydrogen abstraction is found to have the highest reaction barrier. Tetrahydrocarbazole has four hydrogens which can, in principle, be cleaved initially; however, the theory and experiment show that the reaction is dominated by the cleavage of the carbon hydrogens at the carbon atoms in positions 1 and 4. The two pathways originating from these two C-H bond cleavage processes are found to have similar reaction energy profiles and both contribute to the overall reaction.

Au on (111) and (110) surfaces of CeO2: A density-functional theory study

The atomic structures of gold supported on (111) and (110) surfaces of CeO2 have been studied using density-functional theory calculations. A single Au atom is placed on three adsorption sites on the surfaces; the stoichiometric surfaces, an oxygen vacancy and a Ce-vacancy. It is found that (i) the Au adsorption energies are in the following order: E-ad(Ce-vacancy) > E-ad(O-vacancy) > E-ad(stoichiometric surface); and (ii) the Au atom adsorption on the Ce-vacancy activates O atoms nearby. One 0 atom is less stable than that in O-2 in the gas phase and another O atom is much easier to remove compared to that of the stoichiometric surfaces. These results suggest that the Au adsorption on Ce-vacancies not only creates an O-vacancy but also activates an O atom nearby. This provides a piece of direct evidence that Au adsorption on a Ce-vacancy may be responsible for some unique catalytic properties of Au/CeO2. (C) 2008 Elsevier B.V. All rights reserved.

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 sulfur poisoning

Density functional theory calculations have been used to investigate the chemisorption of H, S, SH, and H2S as well as the hydrogenation reactions S+H and SH+H on a Rh surface with steps, Rh(211), aiming to explain sulfur poisoning effect. In the S hydrogenation from S to H2S, the transition state of the first step S+H-->SH is reached when the S moves to the step-bridge and H is on the off-top site. In the second step, SH+H-->H2S, the transition state is reached when SH moves to the top site and H is close to another top site nearby. Our results show that it is difficult to hydrogenate S and they poison defects such as steps. In order to address why S is poisoning, hydrogenation of C, N, and O on Rh(211) has also been calculated and has been found that the reverse and forward reactions possess similar barriers in contrast to the S hydrogenation. The physical origin of these differences has been analyzed and discussed. (C) 2005 American Institute of Physics.

General insight into CO oxidation: A density functional theory study of the reaction mechanism on platinum oxides

CO oxidation on PtO2(110) has been studied using density functional theory calculations. Four possible reaction mechanisms were investigated and the most feasible one is the following: (i) the O at the bridge site of PtO2(110) reacts with CO on the coordinatively unsaturated site (CUS) with a negligible barrier; (ii) O-2 adsorbs on the bridge site and then interacts with CO on the CUS to form an OO-CO complex; (iii) the bond of O-OCO breaks to produce CO2 with a small barrier (0.01 eV). The CO oxidation mechanisms on metals and metal oxides are rationalized by a simple model: The O-surface bonding determines the reactivity on surfaces; it also determines whether the atomic or molecular mechanism is preferred. The reactivity on metal oxides is further found to be related to the 3rd ionization energy of the metal atom.

CO oxidation and NO reduction on metal surfaces: density functional theory investigations

This article reviews the accumulated theoretical results, in particular density functional theory calculations, on two catalytic processes, CO oxidation and NO reduction on metal surfaces. Owing to their importance in automotive emission control, these two reactions have generated a lot of interest in the last 20 years. Here the pathways and energetics of the involved elementary reactions under different catalytic conditions are described in detail and the understanding of the reactions is generalized. It is concluded that density functional theory calculations can be applied to catalysis to elucidate mechanisms of complex surface reactions and to understand the electronic structure of chemical processes in general. The achieved molecular knowledge of chemical reactions is certainly beneficial to new catalyst design.

Catalytic role of gold in gold-based catalysts: A density functional theory study on the CO oxidation on gold

Gold-based catalysts have been of intense interests in recent years, being regarded as a new generation of catalysts due to their unusually high catalytic performance. For example, CO oxidation on Au/TiO2 has been found to occur at a temperature as low as 200 K. Despite extensive studies in the field, the microscopic mechanism of CO oxidation on Au-based catalysts remains controversial. Aiming to provide insight into the catalytic roles of Au, we have performed extensive density functional theory calculations for the elementary steps in CO oxidation on Au surfaces. O atom adsorption, CO adsorption, O-2 dissociation, and CO oxidation on a series of Au surfaces, including flat surfaces, defects and small clusters, have been investigated in detail. Many transition states involved are located, and the lowest energy pathways are determined. We find the following: (i) the most stable site for O atom on Au is the bridge site of step edge, not a kink site; (ii) O-2 dissociation on Au (O-2-->20(ad)) is hindered by high barriers with the lowest barrier being 0.93 eV on a step edge; (iii) CO can react with atomic O with a substantially lower barrier, 0.25 eV, on Au steps where CO can adsorb; (iv) CO can react with molecular O-2 on Au steps with a low barrier of 0.46 eV, which features an unsymmetrical four-center intermediate state (O-O-CO); and (v) O-2 can adsorb on the interface of Au/TiO2 with a reasonable chemisorption energy. On the basis of our calculations, we suggest that (i) O-2 dissociation on Au surfaces including particles cannot occur at low temperatures; (ii) CO oxidation on Au/inactive-materials occurs on Au steps via a two-step mechanism: CO+O-2-->CO2+O, and CO+O-->CO2; and (iii) CO oxidation on Au/active-materials also follows the two-step mechanism with reactions occurring at the interface.

A density functional theory study of stepwise addition reactions in ammonia synthesis on Ru(0001)

To shed light on stepwise addition reactions in ammonia synthesis, density functional theory calculations are carried out to investigate NHx (x = 1-3) formation on Ru(0001). The reactions on a flat surface are first examined. Transition states and reaction barriers are determined. It is found that the reaction barriers for these 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 N-2 dissociation. One of the stepwise addition reactions, NH + H --> NH2, on a stepped surface is also considered. Interestingly, the reaction barrier is found to be significantly lower than that on the flat surface, but is considerably higher than that of N-2 dissociation on the same stepped surface. In addition, the coverage effect on the reaction energetics is also addressed. (C) 2001 Published by Elsevier Science B.V.