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

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.

Electronic structure of the nitrogen-vacancy center in diamond from first-principles theory

The nitrogen-vacancy (NV) center is a paramagnetic defect in diamond with applications as a qubit. Here, we investigate its electronic structure by using ab initio density functional theory for five different NV center models of two different cluster sizes. We describe the symmetry and energetics of the low-lying states and compare the optical frequencies obtained to experimental results. We compute the major transition of the negatively charged NV centers to within 25–100 meV accuracy and find that it is energetically favorable for substitutional nitrogens to donate an electron to NV0. The excited state of the major transition and the NV0 state with a neutral donor nitrogen are found to be close in energy.

Density-potential mapping in time-dependent density-functional theory

The key questions of uniqueness and existence in time-dependent density-functional theory are usually formulated only for potentials and densities that are analytic in time. Simple examples, standard in quantum mechanics, lead, however, to nonanalyticities. We reformulate these questions in terms of a nonlinear Schroedinger equation with a potential that depends nonlocally on the wave function.

A complete analytic theory for structure and dynamics of populations and communities spanning wide ranges in body size

The prediction and management of ecosystem responses to global environmental change would profit from a clearer understanding of the mechanisms determining the structure and dynamics of ecological communities. The analytic theory presented here develops a causally closed picture for the mechanisms controlling community and population size structure, in particular community size spectra, and their dynamic responses to perturbations, with emphasis on marine ecosystems. Important implications are summarised in non-technical form. These include the identification of three different responses of community size spectra to size-specific pressures (of which one is the classical trophic cascade), an explanation for the observed slow recovery of fish communities from exploitation, and clarification of the mechanism controlling predation mortality rates. The theory builds on a community model that describes trophic interactions among size-structured populations and explicitly represents the full life cycles of species. An approximate time-dependent analytic solution of the model is obtained by coarse graining over maturation body sizes to obtain a simple description of the model steady state, linearising near the steady state, and then eliminating intraspecific size structure by means of the quasi-neutral approximation. The result is a convolution equation for trophic interactions among species of different maturation body sizes, which is solved analytically using a novel technique based on a multiscale expansion.

The electronic band structure of indium selenide:Photoemission and theory

The electronic band structure of vacuum cleaved single-crystal indium selenide has been investigated by X-ray and ultraviolet photoelectron spectroscopy. The valence band consists of three well separated groups, one derived from the Se 4s levels, and two derived from p-like wavefunctions. The band structure and valence band density of states has been calculated using a tight-binding single-layer approximation and all the major features in the experimental spectra are well accounted for. The spin-orbit splitting and electron loss structure associated with the In 4d core level is also reported.

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.

An insight into alkali promotion: a density functional theory study of CO dissociation on K/Rh(111)

The important role of alkali additives in heterogeneous catalysis is, to a large extent, related to the high promotion effect they have on many fundamental reactions. The wide application of alkali additives in industry does not, however, reflect a thorough understanding of the mechanism of their promotional abilities. To investigate the physical origin of the alkali promotion effect, we have studied CO dissociation on clean Rh(111) and K-covered Rh(111) surfaces using density functional theory. By varying the position of potassium atoms relative to a dissociating CO, we have mapped out the importance of different K effects on the CO dissociation reactions. The K-induced changes in the reaction pathways and reaction barriers have been determined; in particular, a large reduction of the CO dissociation barrier has been identified. A thorough analysis of this promotion effect allows us to rationalize both the electronic and the geometrical factors that govern alkali promotion effect: (i) The extent of barrier reductions depends strongly on how close K is to the dissociating CO. (ii) Direct K-O bonding that is in a very short range plays a crucial role in reducing the barrier. (iii) K can have a rather long-range effect on the TS structure, which could reduce slightly the barriers.

Bandstructure approach to near edge structure

We review the current state of the art in EELS fingerprinting by computer simulation, focusing on the bandstructure approach to the problem. Currently calculations are made using a one electron theory, but we describe in principle the way to go beyond this to include final state effects. We include these effects within the one electron framework using the Slater transition state formula and assess the errors involved. Two examples are then given which illustrate the use of the one electron approximation within density functional theory. Our approach is to combine predicted atomic structure with predicted electronic structure to assist in fingerprinting of complex crystal structures.

Effect of relaxation on the oxygen K-edge electron energy-loss near-edge structure in yttria-stabilized zirconia

The electron energy-loss near-edge structure (ELNES) at the oxygen K-edge has been investigated in a range of yttria-stabilized zirconia (YSZ) materials. The electronic structure of the three polymorphs of pure ZrO2 and of the doped YSZ structure close to the 33 mol %Y2O3 composition have been calculated using a full-potential linear muffin-tin orbital method (NFP-LMTO) as well as a pseudopotential based technique. Calculations of the ELNES dipole transition matrix elements in the framework of the NFP-LMTO scheme and inclusion of core hole screening within Slater's transition state theory enable the ELNES to be computed. Good agreement between the experimental and calculated ELNES is obtained for pure monoclinic ZrO2. The agreement is less good with the ideal tetragonal and cubic structures. This is because the inclusion of defects is essential in the calculation of the YSZ ELNES. If the model used contains ordered defects such as vacancies and metal Y planes, agreement between the calculated and experimental O K-edges is significantly improved. The calculations show how the five different O environments of Zr,Y,O, are connected with the features observed in the experimental spectra and demonstrate clearly the power of using ELNES to probe the stabilization mechanism in doped metal oxides.

The near-edge structure in energy-loss spectroscopy: many-electron and magnetic effects in transition metal nitrides and carbides

We investigate the ability of the local density approximation (LDA) in density functional theory to predict the near-edge structure in electron energy-loss spectroscopy in the dipole approximation. We include screening of the core hole within the LDA using Slater's transition state theory. We find that anion K-edge threshold energies are systematically overestimated by 4.22 +/- 0.44 eV in twelve transition metal carbides and nitrides in the rock-salt (B1) structure. When we apply this 'universal' many-electron correction to energy-loss spectra calculated within the transition state approximation to LDA, we find quantitative agreement with experiment to within one or two eV for TiC, TiN and VN. We compare our calculations to a simpler approach using a projected Mulliken density which honours the dipole selection rule, in place of the dipole matrix element itself. We find remarkably close agreement between these two approaches. Finally, we show an anomaly in the near-edge structure in CrN to be due to magnetic structure. In particular, we find that the N K edge in fact probes the magnetic moments and alignments of ther sublattice.

Simplified theory of wideband spectra of liquid H2O and D2O (from 0 to 1000 cm(-1)) due to reorienting polar and vibrating H-bonded water molecules

A semi-phenomenological molecular model is presented, which is capable of describing with the use of analytical formulae, the wideband dielectric(1) and far-infrared spectra of ordinary and heavy water. In the model the vector of a dipole moment is presented as a sum of two components. The absolute value of the first one is constant; the second one changes harmonically with time. The key aspect of this work is consideration of FIR spectra due to the second component. In the context of the modified hybrid model presented in the work, reorientation of the dipoles in the rectangular potential well is considered, as a result of which the librational (near 700 cm (-1)) and translational (near 200 cm (-1)) absorption bands and the microwave Debye relaxation spectrum arise. It is shown that the time-dependent part of a dipole moment contributes most to the translational band, the relevant mechanism is taken to be stretching vibration of the H-bonded molecules. Previous linear-response molecular models were unsuccessful in describing this band (in heavy water) in terms of the complex dielectric permittivity. The spatial and time scales characteristic of water are estimated. (C) 2002 Elsevier Science B.V. All rights reserved.