Calculations of One-Electron Redox Potentials of Oxoiron(IV) Porphyrin Complexes

Mechanism of Methyl Transfer Reaction between CH3Co(dmgBF2)2py and PPh3Ni(Triphos)

DFT calculations were performed for the methyl group transfer reaction between CH3Co (dmgBF2)py and PPh3Ni(Triphos). The reaction mechanism and its energetics were investigated. This reaction is relevant to the catalytic mechanism of the enzyme acetyl coenzyme A synthase. BP86 and PBE functionals and dispersion corrections were used. It was found that intermolecular interactions are very important for this reaction. The influence of the solvent on the reaction was studied.

Electrochemical generation of high-valent oxo-manganese complexes featuring an anionic N5 ligand and their role in O-O bond formation

Generation of high-valent oxomanganese complexes through controlled removal of protons and electrons from low-valent congeners is a crucial step toward the synthesis of functional analogues of the native oxygen evolving complex (OEC). In-depth studies of the water oxidation activity of such biomimetic compounds help in understanding the mechanism of O-O bond formation presumably occurring in the last step of the photosynthetic cycle. Scarce reports of reactive high-valent oxomanganese complexes underscore the impetus for the present work, wherein we report the electrochemical generation of the non-heme oxomanganese(IV) species [(dpaq)MnIV(O)]+ (2) through a proton-coupled electron transfer (PCET) process from the hydroxomanganese complex [(dpaq)MnIII(OH)]ClO4 (1). Controlled potential spectroelectrochemical studies of 1 in wet acetonitrile at 1.45 V vs. NHE revealed quantitative formation of 2 within 10 min. The high-valent oxomanganese(IV) transient exhibited remarkable stability and could be reverted to the starting complex (1) by switching the potential to 0.25 V vs. NHE. The formation of 2via PCET oxidation of 1 demonstrates an alternate pathway for the generation of the oxomanganese(IV) transient (2) without the requirement of redox-inactive metal ions or acid additives as proposed earlier. Theoretical studies predict that one-electron oxidation of [(dpaq)MnIV(O)]+ (2) forms a manganese(V)-oxo (3) species, which can be oxidized further by one electron to a formal manganese(VI)-oxo transient (4). Theoretical analyses suggest that the first oxidation event (2 to 3) takes place at the metal-based d-orbital, whereas, in the second oxidation process (3 to 4), the electron eliminates from an orbital composed of equitable contribution from the metal and the ligand, leaving a single electron in the quinoline-dominant orbital in the doublet ground spin state of the manganese(VI)-oxo species (4). This mixed metal-ligand (quinoline)-based oxidation is proposed to generate a formal Mn(VI) species (4), a non-heme analogue of the species 'compound I', formed in the catalytic cycle of cytochrome P-450. We propose that the highly electrophilic species 4 catches water during cyclic voltammetry experiments and results in O-O bond formation leading to electrocatalytic oxidation of water to hydrogen peroxide.

Sterically attenuated electronic communication in cobalt complexes of meridional isoquinoline-derived ligands for applications in electrocatalysis

The ability to synthetically tune the ligand frameworks of redox-active molecules is of critical importance to the economy of solar fuels because manipulating their redox properties can afford control over the operating potentials of sustained electrocatalytic or photoelectrocatalytic processes. The electronic and steric properties of 2,2':6',2″-terpyridine (Terpy) ligand frameworks can be tuned by functional group substitution on ligand backbones, and these correlate strongly to their Hammett parameters. The synthesis of a new series of tridentate meridional ligands of 2,4,6-trisubstituted pyridines that engineers the ability to finely tune the redox potentials of cobalt complexes to more positive potentials than that of their Terpy analogs is achieved by aryl-functionalizing at the four-position and by including isoquinoline at the two- and six-positions of pyridine (Aryl-DiQ). Their cobalt complex syntheses, their electronic properties, and their catalytic activity for carbon dioxide (CO2) reduction are reported and compared to their Terpy analogs. The cobalt derivatives generally experience a positive shift in their redox features relative to the Terpy-based analogs, covering a complementary potential range. Although those evaluated fail to produce any quantifiable products for the reduction of CO2 and suffer from long-term instability, these results suggest possible alternate strategies for stabilizing these compounds during catalysis. We speculate that lower equilibrium association constants to the cobalt center are intrinsic to these ligands, which originate from a steric interaction between protons on the pyridine and isoquinoline moieties. Nevertheless, the new Aryl-DiQ ligand framework has been engineered to selectively tune homoleptic cobalt complexes' redox potentials.

Factors controlling the reactivity of synthetic compound-I Analogs

A high-valent iron(IV)-oxo porphyrin radical cation (FeIV(O)(porph+•) serves as a key, reactive intermediate for a range of heme enzymes, including cytochrome P450 (CYP), horseradish peroxidase (HRP), and catalase (CAT). Synthetic analogs of this intermediate, known as Compound-I (Cpd-I) in the heme enzyme literature, have been generated with different tetrapyrrolic, macrocyclic ligands, including porphyrin derivatives, and the closely related ring-contracted macrocycles, corroles and corrolazines. These synthetic analogs have been useful for assigning and understanding structural and spectroscopic features and examining the reactivity of Cpd-I-like species in controlled and well-defined environments. This review focuses on summarizing recent developments in the synthesis and reactivity of high-valent iron-oxo porphyrinoid complexes in two main classes of reactions, proton-coupled electron transfer (PCET) and oxygen atom transfer (OAT). The relationship between the structure of the complexes and their reactivity is emphasized, including the influence of axial ligation and peripheral macrocyclic substitution, as well as the effects of solvent and secondary coordination spheres on the reactivity of the Cpd-I analogs. In bringing together the latest findings on Cpd-I analogs, this review intends to broaden our current understanding of the factors that control the stability and reactivity of Cpd-I species. This new knowledge should, in turn, point toward new synthetic strategies for constructing catalysts that rely on Cpd-I-like reactive intermediates.

Ferrocene/ferrocenium, cobaltocene/cobaltocenium and nickelocene/nickelocenium: from gas phase ionization energy to one-electron reduction potential in solvated medium

We propose a theoretical procedure for accurate determination of reduction potentials for three metallocene couples, Cp₂M⁺/Cp₂M, where M = Fe, Co and Ni. This procedure first computes the gas phase ionization energy (IE) using the explicitly correlated CCSD(T)-F12 method and includes the zero-point energy correction, core-valence electronic correlation, and relativistic and spin-orbit coupling effects. By means of Born-Haber thermochemical cycle, the one-electron reduction potential is obtained as the sum of the gas phase IE and the corresponding Gibbs free energies of solvation (ΔG_solv) for both the neutral and cationic species. Among the three solvent models (PCM, SMD and uESE) investigated here, it turns out that only the SMD model (computed at the DFT level) gives the best estimation of the value for "ΔG_solv(cation) - ΔG_solv(neutral)" and thus, combining with the accurate IE values, the theoretical protocol is capable of yielding reliable values (in V) for , and . These predictions compare favorably with the available experimental data (in V): , , and . We show that our theoretical procedure is reliable for accurate reduction potential predictions of Cp₂Fe⁺/Cp₂Fe, Cp₂Co⁺/Cp₂Co and Cp₂Ni⁺/Cp₂Ni redox couples in aqueous and non-aqueous media; the maximum absolute deviation is as small as ≈120 mV, which outperforms those of the existing theoretical methods.

Geometrical benchmarking and analysis of redox potentials of copper(I/II) guanidine-quinoline complexes: Comparison of semi-empirical tight-binding and DFT methods and the challenge of describing the entatic state (part III)

Copper guanidine-quinoline complexes are an important class of bioinorganic complexes that find utilization in electron and atom transfer processes. By substitution of functional groups on the quinoline moiety the electron transfer abilities of these complexes can be tuned. In order to explore the full substitution space by simulations, the accurate theoretical description of the effect of functional groups is essential. In this study, we compare three different methods for the theoretical description of the structures. We use the semi-empirical tight-binding method GFN2-xTB, the density functional TPSSh and the double-hybrid functional B2PLYP. We evaluate the methods on five different complex pairs (Cu(I) and Cu(II) complexes), and compare how well calculated energies can predict the redox potentials. We find even though B2PLYP and TPSSh yield better accordance with the experimental structures. GFN2-xTB performs surprisingly well in the geometry optimization at a fraction of the computational cost. TPSSh offers a good compromise between computational cost and accuracy of the redox potential for real-life complexes.

Dimerization of [FeIII(bpy)3]3+ in Aqueous Solutions: Elucidating a Mechanism Based on Historical Proposals, Electrochemical Data, and Computational Free Energy Analysis

Iron(II) tris-bipyridine, [Fe^II(bpy)₃]²⁺, is a historically significant organometallic coordination complex with attractive redox and photophysical properties. With respect to energy storage, it is a low-cost, high-redox potential complex and thus attractive for use as a catholyte in aqueous redox flow batteries. Despite these favorable characteristics, its oxidized Fe(III) form undergoes dimerization to form μ-O-[Fe^III(bpy)₂(H₂O)]₂⁴⁺, leading to a dramatic ∼0.7 V decrease during battery discharge. To date, the energetics and complete mechanism of this slow, sequential electrochemical-chemical (EC) process, which includes electron transfer, nucleophilic attack, ligand cleavage, μ-oxo bond formation, and spin state transition, have not been elucidated. Using cyclic voltammetry, redox flow battery data, and density functional theory calculations guided by previously proposed mechanisms, we modeled more than 100 complexes and performed more than 50 geometry scans to resolve the key steps dictating these complex chemical processes. Quantitative free energy surfaces are developed to model the mechanism of dimerization accounting for the spins and identities of any possible Fe(II), Fe(III), or Fe(IV) intermediates. Electrochemical reduction of the dimer regenerates [Fe^II(bpy)₃]²⁺ in an overall reversible process. Computational electrochemistry interrogates the influence of spin state, coordination environment, and molecular conformation at the electrode-electrolyte interface through a proposed stepwise dimer reduction process. Experimentally, we show that the considerable overpotential associated with this event can be catalytically mitigated with disparate materials, including platinum, copper hexacyanoferrate, and activated carbon. The findings are of fundamental and applied significance and could elevate [Fe^II(bpy)₃]²⁺ and its derivatives to play a vital role in the burgeoning renewable energy economy.

Fullerene binding effects in Al(III)/Zn(II) Porphyrin/Phthalocyanine photophysical properties and charge transport

We evaluate the fullerene C₆₀ binding effect; through the metal (Al) and through the ligand (Pc,TPP), on the photophysical and charge transport properties of M-porphyrin(TPP)/phthalocyanine(Pc) (M = Al(III), Zn(II)). We perform density functional theory (DFT) and time-dependent DFT calculations for the macrocycle-C₆₀ dyads, showing that all systems studied are thermodynamically favorable. The C₆₀ binding effect on the absorption spectrum is a red-shift of the Q and Soret (B) bands of TPPs and Pcs. The Pc-dyads show longer λ for Q bands (673 nm) than those with TPP (568 nm). AlTPP-C₆₀ and ZnTPP-C₆₀ show a more favorable electron injection to TiO₂ than the analogs Pcs, and the regeneration of the dye is preferred in AlTPP-C₆₀ and AlPc-C₆₀. Zero-bias conductance is computed (10^-4-10^-7 G₀) for the dyads using molecular junctions with Au(111)-based electrodes. When a bias voltage of around 0.6 V up to 1 V is applied, an increase in current is obtained for AlTPP-C₆₀ (10^-7 A), ZnTPP-C₆₀ (10^-7 A), and AlPc-C₆₀ (10^-8 A). Although there is not a unique trend in the behavior of the dyads, Pcs have better photophysical properties than TPPs and the latter are better in the charge transport. We conclude that AlTPP(ZnTPP)-C₆₀ dyads are an excellent alternative for designing new materials for dye-sensitized solar cells or optoelectronic devices.

Theoretical Studies of Acetyl-CoA Synthase Catalytic Mechanism

DFT calculations were performed for the A-cluster from the enzyme Acetyl-CoA synthase (ACS). The acid constants (pKa), reduction potentials, and pH-dependent reduction potential for the A-cluster with different oxidation states and ligands were calculated. Good agreement of the reduction potentials, dependent on pH in the experiment, was obtained. On the basis of the calculations, a mechanism for the methylation reaction involving two–electron reduction and protonation on the proximal nickel atom of the reduced A-cluster is proposed.

Progress on photocatalytic semiconductor hybrids for bacterial inactivation

Due to its use of green and renewable energy and negligible bacterial resistance, photocatalytic bacterial inactivation is to be considered a promising sterilization process. Herein, we explore the relevant mechanisms of the photoinduced process on the active sites of semiconductors with an emphasis on the active sites of semiconductors, the photoexcited electron transfer, ROS-induced toxicity and interactions between semiconductors and bacteria. Pristine semiconductors such as metal oxides (TiO₂ and ZnO) have been widely reported; however, they suffer some drawbacks such as narrow optical response and high photogenerated carrier recombination. Herein, some typical modification strategies will be discussed including noble metal doping, ion doping, hybrid heterojunctions and dye sensitization. Besides, the biosafety and biocompatibility issues of semiconductor materials are also considered for the evaluation of their potential for further biomedical applications. Furthermore, 2D materials have become promising candidates in recent years due to their wide optical response to NIR light, superior antibacterial activity and favorable biocompatibility. Besides, the current research limitations and challenges are illustrated to introduce the appealing directions and design considerations for the future development of photocatalytic semiconductors for antibacterial applications.

Fine-Tuning Metal and Ligand-Centered Redox Potentials of Homoleptic Bis-Terpyridine Complexes with 4'-Aryl Substituents

Homoleptic transition-metal complexes of 2,2':6',2″-terpyridine (terpy) and substituted derivatives of the form [M(R-terpy)₂]²⁺ display a wide range of redox potentials that correlate well to the Hammett parameter of the terpy substituents. Less is known about the impact of incorporating a phenyl spacer between the functional group responsible for controlling the electron density of terpy and how that translates to metal complexes of the form [M(4'-aryl-terpy)₂]²⁺, where M = Mn, Fe, Co, Ni, and Zn. Herein, we report our studies on these complexes revealed a good correlation of redox potentials of both metal- and ligand-centered events with the Hammett parameters of the aryl substituents, regardless of aryl-substitution pattern (i.e., the presence of multiple functional groups, combinations of withdrawing and donating functional groups). The phenyl spacer results in 60-80% attenuation of electron density as compared to the 4'-substituted terpy analogue, depending on the metal and redox couple analyzed. Density functional theory calculations performed on a simple model system revealed a strong correlation between the Hammett parameters and lowest unoccupied molecular orbital energies of the corresponding substituted pyridine models, thus serving as an inexpensive predictive tool when coupled with electrochemical data. Overall, these data suggest that such ligand modifications may be used in combination with previous approaches to further fine-tune the redox potentials of homoleptic transition-metal complexes, which may have applications in photochemical and electrochemical catalytic processes.

Directed Design of a AuI Complex with a Reduced Mesoionic Carbene Radical Ligand: Insights from 1,2,3-Triazolylidene Selenium Adducts and Extensive Electrochemical Investigations

Carbene-based radicals are important for both fundamental and applied chemical research. Herein, extensive electrochemical investigations of nine different 1,2,3-triazolylidene selenium adducts are reported. It is found that the half-wave potentials of the first reduction of the selones correlate with their calculated LUMO levels and the LUMO levels of the corresponding triazolylidene-based mesoionic carbenes (MICs). Furthermore, unexpected quasi-reversibility of the reduction of two triazoline selones, exhibiting comparable reduction potentials, was discovered. Through UV/Vis/NIR and EPR spectroelectrochemical investigations supported by DFT calculations, the radical anion was unambiguously assigned to be triazoline centered. This electrochemical behavior was transferred to a triazolylidene-type MIC-gold phenyl complex resulting in a MIC-radical coordinated AuI species. Apart from UV-Vis-NIR and EPR spectroelectrochemical investigations of the reduction, the reduced gold-coordinated MIC radical complex was also formed in situ in the bulk through chemical reduction. This is the first report of a monodentate triazolylidene-based MIC ligand that can be reduced to its anion radical in a metal complex. The results presented here provide design principles for stabilizing radicals based on MICs.

Ab Initio Evaluation of the Redox Potential of Cytochrome c

Various biochemical activities of metabolism and biosynthesis are fulfilled by redox processes with explicit electron exchange, which furnish redox enzymes with high chemical reactivity. However, theoretical investigation of a redox process, which simultaneously involves a complex electronic change at a redox metal center and conformational reorganization of the surrounding protein environment coupled to the electronic change, requires computationally conflicting approaches, highly accurate quantum chemical calculations, and long-time molecular dynamics (MD) simulations, limiting the physicochemical understanding of biological redox processes. Here, we theoretically examined a redox process of cytochrome c by means of a hybrid molecular simulation technique, which enables one to consistently treat the redox center at the ab initio quantum chemistry level of theory and the protein reorganization with long-time MD simulations on the microsecond timescale. The calculations successfully evaluated a large absolute redox potential, 4.34 eV, with errors of only 0.03 to 0.34 eV to the experimental ones without any problem-specific empirical parameters. Through the long-time MD sampling, large and nonlinear reorganization of the protein environment was unveiled and the molecular determinants for the redox potential were identified. The present ab initio approach significantly expands the applicability of theoretical investigation to biological redox systems with more electronically complicated redox centers such as polynuclear transition metal complexes.

The bimetallic and the anchoring group effects on both optical and charge transport properties of hexaphyrin amethyrin

Bimetallic Cu(ii)-hexaphyrin amethyrin proposed as a molecular switch operated by the application of an external magnetic field.

Oxidative ability of organic iodine(iii) reagents: a theoretical assessment

Computational prediction of the oxidative ability of hypervalent iodine reagents was performed based on redox potentials, calculated at the B3LYP/6-311+G(d,p) level of theory with the SDD ECP46MDF pseudopotential for the iodine atom.

Computational electrochemistry of a novel ferrocene derivative

In this study, the structural and redox properties of a novel ferrocene derivative in dichlomethane solvent were investigated. For this aim, various exchange-correlation functionals and basis sets in gas phase with different continuum solvation models and cavities in liquid phase were applied. The results indicated that UM06/6-31++G(d,p)/SDD level of theory successfully calculated bond lengths and angles with MADs = 0.02 Å and 0.78 deg., respectively. Also, its combination with CPCM-Pauling-UHF/6-31+G(d)/SDD level of theory in liquid phase effectively computed the redox potential with 0.06 V deviation from the experimental value. Moreover, transferability of the proposed method was studied through ferrocene molecule and its new synthesized derivative.

Oxygen Activation and Radical Transformations in Heme Proteins and Metalloporphyrins

As a result of the adaptation of life to an aerobic environment, nature has evolved a panoply of metalloproteins for oxidative metabolism and protection against reactive oxygen species. Despite the diverse structures and functions of these proteins, they share common mechanistic grounds. An open-shell transition metal like iron or copper is employed to interact with O2 and its derived intermediates such as hydrogen peroxide to afford a variety of metal-oxygen intermediates. These reactive intermediates, including metal-superoxo, -(hydro)peroxo, and high-valent metal-oxo species, are the basis for the various biological functions of O2-utilizing metalloproteins. Collectively, these processes are called oxygen activation. Much of our understanding of the reactivity of these reactive intermediates has come from the study of heme-containing proteins and related metalloporphyrin compounds. These studies not only have deepened our understanding of various functions of heme proteins, such as O2 storage and transport, degradation of reactive oxygen species, redox signaling, and biological oxygenation, etc., but also have driven the development of bioinorganic chemistry and biomimetic catalysis. In this review, we survey the range of O2 activation processes mediated by heme proteins and model compounds with a focus on recent progress in the characterization and reactivity of important iron-oxygen intermediates. Representative reactions initiated by these reactive intermediates as well as some context from prior decades will also be presented. We will discuss the fundamental mechanistic features of these transformations and delineate the underlying structural and electronic factors that contribute to the spectrum of reactivities that has been observed in nature as well as those that have been invented using these paradigms. Given the recent developments in biocatalysis for non-natural chemistries and the renaissance of radical chemistry in organic synthesis, we envision that new enzymatic and synthetic transformations will emerge based on the radical processes mediated by metalloproteins and their synthetic analogs.

Toward Quantitatively Accurate Calculation of the Redox-Associated Acid-Base and Ligand Binding Equilibria of Aquacobalamin

Redox processes in complex transition metal-containing species are often intimately associated with changes in ligand protonation states and metal coordination number. A major challenge is therefore to develop consistent computational approaches for computing pH-dependent redox and ligand dissociation properties of organometallic species. Reduction of the Co center in the vitamin B12 derivative aquacobalamin can be accompanied by ligand dissociation, protonation, or both, making these properties difficult to compute accurately. We examine this challenge here by using density functional theory and continuum solvation to compute Co-ligand binding equilibrium constants (Kon/off), pKas, and reduction potentials for models of aquacobalamin in aqueous solution. We consider two models for cobalamin ligand coordination: the first follows the hexa, penta, tetra coordination scheme for Co(III), Co(II), and Co(I) species, respectively, and the second model features saturation of each vacant axial coordination site on Co(II) and Co(I) species with a single, explicit water molecule to maintain six directly interacting ligands or water molecules in each oxidation state. Comparing these two coordination schemes in combination with five dispersion-corrected density functionals, we find that the accuracy of the computed properties is largely independent of the scheme used, but including only a continuum representation of the solvent yields marginally better results than saturating the first solvation shell around Co throughout. PBE performs best, displaying balanced accuracy and superior performance overall, with RMS errors of 80 mV for seven reduction potentials, 2.0 log units for five pKas and 2.3 log units for two log Kon/off values for the aquacobalamin system. Furthermore, we find that the BP86 functional commonly used in corrinoid studies suffers from erratic behavior and inaccurate descriptions of Co-axial ligand binding, leading to substantial errors in predicted pKas and Kon/off values. These findings demonstrate the effectiveness of the present approach for computing electrochemical and thermodynamic properties of a complex transition metal-containing cofactor.

Ionization Energies and Aqueous Redox Potentials of Organic Molecules: Comparison of DFT, Correlated ab Initio Theory and Pair Natural Orbital Approaches

The calculation of redox potentials involves large energetic terms arising from gas phase ionization energies, thermodynamic contributions, and solvation energies of the reduced and oxidized species. In this work we study the performance of a wide range of wave function and density functional theory methods for the prediction of ionization energies and aqueous one-electron oxidation potentials of a set of 19 organic molecules. Emphasis is placed on evaluating methods that employ the computationally efficient local pair natural orbital (LPNO) approach, as well as several implementations of coupled cluster theory and explicitly correlated F12 methods. The electronic energies are combined with implicit solvation models for the solvation energies. With the exception of MP2 and its variants, which suffer from enormous errors arising at least partially from the poor Hartree-Fock reference, ionization energies can be systematically predicted with average errors below 0.1 eV for most of the correlated wave function based methods studies here, provided basis set extrapolation is performed. LPNO methods are the most efficient way to achieve this type of accuracy. DFT methods show in general larger errors and suffer from inconsistent behavior. The only exception is the M06-2X functional which is found to be competitive with the best LPNO-based approaches for ionization energies. Importantly, the limiting factor for the calculation of accurate redox potentials is the solvation energy. The errors in the predicted solvation energies by all continuum solvation models tested in this work dominate the final computed reduction potential, resulting in average errors typically in excess of 0.3 V and hence obscuring the gains that arise from choosing a more accurate electronic structure method.

Transition metal complexes and the activation of dioxygen

In order to address how diverse metalloprotein active sites, in particular those containing iron and copper, guide O₂binding and activation processes to perform diverse functions, studies of synthetic models of the active sites have been performed. These studies have led to deep, fundamental chemical insights into how O₂coordinates to mono- and multinuclear Fe and Cu centers and is reduced to superoxo, peroxo, hydroperoxo, and, after O-O bond scission, oxo species relevant to proposed intermediates in catalysis. Recent advances in understanding the various factors that influence the course of O₂activation by Fe and Cu complexes are surveyed, with an emphasis on evaluating the structure, bonding, and reactivity of intermediates involved. The discussion is guided by an overarching mechanistic paradigm, with differences in detail due to the involvement of disparate metal ions, nuclearities, geometries, and supporting ligands providing a rich tapestry of reaction pathways by which O₂is activated at Fe and Cu sites.

Computational electrochemistry: prediction of liquid-phase reduction potentials

This article reviews recent developments and applications in the area of computational electrochemistry. Our focus is on predicting the reduction potentials of electron transfer and other electrochemical reactions and half-reactions in both aqueous and nonaqueous solutions. Topics covered include various computational protocols that combine quantum mechanical electronic structure methods (such as density functional theory) with implicit-solvent models, explicit-solvent protocols that employ Monte Carlo or molecular dynamics simulations (for example, Car-Parrinello molecular dynamics using the grand canonical ensemble formalism), and the Marcus theory of electronic charge transfer. We also review computational approaches based on empirical relationships between molecular and electronic structure and electron transfer reactivity. The scope of the implicit-solvent protocols is emphasized, and the present status of the theory and future directions are outlined.

Applications of density functional theory to iron-containing molecules of bioinorganic interest

The past decades have seen an explosive growth in the application of density functional theory (DFT) methods to molecular systems that are of interest in a variety of scientific fields. Owing to its balanced accuracy and efficiency, DFT plays particularly useful roles in the theoretical investigation of large molecules. Even for biological molecules such as proteins, DFT finds application in the form of, e.g., hybrid quantum mechanics and molecular mechanics (QM/MM), in which DFT may be used as a QM method to describe a higher prioritized region in the system, while a MM force field may be used to describe remaining atoms. Iron-containing molecules are particularly important targets of DFT calculations. From the viewpoint of chemistry, this is mainly because iron is abundant on earth, iron plays powerful (and often enigmatic) roles in enzyme catalysis, and iron thus has the great potential for biomimetic catalysis of chemically difficult transformations. In this paper, we present a brief overview of several recent applications of DFT to iron-containing non-heme synthetic complexes, heme-type cytochrome P450 enzymes, and non-heme iron enzymes, all of which are of particular interest in the field of bioinorganic chemistry. Emphasis will be placed on our own work.

A DFT study of a new class of gold nanocluster-photochrome multi-functional switches

The structural and electronic properties of dithienylethene photochromic molecules grafted onto a Au₂₅ nanocluster are reviewed and electron/energy transfers are discussed with the help of (TD-)DFT calculations.