Theoretical studies on the redox potentials of Fe dinuclear complexes as models for hydrogenase

Electrocatalytic Reduction of CO2 and H2O with Zn(II) Complexes Through Metal-Ligand Cooperation

Air and water-stable zinc (II) complexes of neutral pincer bis(diphenylphosphino)-2,6-di(amino)pyridine ("PN3P") ligands are reported. These compounds, [Zn(κ2-2,6-{Ph2PNR}2(NC5H3))Br2] (R=Me, 1; R=H, 2), were shown to be capable of electrocatalytic reduction of CO2 at -2.3 V vs. Fc+/0 to selectively yield CO in mixed water/acetonitrile solutions. These complexes also electrocatalytically generate H2 from water in acetonitrile solutions, at the same potential, with Faradaic efficiencies of up to 90 %. DFT computations support a proposed mechanism involving the first reduction of 1 or 2 occurring at the PN3P ligand. Furthermore, computational analysis suggested a mechanism involving metal-ligand cooperation of a Lewis acidic Zn(II) and a basic ligand.

Complexation and redox chemistry of neptunium, plutonium and americium with a hydroxylaminato ligand

There is significant interest in ligands that can stabilize actinide ions in oxidation states that can be exploited to chemically differentiate 5f and 4f elements. Applications range from developing large-scale actinide separation strategies for nuclear industry processing to carrying out analytical studies that support environmental monitoring and remediation efforts. Here, we report syntheses and characterization of Np(iv), Pu(iv) and Am(iii) complexes with N-tert-butyl-N-(pyridin-2-yl)hydroxylaminato, [2-(^tBuNO)py]⁻(interchangeable hereafter with [(^tBuNO)py]⁻), a ligand which was previously found to impart remarkable stability to cerium in the +4 oxidation state. An[(^tBuNO)py]₄ (An = Pu, 1; Np, 2) have been synthesized, characterized by X-ray diffraction, X-ray absorption, ¹H NMR and UV-vis-NIR spectroscopies, and cyclic voltammetry, along with computational modeling and analysis. In the case of Pu, oxidation of Pu(iii) to Pu(iv) was observed upon complexation with the [(^tBuNO)py]⁻ ligand. The Pu complex 1 and Np complex 2 were also isolated directly from Pu(iv) and Np(iv) precursors. Electrochemical measurements indicate that a Pu(iii) species can be accessed upon one-electron reduction of 1 with a large negative reduction potential (E_1/2 = −2.26 V vs. Fc^+/0). Applying oxidation potentials to 1 and 2 resulted in ligand-centered electron transfer reactions, which is different from the previously reported redox chemistry of U^IV[(^tBuNO)py]₄ that revealed a stable U(v) product. Treatment of an anhydrous Am(iii) precursor with the [(^tBuNO)py]⁻ ligand did not result in oxidation to Am(iv). Instead, the dimeric complex [Am^III(μ₂-(^tBuNO)py)((^tBuNO)py)₂]₂ (3) was isolated. Complex 3 is a rare example of a structurally characterized non-aqueous Am-containing molecular complex prepared using inert atmosphere techniques. Predicted redox potentials from density functional theory calculations show a trivalent accessibility trend of U(iii) < Np(iii) < Pu(iii) and that the higher oxidation states of actinides (i.e., +5 for Np and Pu and +4 for Am) are not stabilized by [2-(^tBuNO)py]⁻, in good agreement with experimental observations.

The coordination modes and electronic properties of a strongly coordinating hydroxylaminato ligand with Np, Pu and Am were investigated.Complexes were characterized by a range of experimental and computational techniques.

Electrocatalytic Behavior of Tetrathiafulvalene (TTF) and Extended Tetrathiafulvalene (exTTF) [FeFe] Hydrogenase Mimics

TTF- and exTTF-containing [(μ-S₂)Fe₂(CO)₆] complexes have been prepared by the photochemical reaction of TTF or exTTF and [(μ-S₂)Fe₂(CO)₆]. These complexes are able to interact with PAHs. In the absence of air and in acid media an electrocatalytic dihydrogen evolution reaction (HER) occurs, similarly to analogous [(μ-S₂)Fe₂(CO)₆] complexes. However, in the presence of air, the TTF and exTTF organic moieties strongly influence the electrochemistry of these systems. The reported data may be valuable in the design of [FeFe] hydrogenase mimics able to combine the HER properties of the [FeFe] cores with the unique TTF properties.

A Comparative Review of Metal-Based Charge Carriers in Nonaqueous Flow Batteries

Energy storage is becoming the chief barrier to the utilization of more renewable energy sources on the grid. With independent service operators aiming to acquire gigawatts in the next 10-20 years, there is a large need to develop a suite of new storage technologies. Redox flow batteries (RFB) may be part of the solution if certain key barriers are overcome. This Review focuses on a particular kind of RFB based on nonaqueous media that promises to meet the challenge through higher voltages than the organic and aqueous variants. This class of RFB is divided into three groups: molecular, macromolecular, and redox-targeted systems. The growing field of theoretical modeling is also reviewed and discussed.

Electronic and geometric structure effects on one-electron oxidation of first-row transition metals in the same ligand framework

Developing new transition metal catalysts requires understanding of how both metal and ligand properties determine reactivity. Since metal complexes bearing ligands of the Py5 family (2,6-bis-[(2-pyridyl)methyl]pyridine) have been employed in many fields in the past 20 years, we set out here to understand their redox properties by studying a series of base metal ions (M = Mn, Fe, Co, and Ni) within the Py5OH (pyridine-2,6-diylbis[di-(pyridin-2-yl)methanol]) variant. Both reduced (M^II) and the one-electron oxidized (M^III) species were carefully characterized using a combination of X-ray crystallography, X-ray absorption spectroscopy, cyclic voltammetry, and density-functional theory calculations. The observed metal-ligand interactions and electrochemical properties do not always follow consistent trends along the periodic table. We demonstrate that this observation cannot be explained by only considering orbital and geometric relaxation, and that spin multiplicity changes needed to be included into the DFT calculations to reproduce and understand these trends. In addition, exchange reactions of the sixth ligand coordinated to the metal, were analysed. Finally, by including published data of the extensively characterised Py5OMe (pyridine-2,6-diylbis[di-(pyridin-2-yl)methoxymethane])complexes, the special characteristics of the less common Py5OH ligand were extracted. This comparison highlights the non-innocent effect of the distal OH functionalization on the geometry, and consequently on the electronic structure of the metal complexes. Together, this gives a complete analysis of metal and ligand degrees of freedom for these base metal complexes, while also providing general insights into how to control electrochemical processes of transition metal complexes.

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.

A Density Functional Theory-Based Scheme to Compute the Redox Potential of a Transition Metal Complex: Applications to Heme Compound

We estimated the redox potential of a model heme compound by using the combination of our density functionals with a computational scheme, which corrects the solvation energy to the normal solvent model. Among many density functionals, the LC-BOP12 functional gave the smallest mean absolute error of 0.16 V in the test molecular sets. The application of these methods revealed that the redox potential of a model heme can be controlled within 200 mV by changing the protonation state and even within 20 mV by the flipping of the ligand histidine. In addition, the redox potential depends on the inverse of the dielectric constant, which controls the surroundings. The computational results also imply that a system with a low dielectric constant avoids the charged molecule by controlling either the redox potential or the protonation system.

Ultrafast Photodynamics of Cyano-Functionalized [FeFe] Hydrogenase Model Compounds

[FeFe] hydrogenases are efficient enzymes that produce hydrogen gas under mild conditions. Synthetic model compounds containing all CO or mixed CO/PMe3 ligands were previously studied by us and others with ultrafast ultraviolet or visible pump-infrared probe spectroscopy in an effort to better understand the function and interactions of the active site with light. Studies of anionic species containing cyano groups, which more closely match the biological active site, have been elusive. In this work, two model compounds dissolved in room-temperature acetonitrile solution were examined: [Fe2(μ-S2C3H6)(CO)4(CN)2]2- (1) and [Fe2(μ-S2C2H4)(CO)4(CN)2]2- (2). These species exhibit long-lived transient signals consistent with loss of one CO ligand with potential isomerization of newly formed ground electronic state photoproducts, as previously observed with all-CO and CO/PMe3-containing models. We find no evidence for fast (ca. 150 ps) relaxation seen in the all-CO and CO/PMe3 compounds because of the absence of the metal-to-metal charge transfer band in the cyano-functionalized models. These results indicate that incorporation of cyano ligands may significantly alter the electronic properties and photoproducts produced immediately after photoexcitation, which may influence the catalytic activity of model compounds when attached to photosensitizers.

Using theoretical calculations to predict the redox potential of mononuclear manganese complexes

Prediction of redox potential allows chemists to rationally design metal complexes with a desired redox activity.

Mechanistic Insight into Electrocatalytic H2 Production by [Fe2(CN){μ-CN(Me)2}(μ-CO)(CO)(Cp)2]: Effects of Dithiolate Replacement in [FeFe] Hydrogenase Models

DFT has been used to investigate viable mechanisms of the hydrogen evolution reaction (HER) electrocatalyzed by [Fe₂(CN){μ-CN(Me)₂}(μ-CO)(CO)(Cp)₂] (1) in AcOH. Molecular details underlying the proposed ECEC electrochemical sequence have been studied, and the key functionalities of CN^- and amino-carbyne ligands have been elucidated. After the first reduction, CN^- works as a relay for the first proton from AcOH to the carbyne, with this ligand serving as the main electron acceptor for both reduction steps. After the second reduction, a second protonation occurs at CN^- that forms a Fe(CNH) moiety: i.e., the acidic source for the H₂ generation. The hydride (formally 2e/H⁺), necessary to the heterocoupling with H⁺ is thus provided by the μ-CN(Me)₂ ligand and not by Fe centers, as occurs in typical L₆Fe₂S₂ derivatives modeling the hydrogenase active site. It is remarkable, in this regard, that CN^- plays a role more subtle than that previously expected (increasing electron density at Fe atoms). In addition, the role of AcOH in shuttling protons from CN^- to CN(Me)₂ is highlighted. The incompetence for the HER of the related species [Fe₂{μ-CN(Me)₂}(μ-CO)(CO)₂(Cp)₂]⁺ (2⁺) has been investigated and attributed to the loss of proton responsiveness caused by CN^- replacement with CO. In the context of hydrogenase mimicry, an implication of this study is that the dithiolate strap, normally present in all synthetic models, can be removed from the Fe₂ core without loss of HER, but the redox and acid-base processes underlying turnover switch from a metal-based to a ligand-based chemistry. The versatile nature of the carbyne, once incorporated in the Fe₂ scaffold, could be exploited to develop more active and robust catalysts for the HER.

Molecular Cobalt Complexes with Pendant Amines for Selective Electrocatalytic Reduction of Carbon Dioxide to Formic Acid

We report here on a new series of CO₂-reducing molecular catalysts based on Earth-abundant elements that are very selective for the production of formic acid in dimethylformamide (DMF)/water mixtures (Faradaic efficiency of 90 ± 10%) at moderate overpotentials (500-700 mV in DMF measured at the middle of the catalytic wave). The [CpCo(P^R₂N^R'₂)I]⁺ compounds contain diphosphine ligands, P^R₂N^R'₂, with two pendant amine residues that act as proton relays during CO₂-reduction catalysis and tune their activity. Four different P^R₂N^R'₂ ligands with cyclohexyl or phenyl substituents on phosphorus and benzyl or phenyl substituents on nitrogen were employed, and the compound with the most electron-donating phosphine ligand and the most basic amine functions performs best among the series, with turnover frequency >1000 s^-1. State-of-the-art benchmarking of catalytic performances ranks this new class of cobalt-based complexes among the most promising CO₂-to-formic acid reducing catalysts developed to date; addressing the stability issues would allow further improvement. Mechanistic studies and density functional theory simulations confirmed the role of amine groups for stabilizing key intermediates through hydrogen bonding with water molecules during hydride transfer from the Co center to the CO₂ molecule.

Hydrogenase Enzymes and Their Synthetic Models: The Role of Metal Hydrides

Hydrogenase enzymes efficiently process H2 and protons at organometallic FeFe, NiFe, or Fe active sites. Synthetic modeling of the many H2ase states has provided insight into H2ase structure and mechanism, as well as afforded catalysts for the H2 energy vector. Particularly important are hydride-bearing states, with synthetic hydride analogues now known for each hydrogenase class. These hydrides are typically prepared by protonation of low-valent cores. Examples of FeFe and NiFe hydrides derived from H2 have also been prepared. Such chemistry is more developed than mimicry of the redox-inactive monoFe enzyme, although functional models of the latter are now emerging. Advances in physical and theoretical characterization of H2ase enzymes and synthetic models have proven key to the study of hydrides in particular, and will guide modeling efforts toward more robust and active species optimized for practical applications.

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.

Counterion influence on the N-I-N halogen bond

A detailed investigation of the influence of counterions on the [N-I-N]+ halogen bond in solution, in the solid state and in silico is presented. Translational diffusion coefficients indicate close attachment of counterions to the cationic, three-center halogen bond in dichloromethane solution. Isotopic perturbation of equilibrium NMR studies performed on isotopologue mixtures of regioselectively deuterated and nondeuterated analogues of the model system showed that the counterion is incapable of altering the symmetry of the [N-I-N]+ halogen bond. This symmetry remains even in the presence of an unfavorable geometric restraint. A high preference for the symmetric geometry was found also in the solid state by single crystal X-ray crystallography. Molecular systems encompassing weakly coordinating counterions behave similarly to the corresponding silver(i) centered coordination complexes. In contrast, systems possessing moderately or strongly coordinating anions show a distinctly different behavior. Such silver(i) complexes are converted into multi-coordinate geometries with strong Ag-O bonds, whereas the iodine centered systems remain linear and lack direct charge transfer interaction with the counterion, as verified by 15N NMR and DFT computation. This suggests that the [N-I-N]+ halogen bond may not be satisfactorily described in terms of a pure coordination bond typical of transition metal complexes, but as a secondary bond with a substantial charge-transfer character.

Necroptosis-inducing rhenium(V) oxo complexes

Rhenium(V) oxo complexes of general formula [ReO(OMe)(N^N)Cl2], where N^N = 4,7-diphenyl-1,10-phenanthroline, 1, or 3,4,7,8-tetramethyl-1,10-phenanthroline, 2, effectively kill cancer cells by triggering necroptosis, a non-apoptotic form of cell death. Both complexes evoke necrosome (RIP1-RIP3)-dependent intracellular reactive oxygen species (ROS) production and propidium iodide uptake. The complexes also induce mitochondrial membrane potential depletion, a possible downstream effect of ROS production. Apparently, 1 and 2 are the first rhenium complexes to evoke cellular events consistent with programmed necrosis in cancer cells. Furthermore, 1 and 2 display low acute toxicity in C57BL/6 mice and reasonable stability in fresh human blood.

Accurate standard hydrogen electrode potential and applications to the redox potentials of vitamin C and NAD/NADH

We computationally evaluated the standard hydrogen electrode (SHE) potential in aqueous phase and the Gibbs energy of a proton from the experimental pKa values of alcohol molecules. From the "golden standard" CCSD(T)/aug-cc-pVTZ level calculation, we estimated the SHE potential as 4.48 V, which is very close to the IUPAC-recommended experimental value of 4.44 V. As applications to the Gaussian-3 (G3) methods, which also reproduce the "golden standard" level calculations, we computed various pKa values and redox potentials for a vitamin series. For vitamin C, we support the experimental result of +0.35 V and predict the pKa value of d-ascorbic acid to be 3.7-3.9. Using a model molecule for nicotinamide adenine dinucleotide (NAD), we reproduced the redox potential and determined the order of the proton/electron addition, based on both the proton affinity and redox potential.

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.