Identifying the distinct features of geometric structures for hole trapping to generate radicals on rutile TiO2(110) in photooxidation using density functional theory calculations with hybrid functional

Using density functional theory calculations with HSE 06 functional, we obtained the structures of spin-polarized radicals on rutile TiO2(110), which is crucial to understand the photooxidation at the atomic level, and further calculate the thermodynamic stabilities of these radicals. By analyzing the results, we identify the structural features for hole trapping in the system, and reveal the mutual effects among the geometric structures, the energy levels of trapped hole states and their hole trapping capacities. Furthermore, the results from HSE 06 functional are compared to those from DFT + U and the stability trend of radicals against the number of slabs is tested. The effect of trapped holes on two important steps of the oxygen evolution reaction, i.e. water dissociation and the oxygen removal, is investigated and discussed.

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.

Density Functional Theory Study of Iron and Cobalt Carbides for Fischer-Tropsch Synthesis

Carbides are important phases in heterogeneous catalysis. However, the understanding of carbide phases is inadequate: Fe and Co are the two commercial catalysts for Fischer-Tropsch (FT) synthesis, and experimental work showed that Fe carbide is the active phase in FT synthesis, whereas the appearance of Co carbide is considered as a possible deactivation cause, TO understand very different catalytic roles of carbides, all the key elementary steps in FT synthesis, that is, CO dissociation, C(1) hydrogenation, and C(1)+C(1) coupling, are extensively investigated on both carbide surfaces using first principles calculations. In particular, the most important issues in FT synthesis, the activity and methane selectivity, on the carbide surfaces are quantitatively determined and analyzed. They are also discussed together with metallic Fe and Co surfaces. It is found that (i) Fe carbide is more active than metallic Fe and has similar methane selectivity to Fe, being consistent with the experiments; and (ii) Co carbide is less active than Co and has higher methane selectivity, providing evidence on the molecular level to support the suggestion that the formation of Co carbide is a cause of relatively high methane selectivity and deactivation on Co catalysts.

An Energy Descriptor To Quantify Methane Selectivity in Fischer-Tropsch Synthesis:A Density Functional Theory Study

Selectivity is a fundamental issue in heterogeneous catalysis. In this study, the CH(4) selectivity in Fischer-Tropsch synthesis is chosen to be investigated: CH4 selectivity on Rh, Co, Ru, Fe, and Re surfaces is computed by first-principles methods. In conjunction with kinetic analyses, we are able to derive the effective barrier difference between methane formation and chain growth (Delta E(eff)) to quantify the CH(4) selectivity. By using this energy descriptor, the ranking of methane selectivity predicted from density functional theory (DFT) calculations is consistent with experimental work. Moreover, a linear correlation between Delta E(eff) and the chemisorption energy of C + 4H (Delta H) is found. This fundamental finding possesses the following significance: (i) it shows that the selectivity, which appears to have kinetic characteristics, is largely determined by thermodynamic properties; and (ii) it suggests that an increase of the binding strength of C + 4H will suppress methane selectivity.

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.

Helium in polygonal nanopores at zero temperature: Density functional theory calculations

We investigate adsorption of helium in nanoscopic polygonal pores at zero temperature using a finite-range density functional theory. The adsorption potential is computed by means of a technique denoted as the elementary source method. We analyze a rhombic pore with Cs walls, where we show the existence of multiple interfacial configurations at some linear densities, which correspond to metastable states. Shape transitions and hysterectic loops appear in patterns which are richer and more complex than in a cylindrical tube with the same transverse area.

Coherent tunnelling in Cu2+- and Ag2+-doped MgO and CaO: Cu2+ explored through ab initio calculations

The observation of coherent tunnelling in Cu2+ - and Ag2+ -doped MgO and CaO:Cu2+ was a crucial discovery in the realm of the Jahn-Teller (JT) effect. The main reasons favoring this dynamic behavior are now clarified through ab initio calculations on Cu2+ - and Ag2+ -doped cubic oxides. Small JT distortions and an unexpected low anharmonicity of the eg JT mode are behind energy barriers smaller than 25 cm-1 derived through CASPT2 calculations for Cu2+ - and Ag2+ -doped MgO and CaO:Cu2+ . The low anharmonicity is shown to come from a strong vibrational coupling of MO610- units (M=Cu,Ag) to the host lattice. The average distance between the d9 impurity and ligands is found to vary significantly on passing from MgO to SrO following to a good extent the lattice parameter.

Relativistic ab initio calculations for ion-atom collisions

Within the independent particle model we solve the time-dependent single-particle equation using ab initio SCF-DIRAC-FOCK-SLATER wavefunctions as a basis. To reinstate the many-particle aspect of the collision system we use the inclusive probability formalism to answer experimental questions. As an example we show an application to the case of S{^15+} on Ar where experimental data on the K-K charge transfer are available for a wide range of impact energies from 4.7 to 90 MeV. Our molecular adiabatic calculations and the evaluation using the inclusive probability formalism show good results in the low energy range from 4.7 to 16 MeV impact energy.

Exploring chromium (VI) dioxodihalides chemistry: is density functional theory the most suitable tool?

A comparative systematic study of the CrO2F2 compound has been performed using different conventional ab initio methodologies and density functional procedures. Two points have been analyzed: first, the accuracy of results yielded by each method under study, and second, the computational cost required to reach such results. Weighing up both aspects, density functional theory has been found to be more appropriate than the Hartree-Fock (HF) and the analyzed post-HF methods. Hence, the structural characterization and spectroscopic elucidation of the full CrO2X2 series (X=F,Cl,Br,I) has been done at this level of theory. Emphasis has been given to the unknown CrO2I2 species, and specially to the UV/visible spectra of all four compounds. Furthermore, a topological analysis in terms of charge density distributions has revealed why the valence shell electron pair repulsion model fails in predicting the molecular shape of such CrO2X2 complexes

Low-lying excited states and primary photoproducts of [Os-3(CO)(10)-(s-cis-L)] (L = cyclohexa-1,3-diene, buta-1,3-diene)] clusters studied by picosecond time-resolved UV/Vis and IR spectroscopy and by density functional theory

Combined picosecond transient absorption and time-resolved infrared studies were performed, aimed at characterising low-lying excited states of the cluster [Os-3(CO)(10)(s-cis-L)] (L= cyclohexa-1,3-diene, 1) and monitoring the formation of its photoproducts. Theoretical (DFT and TD-DFT) calculations on the closely related cluster with L=buta-1,3-diene (2') have revealed that the low-lying electronic transitions of these [Os-3(CO)(10)(s-cis-1,3-diene)] clusters have a predominant sigma(core)pi*(CO) character. From the lowest sigmapi* excited state, cluster 1 undergoes fast Os-Os(1,3-diene) bond cleavage (tau=3.3 ps) resulting in the formation of a coordinatively unsaturated primary photoproduct (1a) with a single CO bridge. A new insight into the structure of the transient has been obtained by DFT calculations. The cleaved Os-Os(1,3-diene) bond is bridged by the donor 1,3-diene ligand, compensating for the electron deficiency at the neighbouring Os centre. Because of the unequal distribution of the electron density in transient la, a second CO bridge is formed in 20 ps in the photoproduct [Os-3(CO)(8)(mu-CO)(2)- (cyclohexa-1,3-diene)] (1b). The latter compound, absorbing strongly around 630 nm, mainly regenerates the parent cluster with a lifetime of about 100 ns in hexane. Its structure, as suggested by the DFT calculations, again contains the 1,3-diene ligand coordinated in a bridging fashion. Photoproduct 1b can therefore be assigned as a high-energy coordination isomer of the parent cluster with all Os-Os bonds bridged.

Relationship between phonons and thermal expansion in Zn(CN)2 and Ni(CN)2 from inelastic neutron scattering and ab initio calculations

Zn(CN)2 and Ni(CN)2 are known for exhibiting anomalous thermal expansion over a wide temperature range. The volume thermal expansion coefficient for the cubic, three dimensionally connected material, Zn(CN)2, is negative (alpha(V) = −51  10(-6) K-1) while for Ni(CN)2, a tetragonal material, the thermal expansion coefficient is negative in the two dimensionally connected sheets (alpha(a) = −7  10(-6) K-1), but the overall thermal expansion coefficient is positive (alpha(V) = 48  10(-6) K-1). We have measured the temperature dependence of phonon spectra in these compounds and analyzed them using ab initio calculations. The spectra of the two compounds show large differences that cannot be explained by simple mass renormalization of the modes involving Zn (65.38 amu) and Ni (58.69 amu) atoms. This reflects the fact that the structure and bonding are quite different in the two compounds. The calculated pressure dependence of the phonon modes and of the thermal expansion coefficient, alpha(V), are used to understand the anomalous behavior in these compounds. Our ab initio calculations indicate that phonon modes of energy approx. 2 meV are major contributors to negative thermal expansion (NTE) in both the compounds. The low-energy modes of approx.8 and 13 meV in Zn(CN)2 also contribute significantly to the NTE in Zn(CN)2 and Ni(CN)2, respectively. The measured temperature dependence of the phonon spectra has been used to estimate the total anharmonicity of both compounds. For Zn(CN)2, the temperature-dependent measurements (total anharmonicity), along with our previously reported pressure dependence of the phonon spectra (quasiharmonic), is used to separate the explicit temperature effect at constant volume (intrinsic anharmonicity).

Electric field gradients in (111)In-doped (Hf/Zr)(3)Al(2) and (Hf/Zr)(4)Al(3) mixed compounds: ab initio calculations, perturbed angular correlation measurements and site preference

The quadrupolar hyperfine interactions of in-diffused (111)In -> (111)Cd probes in polycrystalline isostructural Zr(4)Al(3) and Hf(4)Al(3) samples containing small admixtures of the phases (Zr/Hf)(3)Al(2) were investigated. A strong preference of (111)In solutes for the contaminant (Zr/Hf)(3)Al(2) minority phases was observed. Detailed calculations of the electric field gradient (EFG) at the Cd nucleus using the full-potential augmented plane wave + local orbital formalism allowed us to assign the observed EFG fractions to the various lattice sites in the (Zr/Hf)(3)Al(2) compounds and to understand the preferential site occupation of the minority phases by the (111)In atoms. The effects of the size of the supercell and relaxation around the oversized In and Cd probe atoms were investigated in detail.

Electronic, vibrational and related properties of group IV metal oxides by ab initio calculations

We present our theoretical results for the structural, electronic, vibrational and optical properties of MO(2) (M = Sn, Zr, Hf and Ti) obtained by first-principles calculations. Relativistic effects are demonstrated to be important for a realistic description of the detailed structure of the electronic frequency-dependent dielectric function, as well as of the carrier effective masses. Based on our results, we found that the main contribution of the high values calculated for the oxides dielectric constants arises from the vibrational properties of these oxides, and the vibrational static dielectric constant values diminish with increasing pressure. (c) 2008 Elsevier B.V. All rights reserved.

Detailed analysis of the charge transfer complex N,N-dimethylaniline-SO(2) by Raman spectroscopy and density functional theory calculations

Although the amine sulfur dioxide chemistry was well characterized in the past both experimentally and theoretically, no systematic Raman spectroscopic study describes the interaction between N,N-dimethylaniline (DMA) and sulfur dioxide (SO(2)). The formation of a deep red oil by the reaction of SO(2) with DMA is an evidence of the charge transfer (CT) nature of the DMA-SO(2) interaction. The DMA -SO(2) normal Raman spectrum shows the appearance of two intense bands at 1110 and 1151 cm(-1), which are enhanced when resonance is approached. These bands are assigned to nu(s)(SO(2)) and nu(phi-N) vibrational modes, respectively, confirming the interaction between SO(2) and the amine via the nitrogen atom. The dimethyl group steric effect favors the interaction of SO(2) with the ring pi electrons, which gives rise to a pi-pi* low-energy CT electronic transition, as confirmed by time-dependent density functional theory (TDDFT) calculations. In addition, the calculated Raman DMA-SO(2) spectrum at the B3LYP/6-311++g(3df,3pd) level shows good agreement with the experimental results (vibrational wavenumbers and relative intensities), allowing a complete assignment of the vibrational modes. A better understanding of the intermolecular interactions in this model system can be extremely useful in designing new materials to absorb, detect, or even quantify SO(2). Copyright (C) 2009 John Wiley & Sons, Ltd.

Gaussian basis sets by generator coordinate Hartree-Fock method to ab initio calculations of electron affinities of enolates

The Generator Coordinate Hartree-Fock (GCHF) method is employed to generate uncontracted 15s and 18s11p gaussian basis sets for the H, C and O atoms, respectively. These basis sets are then contracted to 3s and 4s H atom and 6s5p, for C and O atoms by a standard procedure. For quality evaluation of contracted basis sets in molecular calculations, we have accomplished calculations of total and orbital energies in the Hartree-Fock-Roothaaii (HFR) approach for CH, C(2) and CO molecules. The results obtained with the uncontracted basis sets are compared with values obtained with the standard D95, 6-311G basis sets and with values reported in the literature. The 4s and 6s5p basis sets are enriched with polarization and diffuse functions for atoms of the parent neutral systems and of the enolates anions (cycloheptanone enolate, 2,5-dimethyleyelopentanone enolate, 4-heptanone enolate, and di-isopropyl ketone enolate) from the literature, in order to assess their performance in ab initio molecular calculations, and applied for calculations of electron affinities of the enolates. The calculations were performed at the DFT (BLYP and B3LYP) and HF levels and compared with the corresponding experimental values and with those obtained by using other 6-3 1 + +G((*)) and 6-311 + +G((*)) basis sets from literature. For the enolates studied, the differences between the electron affinities obtained with GCHF basis sets, at the B3LYP level, and the experimental values are -0.001, -0,014, -0.001, and -0.001 eV. (C) 2002 Elsevier B.V. B.V. All rights reserved.