A DFT study of the structures, stabilities and redox behaviour of the major surfaces of magnetite Fe3O4

The renewed interest in magnetite (Fe3O4) as a major phase in different types of catalysts has led us to study the oxidation–reduction behaviour of its most prominent surfaces. We have employed computer modelling techniques based on the density functional theory to calculate the geometries and surface free energies of a number of surfaces at different compositions, including the stoichiometric plane, and those with a deficiency or excess of oxygen atoms. The most stable surfaces are the (001) and (111), leading to a cubic Fe3O4 crystal morphology with truncated corners under equilibrium conditions. The scanning tunnelling microscopy images of the different terminations of the (001) and (111) stoichiometric surfaces were calculated and compared with previous reports. Under reducing conditions, the creation of oxygen vacancies in the surface leads to the formation of reduced Fe species in the surface in the vicinity of the vacant oxygen. The (001) surface is slightly more prone to reduction than the (111), due to the higher stabilisation upon relaxation of the atoms around the oxygen vacancy, but molecular oxygen adsorbs preferentially at the (111) surface. In both oxidized surfaces, the oxygen atoms are located on bridge positions between two surface iron atoms, from which they attract electron density. The oxidised state is thermodynamically favourable with respect to the stoichiometric surfaces under ambient conditions, although not under the conditions when bulk Fe3O4 is thermodynamically stable with respect to Fe2O3. This finding is important in the interpretation of the catalytic properties of Fe3O4 due to the presence of oxidised species under experimental conditions.

Surface and interstitial Ti diffusion at the rutile TiO2(110) surface

Diffusion of Ti through the TiO2 (110) rutile surface plays a key role in the growth and reactivity of TiO2. To understand the fundamental aspects of this important process, we present an analysis of the diffusion of Ti adspecies at the stoichiometric TiO2(110) surface using complementary computational methodologies of density functional theory corrected for on-site Coulomb interactions (DFT+U) and a charge equilibration (QEq) atomistic potential to identify minimum energy pathways. We find that diffusion of Ti from the surface to subsurface (and vice versa) follows an intersticialcy exchange mechanism, involving exchange of surface Ti with the 6-fold coordinated Ti below the bridging oxygen rows. Diffusion in the subsurface between layers also follows an interstitialcy mechanism. The diffusion of Ti is discussed in light of continued attempts to understand the re-oxidation of non-stoichiometric TiO2(110) surfaces.

Searching efficiency and the functional response of a pause-travel forager

1. The feeding rates of many predators and parasitoids exhibit type II functional responses, with a decelerating rate of increase to reach an asymptotic value as the density of their prey or hosts increases. Holling's disc equation describes such relationships and predicts that the asymptotic feeding rate at high prey densities is set by handling time, while the rate at which feeding rate increases with increased prey density is determined by searching efficiency. Searching efficiency and handling time are also parameters in other models which describe the functional response. Models which incorporate functional responses in order to make predictions of the effects of food shortage thus rely upon a clear understanding and accurate quantification of searching efficiency and handling time. 2. Blackbird Turdus merula exhibit a type II functional response and use pause-travel foraging, a foraging technique in which animals search for prey while stationary and then move to capture prey. Pause-travel foraging allows accurate direct measurement of feeding rate and both searching efficiency and handling time. We use Blackbirds as a model species to: (i) compare observed measures of both searching efficiency and handling time with those estimated by statistically fitting the disc equation to the observed functional response; and (ii) investigate alternative measures of searching efficiency derived by the established method where search area is assumed to be circular and a new method that we propose where it is not. 3. We find that the disc equation can adequately explain the functional response of blackbirds feeding on artificial prey. However, this depends critically upon how searching efficiency is measured. Two variations on the previous method of measuring search area (a component of searching efficiency) overestimated searching efficiency, and hence predicted feeding rates higher than those observed. Two variations of our alternative approach produced lower estimates of searching efficiency, closer to that estimated by fitting the disc equation, and hence more accurately predicted feeding rate. Our study shows the limitations of the previous method of measuring searching efficiency, and describes a new method for measuring searching efficiency more accurately.

The vibrational spectrum and ultimate modulus of polyethylene

We have performed the first completely ab initio lattice dynamics calculation of the full orthorhombic cell of polyethylene using periodic density functional theory in the local density approximation (LDA) and the generalized gradient approximation (GGA). Contrary to current perceptions, we show that LDA accurately describes the structure whereas GGA fails. We emphasize that there is no parametrization of the results. We then rigorously tested our calculation by computing the phonon dispersion curves across the entire Brillouin zone and comparing them to the vibrational spectra, in particular the inelastic neutron scattering (INS) spectra, of polyethylene (both polycrystalline and aligned) and perdeuteriopolyethylene. The F-point frequencies (where the infrared and Raman active modes occur) are in good agreement with the latest low temperature data. The near-perfect reproduction of the INS spectra, gives confidence in the results and allows Lis to deduce a number of physical properties including the elastic moduli, parallel and perpendicular to the chain. We find that the Young's modulus for an infinitely long, perfectly crystalline polyethylene is 360.2 GPa at 0 K. The highest experimental value is 324 GPa, indicating that current high modulus fibers are similar to 90% of their maximum possible strength.

Gas phase kinetic and quantum chemical studies of the reactions of silylene with the methylsilanes. Absolute rate constants, temperature dependences, RRKM modelling and potential energy surfaces

Time resolved studies of silylene, SiH2, generated by the 193 nm laser. ash photolysis of phenylsilane, have been carried out to obtain rate coefficients for its bimolecular reactions with methyl-, dimethyl- and trimethyl-silanes in the gas phase. The reactions were studied over the pressure range 3 - 100 Torr with SF6 as bath gas and at five temperatures in the range 300 - 625 K. Only slight pressure dependences were found for SiH2 + MeSiH3 ( 485 and 602 K) and for SiH2 + Me2SiH2 ( 600 K). The high pressure rate constants gave the following Arrhenius parameters: [GRAPHICS] These are consistent with fast, near to collision-controlled, association processes. RRKM modelling calculations are consistent with the observed pressure dependences ( and also the lack of them for SiH2 + Me3SiH). Ab initio calculations at both second order perturbation theory (MP2) and coupled cluster (CCSD(T)) levels, showed the presence of weakly-bound complexes along the reaction pathways. In the case of SiH2 + MeSiH3 two complexes, with different geometries, were obtained consistent with earlier studies of SiH2 + SiH4. These complexes were stabilised by methyl substitution in the substrate silane, but all had exceedingly low barriers to rearrangement to product disilanes. Although methyl groups in the substrate silane enhance the intrinsic SiH2 insertion rates, it is doubtful whether the intermediate complexes have a significant effect on the kinetics. A further calculation on the reaction MeSiH + SiH4 shows that the methyl substitution in the silylene should have a much more significant kinetic effect ( as observed in other studies).

Infrared intensities of furan, pyrrole and thiophene: beyond the double harmonic approximation

Infrared intensities of the fundamental, overtone and combination transitions in furan, pyrrole and thiophene have been calculated using the variational normal coordinate code MULTIMODE. We use pure vibrational wavefunctions, and quartic force fields and cubic dipole moment vector surfaces, generated by density functional theory. The results are compared graphically with second-order perturbation calculations and with relative intensities from experiment for furan and pyrrole.

A DFT study of free radicals formed from artemisinin and related compounds

Ab initio calculations using density functional theory have shown that the reactions that occur between artemisinin, 1, a cyclic trioxane active against malaria, and some metal ions and complexes lead to a series of radicals which are probably responsible for its therapeutic activity. In particular it has been shown that the interaction of Fe(H) with artemisinin causes the O-O bond to be broken as indeed does Fe(III) and Cu(I), while Zn(II) does not. Calculations were carried out with Fe(II) in several different forms including the bare ion, [Fe(H2O)(5)](2+) and [FeP(Im)] (P, porphyrin; Im, imadazole) and similar results were obtained. The resulting oxygen-based radicals are readily converted to more stable carbon-based radicals and/or. stable products. Similar radicals and products are also formed from two simple model trioxanes 2 and 3 that show little or no therapeutic action against malaria although some subtle differences were obtained. This suggests that the scaffold surrounding the pharmacophore may be involved in molecular recognition events allowing efficient uptake of this trioxane warhead into the parasite. (C) 2004 Elsevier B.V. All rights reserved.

Synergetic effects of the Cu/Pt{110} surface alloy: enhanced reactivity of water and carbon monoxide

We have used synchrotron-based high-resolution X-ray photoelectron spectroscopy in combination with ab initio density functional theory calculations to investigate the characteristics of water and CO adsorption on the bimetallic Cu/Pt{110}-(2 x 1) surface at a Cu coverage near 0.5 ML. Cu fills the troughs of the reconstructed clean surface forming nanowires, which are stable up to 830 K. Their presence dramatically influences the adsorption of water and CO. Water adsorption changes from intact to partially dissociated while the desorption temperature of CO on this surface increases by up to 27 K with respect to the clean Pt{110} surface. Ab initio calculations and experimental valence band spectra reveal that the Cu 3d-band is narrowed and shifted upward with respect to bulk Cu surfaces. This and electron donation to surface Pt atoms cause the increase in the bond strength between CO and the Pt surface atoms. The pathway for water dissociation occurs via Cu surface atoms. The heat of adsorption of water bonding to Cu surface atoms was calculated to be 0.82 eV, which is significantly higher than on the clean Pt{110} surface; the activation energy for partial dissociation is 0.53 eV (not corrected for zero point energy).

Soluble redox-active polymetallic chains [{Ru-0(CO)(L)(bpy)}(m)](n) (bpy=2,2 '-bipyridine, L = PrCN, Cl-; m=0,-1): electrosynthesis and characterization

Electrochemical and spectroelectrochemical techniques were employed to study in detail the formation and so far unreported spectroscopic properties of soluble electroactive molecular chains with nonbridged metal-metal backbones, namely, [{Ru-0(CO)(PrCN)(bpy)}(m)](n) (m = 0, -1) and [{Ru-0(CO)(bpy)Cl}(m)](n) (m = -1, -2; bpy = 2,2'-bipyridine). The precursors cis-(Cl)-[Ru-II(CO)(MeCN)(bpy)Cl-2] (in PrCN) and mer-[Ru-II(CO)(bpy)Cl-3](-) (in tetrahydrofuran (THF) and PrCN) undergo one-electron reductions to reactive radicals [Ru-II(CO)(MeCN)(bpy(center dot-))Cl-2](-) and [Ru-II(CO)(bpy(center dot-))Cl-3](2-), respectively. Both [bpy(center dot-)]-containing species readily electropolymerize on concomitant dissociation of two chloride ligands and consumption of a second electron. Along this path, mer-to-fac isomerization of the bpy-reduced trichlorido complex (supported by density functional theory calculations) and a concentration-dependent oligomerization process contribute to the complex reactivity pattern. In situ spectroelectrochemistry (IR, UV/vis a has revealed that the charged polymer [{Ru-0(CO)(bpy)Cl}(-)](n) is stable in THF, but in PrCN it converts readily to [Ru-0(CO)(PrCN)(bpy)](n). An excess of chloride ions retards this substitution at low temperatures. Both polymetallic chains are completely soluble in the electrolyte solution and can be reduced reversibly to the corresponding [bpy(center dot-)]-containing species.

Ag(I) and Cu(I) complexes of tetramethyldiphosphinedisulfide: synthesis and structure

Reaction of [M(NCCH3)(4)][PF6] (M = Ag, Cu) with the S2P2Me4 ligand in dichloromethane solution led to substitution of all the nitrile ligands by two molecules of the sulfur ligand, affording the new species [Ag(S2P2Me4)(2)][PF6] (1) and [Cu(S2P2Me4)(2)][PF6] (2). The structures of these complexes were determined by single crystal X-ray diffraction. showing the expected tetrahedral coordination around each metal. Density functional theory (DFT) calculations confirmed the different geometries and energies of the free and coordinated ligand, and provided a very good reproduction of the experimental structures, both for Ag and Cu. The lengths of the S=P bonds are barely affected by coordination, indicating that the pi bond is not important in binding to the metal. (C) 2002 Elsevier Science B.V. All rights reserved.

Intermediate ion stability and regio selectivity polymerization using neutral salicyladiminato in propene nickel(II) and palladium(II) complexes as catalysts

The mechanism of the Heck reaction has been studied with regard to transition metal catalysis of the addition of propene and the formation of unsaturated polymers. The reactivity of nickel and palladium complexes with five different bidentate ligands with O,N donor atoms has been investigated by computational methods involving density functional theory. Hence, it is possible to understand the electronic and steric factors affecting the reaction and their relative importance in determining the products formed in regard of their control of the regiochemistry of the products. Our results show that whether the initial addition of propene is trans to O or to N of the bidentate ligand is of crucial importance to the subsequent reactions. Thus when the propene is trans to 0, 1,2-insertion is favoured, but when the propene is trans to N, then 2,1-insertion is favoured. This difference in the preferred insertion pathway can be related to the charge distribution engendered in the propene moiety when the complex is formed. Indeed charge effects are important for catalytic activity but also for regioselectivity. Steric effects are shown to be of lesser importance even when t-butyl is introduced into the bidentate ligand as a substituent. (C) 2007 Elsevier B.V. All rights reserved.

The crucial importance of agostic interactions in intermediates formed in propylene polymerization using neutral salicyladiminato palladium(II) and nickel complexes as catalysts

The structures of intermediates formed in propylene polymerisation using neutral salicyladiminato palladium(II) and nickel(II) complexes as catalysts have been investigated using density functional theory. Calculations show that all low energy intermediates contain agostic interactions either between the metal and a hydrogen from the added propylene forming four- or five-membered chelate rings, or, when a phenyl ring is present, between the metal and an aromatic C-C bond. The agostic interactions with the metal are concomitant with changes in ligand dimensions and electronic properties. In particular when a metal to hydrogen bond is formed, there is a lengthening of the C-H bond. Significant differences are found for the agostic interactions with palladium and nickel in that for Pd there is a clear preference for specific intermediates but for Ni there are several different structures with similar energies which are likely to lead to a greater variety of products on further polymerisation. (c) 2007 Elsevier B.V. All rights reserved.

A DFT study of propylene polymerization using neutral salicyladiminato nickel(II) and palladium(II) as catalysts

Propylene polymerization using salicyladiminato metal catalalysts has been studied using density functional theory at the B3LYP/LANL2DZ level. In particular, the effects on the reaction mechanisms of changing the metal from Pd(II) to Ni(II) have been investigated. While the reaction mechanisms involving the salicyladiminato Ni(II) catalyst have been found to be similar to those established previously for the salicyladiminato Pd(II) catalyst, the nickel catalyst was found to differentiate the trans-O intermediate from the trans-.N intermediate with an energy difference of 46.63 U mol(-1) significantly more than the palladium catalyst for which the energy difference was calculated as 35.82 kJ mol(-1). The energy difference between the trans-O configuration and the trans-N configuration is decreased significantly when combining a molecule of propylene with the catalyst. For the Ni catalyst, the trans-O isomer is more stable than the trans-N isomer to a greater extent than for Pd, so that the insertion of propylene from 20 is relatively less favoured for Ni than for Pd. It is predicted that the mechanism of isomerization from 20 to 2N through a rotational transition state TS2O2N is more appropriate for the Ni catalyst system. The palladium system shows a larger preference for pi-coordination than its nickel counterpart, although the latter possesses a lower reaction barrier. It was found that the occupation of the trans-O position in the asymmetric salicyladiminato catalyst is also more favored by the alkene as it is by the alkyl so that insertion of the alkene may always start from a particular configuration so that specific products are obtained. (c) 2005 Elsevier B.V. All rights reserved.

Inelastic neutron scattering in spectroscopic studies of hydrogen on carbon-supported catalysts-experimental spectra and computed spectra of model systems

Inelastic neutron scattering spectroscopy has been used to observe and characterise hydrogen on the carbon component of a Pt/C catalyst. INS provides the complete vibration spectrum of coronene, regarded as a molecular model of a graphite layer. The vibrational modes are assigned with the aid of ab initio density functional theory calculations and the INS spectra by the a-CLIMAX program. A spectrum for which the H modes of coronene have been computationally suppressed, a carbon-only coronene spectrum, is a better representation of the spectrum of a graphite layer than is coronene itself. Dihydrogen dosing of a Pt/C catalyst caused amplification of the surface modes of carbon, an effect described as H riding on carbon. From the enhancement of the low energy carbon modes (100-600 cm(-1)) it is concluded that spillover hydrogen becomes attached to dangling bonds at the edges of graphitic regions of the carbon support. (C) 2003 Elsevier Science B.V. All rights reserved.

Charge transfer in Cr adsorption and reaction at the rutile TiO2(110) surface

The rutile TiO2(110) surface has been doped with sub-monolayer metallic Cr, which oxidises and donates charge to specific surface Ti ions. X-Ray and ultra violet photoemission spectroscopy and first principles density functional theory with Hubbard U are used to assign the oxidation states of Cr and surface Ti and we find that Cr2+ forms on bridging oxygen ions and a 5-fold coordinated surface Ti atom is reduced to Ti3+ and the Cr ions readily react with oxygen (to Cr3+), which leads to depletion of surface Ti3+ 3d electrons.