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

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.

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.

Large-scale density functional calculations of the surface properties of the Wigner crystal

The surface properties of the jellium model have been investigated by large supercell computations in the density functional theory-local spin-density (DFT-LSD) approach for planar slabs with up to 1000 electrons. A wide interval of densities has been explored, extending into the stability range of the Wigner crystal. Most computations have been carried out on nominally paramagnetic samples with an equal number of spin-up and spin-down electrons. The results show that within DFT-LSD spontaneous spin polarization and charge localization start nearly simultaneously at the surface for r(s) similar to 20, then, with decreasing density, they progress toward the center of the slab. Electrons are fully localized and spin polarized at r(s) = 30. At this density the charge distribution is the superposition of disjoint charge blobs, each corresponding to one electron. The distribution of blobs displays both regularities and disorder, the first being represented by well-defined planes and simple in-plane geometries, and the latter by a variety of surface defects. The surface energy, surface dipole, electric polarisability, and magnetization pattern have been determined as a function of density. All these quantities display characteristic anomalies at the density of the localization transition. The analysis of the low-frequency electric conductivity shows that in the fluid paramagnetic regime the in-plane current preferentially flows in the central region of the slab and the two spin channels are equally conducting. In the charge localized, spin-polarized regime, conductivity is primarily a surface effect, and an apparent asymmetry is observed in the two spin currents.

Focusing Dynamics of High-Energy Density, Laser-Driven Ion Beams

The dynamics of the focusing of laser-driven ion beams produced from concave solid targets was studied. Most of the ion beam energy is observed to converge at the center of the cylindrical targets with a spot diameter of 30 mu m, which can be very beneficial for applications requiring high beam energy densities. Also, unbalanced laser irradiation does not compromise the focusability of the beam. However, significant filamentation occurs during the focusing, potentially limiting the localization of the energy deposition region by these beams at focus. These effects could impact the applicability of such high-energy density beams for applications, e. g., in proton-driven fast ignition.

The failure of generalized gradient approximations (GGAs) and meta-GGAs for the two-center three-electron bonds in He-2(+), (H2O)(2)(+), and (NH3)(2)(+)

The radical cations He-2(+) (H2O)(2)(+), and (NH3)(2)(+) with two-center three-electron A-A bonds are investigated at the configuration interaction (CI), accurate Kohn-Sham (KS), generalized gradient approximation (GGA), and meta-GGA levels. Assessment of seven different GGA and six meta-GGA methods shows that the A(2)(+) systems remain a difficult case for density functional theory (DFT). All methods tested consistently overestimate the stability of A(2)(+): the corresponding D-e errors decrease for more diffuse valence densities in the series He-2(+) > (H2O)(2)(+) > (NH3)(2)(+). Upon comparison to the energy terms of the accurate Kohn-Sham solutions, the approximate exchange functionals are found to be responsible for the errors of GGA-type methods, which characteristically overestimate the exchange in A(2)(+). These so-called exchange functionals implicitly use localized holes. Such localized holes do occur if there is left-right correlation, i.e., the exchange functionals then also describe nondynamical correlation. However, in the hemibonded A(2)(+) systems the typical molecular (left-right, nondynamical) correlation of the two-electron pair bond is absent. The nondynamical correlation built into the exchange functionals is then spurious and yields too low energies.

Hyperfine interaction in the ground state of the negatively charged nitrogen vacancy center in diamond

The N-14, N-15, and C-13 hyperfine interactions in the ground state of the negatively charged nitrogen vacancy (NV-) center have been investigated using electron-paramagnetic-resonance spectroscopy. The previously published parameters for the N-14 hyperfine interaction do not produce a satisfactory fit to the experimental NV- electron-paramagnetic-resonance data. The small anisotropic component of the NV- hyperfine interaction can be explained from dipolar interaction between the nitrogen nucleus and the unpaired-electron probability density localized on the three carbon atoms neighboring the vacancy. Optical spin polarization of the NV- ground state was used to enhance the electron-paramagnetic-resonance sensitivity enabling detailed study of the hyperfine interaction with C-13 neighbors. The data confirmed the identification of three equivalent carbon nearest neighbors but indicated the next largest C-13 interaction is with six, rather than as previously assumed three, equivalent neighboring carbon atoms.

Density functional theory study on the activation of molecular oxygen on a stepped gold surface in an aqueous environment: a new approach for simulating reactions in solution

The activation of oxygen molecules is an important issue in the gold-catalyzed partial oxidation of alcohols in aqueous solution. The complexity of the solution arising from a large number of solvent molecules makes it difficult to study the reaction in the system. In this work, O-2 activation on an Au catalyst is investigated using an effective approach to estimate the reaction barriers in the presence of solvent. Our calculations show that O-2 can be activated, undergoing OOH* in the presence of water molecules. The OOH* can readily be formed on Au(211) via four possible pathways with almost equivalent free energy barriers at the aqueous-solid interface: the direct or indirect activation of O-2 by surface hydrogen or the hydrolysis of O-2 following a Langmuir-Hinshelwood mechanism or an Eley-Rideal mechanism. Among them, the Eley-Rideal mechanism may be slightly more favorable due to the restriction of the low coverage of surface H on Au(211) in the other mechanisms. The results shed light on the importance of water molecules on the activation of oxygen in gold-catalyzed systems. Solvent is found to facilitate the oxygen activation process mainly by offering extra electrons and stabilizing the transition states. A correlation between the energy barrier and the negative charge of the reaction center is found. The activation barrier is substantially reduced by the aqueous environment, in which the first solvation shell plays the most important role in the barrier reduction. Our approach may be useful for estimating the reaction barriers in aqueous systems.