Computational definition of a mixed valent Fe(II)Fe(I) model of the [FeFe]hydrogenase active site resting state

Computational method on highly efficient D-π-A-π-D-based different molecular acceptors for organic solar cells applications and non-linear optical behaviour

Eight molecular structures (BT-A1 to BT-A8) with high-performance non-fullerene acceptor (NFA) were selected for organic solar cells (OSCs) and non-linear optical (NLO) applications. Their electronic, photovoltaic (PV) and optoelectronic properties were tuned by adding powerful electron-withdrawing groups to the acceptor (A) of the D-π-A-π-D structure. Using time-dependent density functional theory (TD-DFT) techniques, based on the laws of quantum chemical calculations, the absorption spectra, stability of the highest and lowest-energy molecular orbitals (HOMO/LUMOs), electron density, intramolecular charge transfer (ICT), transition density matrix (TDM), were examined. The binding energy (Eb) and density of states (DOS) were probed to realize the optoelectronic analysis of the structures BT-A1 to BT-A8. Noncovalent interactions (NCIs) based on a reduced density gradient (RDG) were used to describe the nature and strength of D-A interactions in the molecules BT-A1 to BT-A8. The new refined molecules BT-A1 to BT-A8 exhibited strong absorbance bands between 408-721 nm and high electron transfer contribution (ETC) ranges between 87-96 %, along with the smallest excitation energies (Ex) between 1.71-3.55 eV in the solvent dichloromethane. Dipolar moment strengths ranging from 0.38 to 4.72 Debye in both the excited and ground states have determined with good solubility properties of BT-A1 to BT-A8 in polar solvent. Highly effective charge mobilities and prevention of charge recombination have been demonstrated by the electron (0.18-0.41 eV) and hole RE values (0.13-0.89 eV) for the new compounds. Power conversion efficiencies (PCE) of BT-A1 to BT-A8 were nearly the same because of better outcomes compared to the molecules in the BT. Compared to poly[4.8-bis(5-(2-ethylhexyl)thiophen-2-yl)benzo[1,2-b: 4,5-b']dithiophene-2,6- diyl-alt-(4-2-ethylhexyl)-3-fluorothieno[3,4-b]thiophene-)-2-carboxylate-2-6-diyl)] (PTB7-Th), the open circuit voltages (Voc) of compounds BT-A1 to BT-A8 were ranged from 1.52 to 2.13 eV. The polarizability (α) and hyperpolarizability (β) of the molecules BT-A1 to BT-A8 were used to determine the non-linear optical (NLO) properties. The results showed that BT-A2, BT-A6 and BT-A7 have good NLO activity. This computational analysis demonstrates the superiority of the molecules with NFA. Hence the compounds are advised for the use in production of high-performance OSCs and NLO activity.

[FeFe]-Hydrogenases: maturation and reactivity of enzymatic systems and overview of biomimetic models

[FeFe]-hydrogenases recieved increasing interest in the last decades. This review summarises important findings regarding their enzymatic reactivity as well as inorganic models applied as electro- and photochemical catalysts.

Catalytic H2 evolution/oxidation in [FeFe]-hydrogenase biomimetics: account from DFT on the interplay of related issues and proposed solutions

A DFT overview on selected issues regarding diiron catalysts related to [FeFe]-hydrogenase biomimetic research, with implications for both energy conversion and storage strategies.

Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities

Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.

Hypothesis: entaticversusecstatic states in metalloproteins

Biological metal sites are optimized for function by conformational properties of the protein macroligand.

Proton-Coupled Electron Transfer: Moving Together and Charging Forward

Proton-coupled electron transfer (PCET) is ubiquitous throughout chemistry and biology. This Perspective discusses recent advances and current challenges in the field of PCET, with an emphasis on the role of theory and computation. The fundamental theoretical concepts are summarized, and expressions for rate constants and kinetic isotope effects are provided. Computational methods for calculating reduction potentials and pKa's for molecular electrocatalysts, as well as insights into linear correlations and non-innocent ligands, are also described. In addition, computational methods for simulating the nonadiabatic dynamics of photoexcited PCET are discussed. Representative applications to PCET in solution, proteins, electrochemistry, and photoinduced processes are presented, highlighting the interplay between theoretical and experimental studies. The current challenges and suggested future directions are outlined for each type of application, concluding with an overall view to the future.

Ligand effects on the electrochemical behavior of [Fe2(CO)5(L){μ-(SCH2)2(Ph)PO}] (L = PPh3, P(OEt)3) hydrogenase model complexes

In this paper we study the influence of substituting one CO ligand in [Fe₂(CO)₆{μ-(SCH₂)₂(Ph)PO}] (1) by better σ-donors (PPh₃(2) and P(OMe)₃(3)) in relation to the electrochemical behavior.

Computational investigation of [FeFe]-hydrogenase models: characterization of singly and doubly protonated intermediates and mechanistic insights

The [FeFe]-hydrogenase enzymes catalyze hydrogen oxidation and production efficiently with binuclear Fe metal centers. Recently the bioinspired H₂-producing model system Fe₂(adt)(CO)₂(dppv)₂ (adt=azadithiolate and dppv=diphosphine) was synthesized and studied experimentally. In this system, the azadithiolate bridge facilitates the formation of a doubly protonated ammonium-hydride species through a proton relay. Herein computational methods are utilized to examine this system in the various oxidation states and protonation states along proposed mechanistic pathways for H₂ production. The calculated results agree well with the experimental data for the geometries, CO vibrational stretching frequencies, and reduction potentials. The calculations illustrate that the NH···HFe dihydrogen bonding distance in the doubly protonated species is highly sensitive to the effects of ion-pairing between the ammonium and BF₄^– counterions, which are present in the crystal structure, in that the inclusion of BF₄^– counterions leads to a significantly longer dihydrogen bond. The non-hydride Fe center was found to be the site of reduction for terminal hydride species and unsymmetric bridging hydride species, whereas the reduced symmetric bridging hydride species exhibited spin delocalization between the Fe centers. According to both experimental measurements and theoretical calculations of the relative pK_a values, the Fe_d center of the neutral species is more basic than the amine, and the bridging hydride species is more thermodynamically stable than the terminal hydride species. The calculations implicate a possible pathway for H₂ evolution that involves an intermediate with H₂ weakly bonded to one Fe, a short H₂ distance similar to the molecular bond length, the spin density delocalized over the two Fe centers, and a nearly symmetrically bridged CO ligand. Overall, this study illustrates the mechanistic roles of the ammonium-hydride interaction, flexibility of the bridging CO ligand, and intramolecular electron transfer between the Fe centers in the catalytic cycle. Such insights will assist in the design of more effective bioinspired catalysts for H₂ production.

Theoretical calculations in conjunction with supporting experimental data are used to analyze the mechanistic pathway for hydrogen evolution catalyzed by the bioinspired model Fe₂(adt)(CO)₂(dppv)₂. This study elucidates the site of reduction and the pK_a values associated with formation of the singly and doubly protonated species, as well as the roles of the ammonium-hydride interaction, flexibility of the bridging CO ligand, and intramolecular electron transfer between the Fe centers in the catalytic cycle for H₂ production.

Bioinspired Hydrogenase Models: The Mixed-Valence Triiron Complex [Fe3(CO)7(μ-edt)2] and Phosphine Derivatives [Fe3(CO)7-x (PPh3) x (μ-edt)2] (x = 1, 2) and [Fe3(CO)5(κ2-diphosphine)(μ-edt)2] as Proton Reduction Catalysts

The mixed-valence triiron complexes [Fe3(CO)7-x (PPh3) x (μ-edt)2] (x = 0-2; edt = SCH2CH2S) and [Fe3(CO)5(κ2-diphosphine)(μ-edt)2] (diphosphine = dppv, dppe, dppb, dppn) have been prepared and structurally characterized. All adopt an anti arrangement of the dithiolate bridges, and PPh3 substitution occurs at the apical positions of the outer iron atoms, while the diphosphine complexes exist only in the dibasal form in both the solid state and solution. The carbonyl on the central iron atom is semibridging, and this leads to a rotated structure between the bridged diiron center. IR studies reveal that all complexes are inert to protonation by HBF4·Et2O, but addition of acid to the pentacarbonyl complexes results in one-electron oxidation to yield the moderately stable cations [Fe3(CO)5(PPh3)2(μ-edt)2]+ and [Fe3(CO)5(κ2-diphosphine)(μ-edt)2]+, species which also result upon oxidation by [Cp2Fe][PF6]. The electrochemistry of the formally Fe(I)-Fe(II)-Fe(I) complexes has been investigated. Each undergoes a quasi-reversible oxidation, the potential of which is sensitive to phosphine substitution, generally occurring between 0.15 and 0.50 V, although [Fe3(CO)5(PPh3)2(μ-edt)2] is oxidized at -0.05 V. Reduction of all complexes is irreversible and is again sensitive to phosphine substitution, varying between -1.47 V for [Fe3(CO)7(μ-edt)2] and around -1.7 V for phosphine-substituted complexes. In their one-electron-reduced states, all complexes are catalysts for the reduction of protons to hydrogen, the catalytic overpotential being increased upon successive phosphine substitution. In comparison to the diiron complex [Fe2(CO)6(μ-edt)], [Fe3(CO)7(μ-edt)2] catalyzes proton reduction at 0.36 V less negative potentials. Electronic structure calculations have been carried out in order to fully elucidate the nature of the oxidation and reduction processes. In all complexes, the HOMO comprises an iron-iron bonding orbital localized between the two iron atoms not ligated by the semibridging carbonyl, while the LUMO is highly delocalized in nature and is antibonding between both pairs of iron atoms but also contains an antibonding dithiolate interaction.

The influence of a S-to-S bridge in diiron dithiolate models on the oxidation reaction: a mimic of the Hairox state of [FeFe]-hydrogenases

A diiron complex replicates the [FeFe]-hydrogenase Hairox state and the reversible interconversion with an Fe^IFe^I model of the H_red state.

Connecting [NiFe]- and [FeFe]-hydrogenases: mixed-valence nickel-iron dithiolates with rotated structures

New mixed-valence iron-nickel dithiolates are described that exhibit structures similar to those of mixed-valence diiron dithiolates. The interaction of tricarbonyl salt [(dppe)Ni(pdt)Fe(CO)(3)]BF(4) ([1]BF(4), where dppe = Ph(2)PCH(2)CH(2)PPh(2) and pdt(2-) = -SCH(2)CH(2)CH(2)S-) with P-donor ligands (L) afforded the substituted derivatives [(dppe)Ni(pdt)Fe(CO)(2)L]BF(4) incorporating L = PHCy(2) ([1a]BF(4)), PPh(NEt(2))(2) ([1b]BF(4)), P(NMe(2))(3) ([1c]BF(4)), P(i-Pr)(3) ([1d]BF(4)), and PCy(3) ([1e]BF(4)). The related precursor [(dcpe)Ni(pdt)Fe(CO)(3)]BF(4) ([2]BF(4), where dcpe = Cy(2)PCH(2)CH(2)PCy(2)) gave the more electron-rich family of compounds [(dcpe)Ni(pdt)Fe(CO)(2)L]BF(4) for L = PPh(2)(2-pyridyl) ([2a]BF(4)), PPh(3) ([2b]BF(4)), and PCy(3) ([2c]BF(4)). For bulky and strongly basic monophosphorus ligands, the salts feature distorted coordination geometries at iron: crystallographic analyses of [1e]BF(4) and [2c]BF(4) showed that they adopt "rotated" Fe(I) centers, in which PCy(3) occupies a basal site and one CO ligand partially bridges the Ni and Fe centers. Like the undistorted mixed-valence derivatives, members of the new class of complexes are described as Ni(II)Fe(I) (S = ½) systems according to electron paramagnetic resonance spectroscopy, although with attenuated (31)P hyperfine interactions. Density functional theory calculations using the BP86, B3LYP, and PBE0 exchange-correlation functionals agree with the structural and spectroscopic data, suggesting that the spin for [1e](+) is mostly localized in a Fe(I)-centered d(z(2)) orbital, orthogonal to the Fe-P bond. The PCy(3) complexes, rare examples of species featuring "rotated" Fe centers, both structurally and spectroscopically incorporate features from homobimetallic mixed-valence diiron dithiolates. Also, when the NiS(2)Fe core of the [NiFe]-hydrogenase active site is reproduced, the "hybrid models" incorporate key features of the two major classes of hydrogenase. Furthermore, cyclic voltammetry experiments suggest that the highly basic phosphine ligands enable a second oxidation corresponding to the couple [(dxpe)Ni(pdt)Fe(CO)(2)L](+/2+). The resulting unsaturated 32e(-) dications represent the closest approach to modeling the highly electrophilic Ni-SI(a) state. In the case of L = PPh(2) (2-pyridyl), chelation of this ligand accompanies the second oxidation.

A unifying structural and electronic concept for Hmd and [FeFe] hydrogenase active sites

The hydrogenases [FeFe] and Hmd feature at first sight rather different active sites. A closer inspection reveals striking similarities, which allow us to define swapped ligand spheres in such a way that the single active iron center of Hmd functions in a first-shell ligand environment resembling the reacting iron atom in [FeFe] hydrogenase and vice versa. These redesigned ligand environments can be conveniently studied with quantum chemical methods and point to general reactivity principles for iron centers with hydrogenase activity.

Redox and structural properties of mixed-valence models for the active site of the [FeFe]-hydrogenase: progress and challenges

The one-electron oxidations of a series of diiron(I) dithiolato carbonyls were examined to evaluate the factors that affect the oxidation state assignments, structures, and reactivity of these low-molecular weight models for the H ox state of the [FeFe]-hydrogenases. The propanedithiolates Fe 2(S 2C 3H 6)(CO) 3(L)(dppv) (L = CO, PMe 3, P i-Pr 3) oxidize at potentials approximately 180 mV milder than the related ethanedithiolates ( Angew. Chem., Int. Ed. 2007, 46, 6152). The steric clash between the central methylene of the propanedithiolate and the phosphine favors the rotated structure, which forms upon oxidation. Electron Paramagnetic Resonance (EPR) spectra for the mixed-valence cations indicate that the unpaired electron is localized on the Fe(CO)(dppv) center in both [Fe 2(S 2C 3H 6)(CO) 4(dppv)]BF 4 and [Fe 2(S 2C 3H 6)(CO) 3(PMe 3)(dppv)]BF 4, as seen previously for the ethanedithiolate [Fe 2(S 2C 2H 4)(CO) 3(PMe 3)(dppv)]BF 4. For [Fe 2(S 2C n H 2 n )(CO) 3(P i-Pr 3)(dppv)]BF 4; however, the spin is localized on the Fe(CO) 2(P i-Pr 3) center, although the Fe(CO)(dppv) site is rotated in the crystalline state. IR and EPR spectra, as well as redox potentials and density-functional theory (DFT) calculations, suggest that the Fe(CO) 2(P i-Pr 3) site is rotated in solution, driven by steric factors. Analysis of the DFT-computed partial atomic charges for the mixed-valence species shows that the Fe atom featuring a vacant apical coordination position is an electrophilic Fe(I) center. One-electron oxidation of [Fe 2(S 2C 2H 4)(CN)(CO) 3(dppv)] (-) resulted in 2e oxidation of 0.5 equiv to give the mu-cyano derivative [Fe (I) 2(S 2C 2H 4)(CO) 3(dppv)](mu-CN)[Fe (II) 2(S 2C 2H 4)(mu-CO)(CO) 2(CN)(dppv)], which was characterized spectroscopically.

Diiron dithiolato carbonyls related to the H(ox)CO state of [FeFe]-hydrogenase

Oxidation of the electron-rich (E(1/2) = -175 vs Ag/AgCl) ethanedithiolato complex Fe2(S2C2H4)(CO)2(dppv)2 (1) under a CO atmosphere yielded [Fe2(S2C2H4)(mu-CO)(CO)2(dppv)2](+) ([1(CO)](+)), a model for the H(ox)(CO) state of the [FeFe]-hydrogenases. This complex exists as two isomers: a kinetically favored unsymmetrical derivative, unsym-[1(CO)](+), and a thermodynamically favored isomer, sym-[1(CO)](+), wherein both diphosphines span apical and basal sites. Crystallographic characterization of sym-[1(CO)](+) confirmed a C2-symmetric structure with a bridging CO ligand and an elongated Fe-Fe bond of 2.7012(14) A, as predicted previously. Oxidation of sym-[1(CO)](+) and unsym-[1(CO)](+) again by 1e(-) oxidation afforded the respective diamagnetic diferrous derivatives where the relative stabilities of the sym and unsym isomers are reversed. DFT calculations indicate that the stabilities of sym and unsym isomers are affected differently by the oxidation state of the diiron unit: the mutually trans CO ligands in the sym isomer are more destabilizing in the mixed-valence state than in the diferrous state. EPR analysis of mixed-valence complexes revealed that, for [1](+), the unpaired spin is localized on a single iron center, whereas for unsym/sym-[1(CO)](+), the unpaired spin was delocalized over both iron centers, as indicated by the magnitude of the hyperfine coupling to the phosphine ligands trans to the Fe-Fe vector. Oxidation of 1 by 2 equiv of acetylferrocenium afforded the dication [1](2+), which, on the basis of low-temperature IR spectrum, is structurally similar to [1](+). Treatment of [1](2+) with CO gives unsym-[1(CO)](2+).