Electronic structure of point defects in controlled self-doping of the TiO2 (110) surface: Combined photoemission spectroscopy and density functional theory study

Point defects in metal oxides such as TiO2 are key to their applications in numerous technologies. The investigation of thermally induced nonstoichiometry in TiO2 is complicated by the difficulties in preparing and determining a desired degree of nonstoichiometry. We study controlled self-doping of TiO2 by adsorption of 1/8 and 1/16 monolayer Ti at the (110) surface using a combination of experimental and computational approaches to unravel the details of the adsorption process and the oxidation state of Ti. Upon adsorption of Ti, x-ray and ultraviolet photoemission spectroscopy (XPS and UPS) show formation of reduced Ti. Comparison of pure density functional theory (DFT) with experiment shows that pure DFT provides an inconsistent description of the electronic structure. To surmount this difficulty, we apply DFT corrected for on-site Coulomb interaction (DFT+U) to describe reduced Ti ions. The optimal value of U is 3 eV, determined from comparison of the computed Ti 3d electronic density of states with the UPS data. DFT+U and UPS show the appearance of a Ti 3d adsorbate-induced state at 1.3 eV above the valence band and 1.0 eV below the conduction band. The computations show that the adsorbed Ti atom is oxidized to Ti2+ and a fivefold coordinated surface Ti atom is reduced to Ti3+, while the remaining electron is distributed among other surface Ti atoms. The UPS data are best fitted with reduced Ti2+ and Ti3+ ions. These results demonstrate that the complexity of doped metal oxides is best understood with a combination of experiment and appropriate computations.

The molecular structure of hexamethyldigermane determined by gas-phase electron diffraction with theoretical calculations for (CH3)3M-M(CH3)3 Where M = C, Si, and Ge

Gas-phase electron diffraction (GED) data together with results from ab initio molecular orbital calculations (HF and MP2/6-311+G(d,p)) have been used to determine the structure of hexamethyldigermane ((CH3)3Ge-Ge(CH3)3). The equilibrium symmetry is D3d, but the molecule has a very low-frequency, largeamplitude, torsional mode (φCGeGeC) that lowers the thermal average symmetry. The effect of this largeamplitude mode on the interatomic distances was described by a dynamic model which consisted of a set of pseudoconformers spaced at even intervals. The amount of each pseudoconformer was obtained from the ab initio calculations (HF/6-311+G(d,p)). The results for the principal distances (ra) and angles (∠h1) obtained from the combined GED/ab initio (with estimated 1σ uncertainties) are r(Ge-Ge) ) 2.417(2) Å, r(Ge-C) ) 1.956(1) Å, r(C-H) ) 1.097(5) Å, ∠GeGeC ) 110.5(2)°, and ∠GeCH ) 108.8(6)°. Theoretical calculations were performed for the related molecules ((CH3)3Si-Si(CH3)3 and (CH3)3C-C(CH3)3).

Ab initio prediction of the infrared-absorption spectrum of the C2Cl radical

The three lowest (1(2)A('), 2(2)A('), and 1(2)A(')) potential-energy surfaces of the C2Cl radical, correlating at linear geometries with (2)Sigma(+) and (2)Pi states, have been studied ab initio using a large basis set and multireference configuration-interaction techniques. The electronic ground state is confirmed to be bent with a very low barrier to linearity, due to the strong nonadiabatic electronic interactions taking place in this system. The rovibronic energy levels of the (CCCl)-C-12-C-12-Cl-35 isotopomer and the absolute absorption intensities at a temperature of 5 K have been calculated, to an upper limit of 2000 cm(-1), using diabatic potential-energy and dipole moment surfaces and a recently developed variational method. The resulting vibronic states arise from a strong mixture of all the three electronic components and their assignments are intrinsically ambiguous. (c) 2005 American Institute of Physics.

Ab initio prediction of the infrared absorption spectrum of the C2Br radical

The first three electronic states (1(2)A', 2(2)A', 1(2)A '') of the C2Br radical, correlating at linear geometries with (2)Sigma(+) and (2)Pi states, have been studied ab initio, using Multi Reference Configuration Interaction techniques. The electronic ground state is found to have a bent equilibrium geometry, R-CC = 1.2621 angstrom, R-CBr = 1.7967 angstrom, < CCBr 156.1 degrees, with a very low barrier to linearity. Similarly to the valence isoelectronic radicals C2F and C2Cl, this anomalous behaviour is attributed to a strong three-state non-adiabatic electronic interaction. The Sigma, Pi(1/2), Pi(3/2) vibronic energy levels and their absolute infrared absorption intensities at a temperature of 5K have been calculated for the (CCBr)-C-12-C-12-Br-79 isotopomer, to an upper limit of 2000 cm(-1), using ab initio diabatic potential energy and dipole moment surfaces and a recently developed variational method.

The molecular structure of propylene sulphide (methylthiirane) by gas-phase electron diffraction and theoretical calculations: a molecule used in MOCVD

Gas-phase electron-diffraction (GED) data together with results from ab initio molecular orbital calculations have been used to determine the structure of propylene sulphide. Values found for the main structural parameters for the molecule are consistent with those obtained from microwave studies and are compared here with those found for similar sulphur containing rings of general formula S(CH2)n (n = 2–5). A high ring strain enthalpy was calculated for propylene sulphide which is consistent with the small C–S–C angle (48.2(6)degrees) and the relatively long C–S bond lengths (ra = 1.831(2) Å). This is thought to account for the ease of ring opening in propylene sulphide observed in MOCVD reactions and the ready polymerisation of the molecule.

Unraveling the electronic structures of low-valent naphthalene and anthracene iron complexes: x-ray, spectroscopic, and density functional theory studies

Naphthalene and anthracene transition metalates are potent reagents, but their electronic structures have remained poorly explored. A study of four Cp*-substituted iron complexes (Cp* = pentamethylcyclopentadienyl) now gives rare insight into the bonding features of such species. The highly oxygen- and water-sensitive compounds [K(18-crown- 6){Cp*Fe(η4-C10H8)}] (K1), [K(18-crown-6){Cp*Fe(η4-C14H10)}] (K2), [Cp*Fe(η4-C10H8)] (1), and [Cp*Fe(η4-C14H10)] (2) were synthesized and characterized by NMR, UV−vis, and 57Fe Mössbauer spectroscopy. The paramagnetic complexes 1 and 2 were additionally characterized by electron paramagnetic resonance (EPR) spectroscopy and magnetic susceptibility measurements. The molecular structures of complexes K1, K2, and 2 were determined by single-crystal X-ray crystallography. Cyclic voltammetry of 1 and 2 and spectroelectrochemical experiments revealed the redox properties of these complexes, which are reversibly reduced to the monoanions [Cp*Fe(η4-C10H8)]− (1−) and [Cp*Fe(η4-C14H10)]− (2−) and reversibly oxidized to the cations [Cp*Fe(η6-C10H8)]+ (1+) and [Cp*Fe(η6-C14H10)]+ (2+). Reduced orbital charges and spin densities of the naphthalene complexes 1−/0/+ and the anthracene derivatives 2−/0/+ were obtained by density functional theory (DFT) methods. Analysis of these data suggests that the electronic structures of the anions 1− and 2− are best represented by low-spin FeII ions coordinated by anionic Cp* and dianionic naphthalene and anthracene ligands. The electronic structures of the neutral complexes 1 and 2 may be described by a superposition of two resonance configurations which, on the one hand, involve a low-spin FeI ion coordinated by the neutral naphthalene or anthracene ligand L, and, on the other hand, a low-spin FeII ion coordinated to a ligand radical L•−. Our study thus reveals the redox noninnocent character of the naphthalene and anthracene ligands, which effectively stabilize the iron atoms in a low formal, but significantly higher spectroscopic oxidation state.

Density functional theory study of rutile VO2 surfaces

We present the results of a density functional theory (DFT) investigation of the surfaces of rutile-like vanadium dioxide, VO2(R). We calculate the surface energies of low Miller index planes, and find that the most stable surface orientation is the (110). The equilibrium morphology of a VO2(R) particle has an acicular shape, laterally confined by (110) planes and topped by (011) planes. The redox properties of the (110) surface are investigated by calculating the relative surface free energies of the non-stoichiometric compositions as a function of oxygen chemical potential. It is found that the VO2(110) surface is oxidized with respect to the stoichiometric composition, not only at ambient conditions but also at the more reducing conditions under which bulk VO2 is stable in comparison with bulk V2O5. The adsorbed oxygen forms surface vanadyl species much more favorably than surface peroxo species.

Algorithm for generating defective graphene sheets

An algorithm is presented for the generation of molecular models of defective graphene fragments, containing a majority of 6-membered rings with a small number of 5- and 7-membered rings as defects. The structures are generated from an initial random array of points in 2D space, which are then subject to Delaunay triangulation. The dual of the triangulation forms a Voronoi tessellation of polygons with a range of ring sizes. An iterative cycle of refinement, involving deletion and addition of points followed by further triangulation, is performed until the user-defined criteria for the number of defects are met. The array of points and connectivities are then converted to a molecular structure and subject to geometry optimization using a standard molecular modeling package to generate final atomic coordinates. On the basis of molecular mechanics with minimization, this automated method can generate structures, which conform to user-supplied criteria and avoid the potential bias associated with the manual building of structures. One application of the algorithm is the generation of structures for the evaluation of the reactivity of different defect sites. Ab initio electronic structure calculations on a representative structure indicate preferential fluorination close to 5-ring defects.

What have we learnt about heavy carbenes through laser flash photolysis studies?

Time resolved gas-phase kinetic studies have contributed a great deal of fundamental information about the reactions and reactivity of heavy carbenes (silylenes, germylenes and stannylenes) during the past two decades. In this article we trace the development of our understanding through the mechanistic themes of intermediate complexes, third body assisted associations, catalysed reactions, non-observed reactions and substituent effects. Ab initio (quantum chemical) calculations have substantially assisted mechanistic interpretation and are discussed where appropriate. Trends in reactivity are identified and some signposts to future studies are indicated. This review, although detailed, is not comprehensive.

Substituent effects in single helix-double helix interconversions in copper(I) complexes

The copper(I) complex of L, the 1:2 condensate of benzil dihydrazone and 2-formylpyridine, exists as single, helical [CuL](+) and double helical [Cu2L2](2+) in dichloromethane solution but crystallizes only as the double helicate [Cu2L2](ClO4)(2). In contrast, earlier [New J Chem, 27 (2003) 193] it has been found that with L', the 1:2 condensate of benzil dihydrazone and 2-acetylpyridine, only the single helical monomeric species [CuL'](+) is isolable as solid. This contrasting behaviour of the copper(I) complexes of L and L' are scrutinised here by density functional calculations.

Metallo-organic domino reactions: C-H versus C-C bond breaking

HL and MeL are prepared by condensing benzil dihydrazone with 2-formylpyridine and 2-acetylpyridine, respectively, in 1:2 molar proportions. While in a reaction with [Ru-(C6H6)Cl-2](2), HL yields the cation [Ru(C6H6){5,6-diphenyl-3-(pyridin-2-yl)- 1,2,4-triazine}Cl](+), MeL gives the cation [Ru(C6H6)(MeL)Cl](+). Both the cations are isolated as their hexafluorophosphate salts and characterised by X-ray crystallography. In the case of HL, double domino electrocyclic/elimination reactions are found to occur. The electrocyclic reaction occurs in a C=N-N=C-C=N fragment of HL and the elimination reaction involves breaking of a C-H bond of HL. Density functional calculations on model complexes indicate that the identified electrocyclic reaction is thermochemically as well as kinetically feasible for both HL and MeL in the gas phase. For a double domino reaction, similar to that operative in HL, to occur for MeL, breaking of a C-C bond would be required in the elimination step. Our model calculations show the energy barrier for this elimination step to be much higher (329.1 kJ mol(-1)) for MeL than that for HL (96.3 kJ mol(-1)). Thus, the domino reaction takes place for HL and not for MeL. This accounts for the observed stability of [Ru(C6H6)-(MeL)Cl](+) under the reaction conditions employed.

Non-stoichiometric oxide and metal interfaces and reactions

We have employed a combination of experimental surface science techniques and density functional calculations to study the reduction of TiO2(110) surfaces through the doping with submonolayer transition metals. We concentrate on the role of Ti adatoms in self doping of rutile and contrast the behaviour to that of Cr. DFT+U calculations enable identification of probable adsorption structures and their spectroscopic characteristics. Adsorption of both metals leads to a broken symmetry and an asymmetric charge transfer localised around the defect site of a mixed localised/delocalised character. Charge transfer creates defect states with Ti 3d character in the band gap at similar to 1-eV binding energy. Cr adsorption, however, leads to a very large shift in the valence-band edge to higher binding energy and the creation of Cr 3d states at 2.8-eV binding energy. Low-temperature oxidation lifts the Ti-derived band-gap states and modifies the intensity of the Cr features, indicative of a change of oxidation state from Cr3+ to Cr4+. Higher temperature processing leads to a loss of Cr from the surface region, indicative of its substitution into the bulk.

Transition metal clusters on graphite

The topic of this work is 3d transition metals deposited on graphite. Spin-polarised density-functional calculations are used to obtain the magnetic moments of deposited adatoms and dimers. Interatomic potentials are also deduced. These are used in molecular dynamics simulations to study cluster formation and to investigate cluster morphology.

On the difficulties of present theoretical models to predict the oxidation state of atomic Au adsorbed on regular sites of CeO2(111)

The electronic structure and oxidation state of atomic Au adsorbed on a perfect CeO2(111) surface have been investigated in detail by means of periodic density functional theory-based calculations, using the LDA+U and GGA+U potentials for a broad range of U values, complemented with calculations employing the HSE06 hybrid functional. In addition, the effects of the lattice parameter a0 and of the starting point for the geometry optimization have also been analyzed. From the present results we suggest that the oxidation state of single Au atoms on CeO2(111) predicted by LDA+U, GGA+U, and HSE06 density functional calculations is not conclusive and that the final picture strongly depends on the method chosen and on the construction of the surface model. In some cases we have been able to locate two well-defined states which are close in energy but with very different electronic structure and local geometries, one with Au fully oxidized and one with neutral Au. The energy difference between the two states is typically within the limits of the accuracy of the present exchange-correlation potentials, and therefore, a clear lowest-energy state cannot be identified. These results suggest the possibility of a dynamic distribution of Au0 and Au+ atomic species at the regular sites of the CeO2(111) surface.