Density Functional Calculations for Prediction of 57Fe Mössbauer Isomer Shifts and Quadrupole Splittings in β-Diketiminate Complexes

A structural and functional model for alkene dioxygenases

In this article, we report sterically-controlled iron sites based on non-chelating bulky imidazole ligands. Adding 6 equiv. of 1,2-dimethylimidazole (1,2-Me2Im) to Fe(OTf)2⋅2CH3CN affords the first example of a 5-coordinate imidazole‑iron complex ([Fe(1,2-Me2Im)5](OTf)2, 1). The structure is distorted square pyramidal (τ5 = 0.41). When an iPr group is substituted for the methyl group at the 2-position on the imidazole (2-iPr-1-MeIm), the 14-electron complex ([Fe(2-iPr-1-MeIm)4](OTf)2, 2) is obtained. This complex exhibits slightly distorted tetrahedral geometry (τ'4 = 0.93) with four N-donors and serves as a 4-His iron structural model complex for carotenoid cleavage dioxygenases (CCD). The electronic structure of 1 and 2 were characterized by Mössbauer spectroscopy. Reactions of 1 and 2 with model olefin substrates (1-R-4-(1-methoxyprop-1-en-2-yl)benzene; R = Me or Br) in the presence of oxygen result in olefin cleavage yielding ketone and aldehyde products, although 2 yields more products than 1. Support for a proposed reaction mechanism for 2 is offered from Density Functional Theory (DFT) calculations.

Extensive reference set and refined computational protocol for calculations of 57Fe Mössbauer parameters

Mössbauer spectroscopy is a powerful technique for probing the local electronic structure of iron compounds, because it reports in an element-selective manner on both the oxidation state and coordination environment of the Fe ion. Computational prediction of the two main Mössbauer parameters, isomer shift (δ) and quadrupole splitting (ΔEQ), has long been targeted by quantum chemical studies, and useful protocols based on density functional theory have been proposed. Here we present an extensive curated reference set of Fe compounds that is considerably larger and more diverse than literature precedents. We make a distinction between low-temperature and high-temperature experimental subgroups. This set is employed for optimizing a refined computational protocol utilizing the scalar version of the exact 2-component (X2C) Hamiltonian with the finite nucleus approximation. Attention is devoted to having an accurate and flexible all-electron basis set for Fe. We assess the performance of several DFT methods that cover all representative families and rungs of functionals and find that hybrid functionals with ca. 25-30% exact exchange offer the best accuracy for isomer shifts. The work establishes a refined general protocol of wide applicability that achieves good performance for the prediction of isomer shifts in a wider variety of systems than before, but the limitations of DFT for quadrupole splittings are also highlighted. Finally, comparison of calculated values with high-temperature experimental results shows that the use of an empirical correction factor is required to account for the second-order Doppler shift and to achieve the same quality of correlation as with the low-temperature data.

A Caged Neutral 17-Valence-Electron Iron(I) Radical [Fe(CO)2(Cl)(P((CH2)10)3P)]•: Synthetic, Structural, Spectroscopic, Redox, and Computational Studies

UV irradiation of yellow CH2Cl2 solutions of trans-Fe(CO)3(P((CH2)10)3P) (2a) and PMe3 (10 equiv) gives, in addition to the previously reported dibridgehead diphosphine P((CH2)10)3P (46%), a green paramagnetic complex that crystallography shows to be the trigonal-bipyramidal iron(I) radical trans-[Fe(CO)2(Cl)(P((CH2)10)3P)]• (1a•; 31% after workup). This is a rare example of an isolable species of the formula [Fe(CO)4-n(L)n(X)]• (n = 0-3, L = two-electron-donor ligand; X = one-electron-donor ligand). Analogous precursors with longer P(CH2)nP segments (n = 12, 14, 16, 18) give only the demetalated diphosphines, and a rationale is proposed. The magnetic susceptibility of 1a•, assayed by Evans' method and SQUID measurements, indicates a spin (S) of 1/2. Cyclic voltammetry shows that 1a• undergoes a partially reversible one-electron oxidation, but no facile reduction. The UV-visible, EPR, and 57Fe Mössbauer spectra are analyzed in detail. Complex 2a is similarly studied, and, despite the extra valence electron, exhibits a comparable oxidation potential (ΔE1/2 ≤ 0.04 V). The crystal structure shows a cage conformation, solvation level, disorder motif, and unit cell parameters essentially identical to those of 1a•. DFT calculations provide much insight regarding the structural, redox, and spectroscopic properties.

Calculation of electric field gradients with the exact two-component (X2C) quasi-relativistic method and its local approximations

When calculating electric field gradients (EFGs), relativistic and electron correlation effects are crucial for obtaining accurate results, and the commonly used density functional methods produce unsatisfactory results, especially for heavy elements and/or strongly correlated systems. In this work, a stand-alone program is presented, which enables calculation of EFGs from the molecular orbitals supplied by an external high accuracy quantum chemical calculation and includes relativistic effects through the exact two-component (X2C) formalism and efficient local approximations to it. Application to BiN and BiP molecules shows that a high precision can be achieved in the calculation of nuclear quadrupole coupling constants of 209Bi by combining advanced ab initio methods with the X2C approach. For seventeen iron compounds, the Mössbauer nuclear quadrupole splittings (NQS) of 57Fe calculated using a double-hybrid functional method are in very good agreement with the experimental values. It is shown that, for strongly correlated molecules, the double-hybrid functionals are much more accurate than the commonly used hybrid functionals. The computer program developed in this study furnishes a useful utility for obtaining EFGs and related nuclear properties with high accuracy.

Axial Ligation Impedes Proton-Coupled Electron-Transfer Reactivity of a Synthetic Compound-I Analogue

The nature of the axial ligand in high-valent iron-oxo heme enzyme intermediates and related synthetic catalysts is a critical structural element for controlling proton-coupled electron-transfer (PCET) reactivity of these species. Herein, we describe the generation and characterization of three new 6-coordinate, iron(IV)-oxo porphyrinoid-π-cation-radical complexes and report their PCET reactivity together with a previously published 5-coordinate analogue, FeIV(O)(TBP8Cz+•) (TBP8Cz = octakis(p-tert-butylphenyl)corrolazinato3-) (2) (Cho, K. A high-valent iron-oxo corrolazine activates C-H bonds via hydrogen-atom transfer. J. Am. Chem. Soc. 2012, 134, 7392-7399). The new complexes FeIV(O)(TBP8Cz+•)(L) (L = 1-methyl imidazole (1-MeIm) (4a), 4-dimethylaminopyridine (DMAP) (4b), cyanide (CN-)(4c)) can be generated from either oxidation of the ferric precursors or by addition of L to the Compound-I (Cpd-I) analogue at low temperatures. These complexes were characterized by UV-vis, electron paramagnetic resonance (EPR), and Mössbauer spectroscopies, and cryospray ionization mass spectrometry (CSI-MS). Kinetic studies using 4-OMe-TEMPOH as a test substrate indicate that coordination of a sixth axial ligand dramatically lowers the PCET reactivity of the Cpd-I analogue (rates up to 7000 times slower). Extensive density functional theory (DFT) calculations together with the experimental data show that the trend in reactivity with the axial ligands does not correlate with the thermodynamic driving force for these reactions or the calculated strengths of the O-H bonds being formed in the FeIV(O-H) products, pointing to non-Bell-Evans-Polanyi behavior. However, the PCET reactivity does follow a trend with the bracketed reduction potential of Cpd-I analogues and calculated electron affinities. The combined data suggest a concerted mechanism (a concerted proton electron transfer (CPET)) and an asynchronous movement of the electron/proton pair in the transition state.

Evaluating Diazene to N2 Interconversion at Iron-Sulfur Complexes

Biological N2 reduction occurs at sulfur-rich multiiron sites, and an interesting potential pathway is concerted double reduction/ protonation of bridging N2 through PCET. Here, we test the feasibility of using synthetic sulfur-supported diiron complexes to mimic this pathway. Oxidative proton transfer from μ-η1 : η1-diazene (HN=NH) is the microscopic reverse of the proposed N2 fixation pathway, revealing the energetics of the process. Previously, Sellmann assigned the purple metastable product from two-electron oxidation of [{Fe2+(PPr3)L1}2(μ-η1 : η1-N2H2)] (L1=tetradentate SSSS ligand) at -78 °C as [{Fe2+(PPr3)L1}2(μ-η1 : η1-N2)]2+, which would come from double PCET from diazene to sulfur atoms of the supporting ligands. Using resonance Raman, Mössbauer, NMR, and EPR spectroscopies in conjunction with DFT calculations, we show that the product is not an N2 complex. Instead, the data are most consistent with the spectroscopically observed species being the mononuclear iron(III) diazene complex [{Fe(PPr3)L1}(η2-N2H2)]+. Calculations indicate that the proposed double PCET has a barrier that is too high for proton transfer at the reaction temperature. Also, PCET from the bridging diazene is highly exergonic as a result of the high Fe3+/2+ redox potential, indicating that the reverse N2 protonation would be too endergonic to proceed. This system establishes the "ground rules" for designing reversible N2/N2H2 interconversion through PCET, such as tuning the redox potentials of the metal sites.

Iron(iv) alkyl complexes: electronic structure contributions to Fe-C bond homolysis and migration reactions that form N-C bonds from N2

High-valent iron alkyl complexes are rare, as they are prone to Fe-C bond homolysis. Here, we describe an unusual way to access formally iron(iv) alkyl complexes through double silylation of iron(i) alkyl dinitrogen complexes to form an NNSi2 group. Spectroscopically validated computations show that the disilylehydrazido(2-) ligand stabilizes the formal iron(iv) oxidation state through a strongly covalent Fe-N π-interaction, in which one π-bond fits an "inverted field" description. This means that the two bonding electrons are localized more on the metal than the ligand, and thus an iron(ii) resonance structure is a significant contributor, similar to the previously-reported phenyl analogue. However, in contrast to the phenyl complex which has an S = 1 ground state, the ground state of the alkyl complex is S = 2, which places one electron in the π* orbital, leading to longer and weaker Fe-N bonds. The reactivity of these hydrazido(2-) complexes is dependent on the steric and electronic properties of the specific alkyl group. When the alkyl group is the bulky trimethylsilylmethyl, the formally iron(iv) species is stable at room temperature and no migration of the alkyl ligand is observed. However, the analogous complex with the smaller methyl ligand does indeed undergo migration of the carbon-based ligand to the NNSi2 group to form a new N-C bond. This migration is followed by isomerization of the hydrazido ligand, and the product exists as two isomers that have distinct η1 and η2 binding of the hydrazido group. Lastly, when the alkyl group is benzyl, the Fe-C bond homolyzes to give a three-coordinate hydrazido(2-) complex which is likely due to the greater stability of a benzyl radical compared to that for methyl or trimethylsilylmethyl. These studies demonstrate the availability of a hydrocarbyl migration pathway at formally iron(iv) centers to form new N-C bonds directly to N2, though product selectivity is highly dependent on the identity of the migrating group.

Spin effect on redox acceleration and regioselectivity in Fe-catalyzed alkyne hydrosilylation

Iron catalysts are ideal transition metal catalysts because of the Earths abundant, cheap, biocompatible features of iron salts. Iron catalysts often have unique open-shell structures that easily undergo spin crossover in chemical transformations, a feature rarely found in noble metal catalysts. Unfortunately, little is known currently about how the open-shell structure and spin crossover affect the reactivity and selectivity of iron catalysts, which makes the development of iron catalysts a low efficient trial-and-error program. In this paper, a combination of experiments and theoretical calculations revealed that the iron-catalyzed hydrosilylation of alkynes is typical spin-crossover catalysis. Deep insight into the electronic structures of a set of well-defined open-shell active formal Fe(0) catalysts revealed that the spin-delocalization between the iron center and the 1,10-phenanthroline ligand effectively regulates the iron center's spin and oxidation state to meet the opposite electrostatic requirements of oxidative addition and reductive elimination, respectively, and the spin crossover is essential for this electron transfer process. The triplet transition state was essential for achieving high regioselectivity through tuning the nonbonding interactions. These findings provide an important reference for understanding the effect of catalyst spin state on reaction. It is inspiring for the development of iron catalysts and other Earth-abundant metal catalysts, especially from the point of view of ligand development.

Valence Delocalization and Metal-Metal Bonding in Carbon-Bridged Mixed-Valence Iron Complexes

The carbide ligand in the iron-molybdenum cofactor (FeMoco) in nitrogenase bridges iron atoms in different oxidation states, yet it is difficult to discern its ability to mediate magnetic exchange interactions due to the structural complexity of the cofactor. Here, we describe two mixed-valent diiron complexes with C-based ketenylidene bridging ligands, and compare the carbon bridges with the more familiar sulfur bridges. The ground state of the [Fe2 (μ-CCO)2 ]+ complex with two carbon bridges (4) is S=1 / 2 ${{ 1/2 }}$ , and it is valence delocalized on the Mössbauer timescale with a small thermal barrier for electron hopping that stems from the low Fe-C force constant. In contrast, one-electron reduction of the [Fe2 (μ-CCO)] complex with one carbon bridge (2) affords a mixed-valence species with a high-spin ground state (S=7 / 2 ${ 7/2 }$ ), and the Fe-Fe distance contracts by 1 Å. Spectroscopic, magnetic, and computational studies of the latter reveal an Fe-Fe bonding interaction that leads to complete valence delocalization. Analysis of near-IR intervalence charge transfer transitions in 5 indicates a very large double exchange constant (B) in the range of 780-965 cm-1 . These results show that carbon bridges are extremely effective at stabilizing valence delocalized ground states in mixed-valent iron dimers.

Mössbauer and Nuclear Resonance Vibrational Spectroscopy Studies of Iron Species Involved in N-N Bond Cleavage

Diketiminate-supported iron complexes are capable of cleaving the strong triple bond of N2 to give a tetra-iron complex with two nitrides (Rodriguez et al., Science, 2011, 334, 780-783). The mechanism of this reaction has been difficult to determine, but a transient green species was observed during the reaction that corresponds to a potential intermediate. Here, we describe studies aiming to identify the characteristics of this intermediate, using a range of spectroscopic techniques, including Mössbauer spectroscopy, electronic absorption spectroscopy, Raman spectroscopy, nuclear magnetic resonance (NMR) spectroscopy, and nuclear resonance vibrational spectroscopy (NRVS) complemented by density functional theory (DFT) calculations. We successfully elucidated the nature of the starting iron(II) species and the bis(nitride) species in THF solution, and in each case, THF breaks up the multiiron species. Various observations on the green intermediate species indicate that it has one N2 per two Fe atoms, has THF associated with it, and has NRVS features indicative of bridging N2. Computational models with a formally diiron(0)-N2 core are most consistent with the accumulated data, and on this basis, a mechanism for N2 splitting is suggested. This work shows the power of combining NRVS, Mössbauer, NMR, and vibrational spectroscopies with computations for revealing the nature of transient iron species during N2 cleavage.

Well-Defined Iron Sites in Crystalline Carbon Nitride

Carbon nitride materials can be hosts for transition metal sites, but Mössbauer studies on iron complexes in carbon nitrides have always shown a mixture of environments and oxidation states. Here we describe the synthesis and characterization of a crystalline carbon nitride with stoichiometric iron sites that all have the same environment. The material (formula C6N9H2Fe0.4Li1.2Cl, abbreviated PTI/FeCl2) is derived from reacting poly(triazine imide)·LiCl (PTI/LiCl) with a low-melting FeCl2/KCl flux, followed by anaerobic rinsing with methanol. X-ray diffraction, X-ray absorption and Mössbauer spectroscopies, and SQUID magnetometry indicate that there are tetrahedral high-spin iron(II) sites throughout the material, all having the same geometry. The material is active for electrocatalytic nitrate reduction to ammonia, with a production rate of ca. 0.1 mmol cm-2 h-1 and Faradaic efficiency of ca. 80% at -0.80 V vs RHE.

Magnetically Induced Binary Ferrocene with Oxidized Iron

Ferrocene is perhaps the most popular and well-studied organometallic molecule, but our understanding of its structure and electronic properties has not changed for more than 70 years. In particular, all previous attempts of chemically oxidizing pure ferrocene by binding directly to the iron center have been unsuccessful, and no significant change in structure or magnetism has been reported. Using a metal organic framework host material, we were able to fundamentally change the electronic and magnetic structure of ferrocene to take on a never-before observed physically stretched/bent high-spin Fe(II) state, which readily accepts O2 from air, chemically oxidizing the iron from Fe(II) to Fe(III). We also show that the binding of oxygen is reversible through temperature swing experiments. Our analysis is based on combining Mößbauer spectroscopy, extended X-ray absorption fine structure, in situ infrared, SQUID, thermal gravimetric analysis, and energy dispersive X-ray fluorescence spectroscopy measurements with ab initio modeling.

Desulfurization and N2 Binding at an Iron Complex Derived from the C-S Activation of Benzothiophene

Metal insertion into the C-S bonds of thiophenes is a facile route to interesting polydentate ligand scaffolds with C and S donors. Here, we describe iron-mediated C-S activation of a diphenylphosphine-functionalized benzothiophene proligand. Metalation of the proligand with "tetrakis(trimethylphosphine)iron" gives an initial five-coordinate, diamagnetic iron(II) species with two PMe3 ligands and a dianionic PCS pincer ligand. Upon one-electron reduction, a reactive anionic iron(I) complex is formed. This species then undergoes deep-seated changes, notably cleavage of C-S and C-P bonds in the supporting ligand. Substantial coordination sphere alterations accompany the ligand C-S bond activation, including loss of a sulfur anion from the S-Fe-C metallacycle and reorganization of the two PMe3 ligands. The resulting desulfurized six-coordinate PCC iron complex also has an N2 ligand trans to the vinyl C. Reducing this complex then cleaves a C-P bond in the appended diphenylphosphine, giving a phosphido arm. These ligand transformations demonstrate novel approaches to pincers with thiolates and phosphides, which would be difficult to synthesize using typical methods through free ligand salts.

Investigating Metal-Metal Bond Polarization in a Heteroleptic Tris-Ylide Diiron System

This article describes the synthesis, characterization, and S-atom transfer reactivity of a series of C3v-symmetric diiron complexes. The iron centers in each complex are coordinated in distinct ligand environments, with one (FeN) bound in a pseudo-trigonal bipyramidal geometry by three phosphinimine nitrogens in the equatorial plane, a tertiary amine, and the second metal center (FeC). FeC is coordinated, in turn, by FeN, three ylidic carbons in a trigonal plane, and, in certain cases, by an axial oxygen donor. The three alkyl donors at FeC form through the reduction of the appended N═PMe3 arms of the monometallic parent complex. The complexes were studied crystallographically, spectroscopically (NMR, UV-vis, and Mössbauer), and computationally (DFT, CASSCF) and found to be high-spin throughout, with short Fe-Fe distances that belie weak orbital overlap between the two metals. Further, the redox nature of this series allowed for the determination that oxidation is localized to the FeC. S-atom transfer chemistry resulted in the formal insertion of a S atom into the Fe-Fe bond of the reduced diiron complex to form a mixture of Fe4S and Fe4S2 products.

Overlooked Transformation of Nitrated Polycyclic Aromatic Hydrocarbons in Natural Waters: Role of Self-Photosensitization

Photochemical transformation is an important process that involves trace organic contaminants (TrOCs) in sunlit surface waters. However, the environmental implications of their self-photosensitization pathway have been largely overlooked. Here, we selected 1-nitronaphthalene (1NN), a representative nitrated polycyclic aromatic hydrocarbon, to study the self-photosensitization process. We investigated the excited-state properties and relaxation kinetics of 1NN after sunlight absorption. The intrinsic decay rate constants of triplet (³1NN*) and singlet (¹1NN*) excited states were estimated to be 1.5 × 10⁶ and 2.5 × 10⁸ s^-1, respectively. Our results provided quantitative evidence for the environmental relevance of ³1NN* in waters. Possible reactions of ³1NN* with various water components were evaluated. With the reduction and oxidation potentials of -0.37 and 1.95 V, ³1NN* can be either oxidized or reduced by dissolved organic matter isolates and surrogates. We also showed that hydroxyl (^•OH) and sulfate (SO₄^•-) radicals can be generated via the ³1NN*-induced oxidation of inorganic ions (OH^- and SO₄^2-, respectively). We further investigated the reaction kinetics of ³1NN* and OH^- forming ^•OH, an important photoinduced reactive intermediate, through complementary experimental and theoretical approaches. The rate constants for the reactions of ³1NN* with OH^- and 1NN with ^•OH were determined to be 4.22 × 10⁷ and 3.95 ± 0.01 × 10⁹ M^-1 s^-1, respectively. These findings yield new insights into self-photosensitization as a pathway for TrOC attenuation and provide more mechanistic details into their environmental fate.

Two Shared Icosahedral Metallacarboranes through Iron: A Joint Experimental and Theoretical Refinement of Mössbauer Spectrum in [Fe(1,2-C2B9H11)2]Cs

Mössbauer and X-ray photoelectron spectroscopies (XPS) are complemented with high-level quantum-chemical computations in the study of the geometric and electronic structure of the paramagnetic salt of the metallacarborane sandwich complex [Fe(1,2-C₂B₉H₁₁)₂]Cs = FeSanCs. Experimental ⁵⁷Fe isomer shifts and quadrupole splitting parameters are compared with the theoretical prediction, with good agreement. The appearance of two sets of Cs(3d) doublets in the XPS spectrum, separated by 2 eV, indicates that Cs has two different chemical environments due to ease of the Cs⁽⁺⁾ cation moving around the sandwich complex with low-energy barriers, as confirmed by quantum-chemical computations. Several minimum-energy geometries of the FeSanCs structure with the corresponding energies and Mössbauer parameters are discussed, in particular the atomic charges and spin population and the surroundings of the Fe atom in the complex. The Mössbauer spectra were taken at different temperatures showing the presence of a low-spin Fe atom with S = 1/2 and thus confirming a paramagnetic Fe^III species.

Mechanism of Nitrogen-Carbon Bond Formation from Iron(IV) Disilylhydrazido Intermediates during N2 Reduction

We recently reported a reaction sequence that activates C-H bonds in simple arenes as well as the N-N triple bond in N₂, delivering the aryl group to N₂ to form a new N-C bond (Nature 2020, 584, 221). This enables the transformation of abundant feedstocks (arenes and N₂) into N-containing organic compounds. The key N-C bond forming step occurs upon partial silylation of N₂. However, the pathway through which reduction, silylation, and migration occurred was unknown. Here, we describe synthetic, structural, magnetic, spectroscopic, kinetic, and computational studies that elucidate the steps of this transformation. N₂ must be silylated twice at the distal N atom before aryl migration can occur, and sequential silyl radical and silyl cation addition is a kinetically competent pathway to a formally iron(IV)-NN(SiMe₃)₂ intermediate that can be isolated at low temperature. Kinetic studies show its first-order conversion to the migrated product, and DFT calculations indicate a concerted transition state for migration. The electronic structure of the formally iron(IV) intermediate is examined using DFT and CASSCF calculations, which reveal contributions from iron(II) and iron(III) resonance forms with oxidized NNSi₂ ligands. The depletion of electron density from the Fe-coordinated N atom makes it electrophilic enough to accept the incoming aryl group. This new pathway for the N-C bond formation offers a method for functionalizing N₂ using organometallic chemistry.

Iron Complexes of a Proton-Responsive SCS Pincer Ligand with a Sensitive Electronic Structure

Sulfur/carbon/sulfur pincer ligands have an interesting combination of strong-field and weak-field donors, a coordination environment that is also present in the nitrogenase active site. Here, we explore the electronic structures of iron(II) and iron(III) complexes with such a pincer ligand, bearing a monodentate phosphine, thiolate S donor, amide N donor, ammonia, or CO. The ligand scaffold features a proton-responsive thioamide site, and the protonation state of the ligand greatly influences the reduction potential of iron in the phosphine complex. The N-H bond dissociation free energy, derived from the Bordwell equation, is 56 ± 2 kcal/mol. Electron paramagnetic resonance (EPR) spectroscopy and superconducting quantum interference device (SQUID) magnetometry measurements show that the iron(III) complexes with S and N as the fourth donors have an intermediate spin (S = 3/2) ground state with a large zero field splitting, and X-ray absorption spectra show a high Fe-S covalency. The Mössbauer spectrum changes drastically with the position of a nearby alkali metal cation in the iron(III) amido complex, and density functional theory calculations explain this phenomenon through a change between having the doubly occupied orbital as dz2 or dyz, as the former is more influenced by the nearby positive charge.

Accurate 57Fe Mössbauer Parameters from General Gaussian Basis Sets

The prediction of isomer shifts in ⁵⁷Fe Mossbauer spectra is typically achieved by building calibration lines using the values of the density at the nuclear position. Using Slater-type orbital basis or large and specific Gaussian-type orbital basis has been thus far mandatory to achieve accurate predictions with density functional theory methods. In this work, we show that replacing the value of the density at the nucleus by the density integrated in a sphere of radius 0.06 au centered on the Fe nuclei yields excellent calibration lines (r² = 0.976) with a high predictive power (q² = 0.975, MAE = 0.055 mm·s^–1) while using the conventional def2-TZVP basis set and X-ray geometrical parameters. Our data set comprises 69 ⁵⁷Fe-containing compounds and 103 signals. We also find B3LYP performing significantly better than the PW91 functional.

Electrochemically switchable polymerization from surface-anchored molecular catalysts

Redox-switchable polymerizations of lactide and epoxides were extended to the solid state by anchoring an iron-based polymerization catalyst to TiO2 nanoparticles. The reactivity of the molecular complexes and their redox-switching characteristics were maintained in the solid-state. These properties resulted in surface-initiated polymerization reactions that produced polymer brushes whose chemical composition is dictated by the oxidation state of the iron-based complex. Depositing the catalyst-functionalized TiO2 nanoparticles on fluorine-doped tin oxide resulted in an electrically addressable surface that could be used to demonstrate spatial control in redox-switchable polymerization reactions. By using a substrate that contained two electrically isolated domains wherein one domain was exposed to an oxidizing potential, patterns of surface-bound polyesters and polyethers were accessible through sequential application of lactide and cyclohexene oxide. The differentially functionalized surfaces demonstrated distinct physical properties that illustrated the promise for using the method to pattern surfaces with multiple, chemically distinct polymer brushes.

57Fe Mössbauer parameters from domain based local pair-natural orbital coupled-cluster theory

We report on applications of the domain based local pair-natural orbital (PNO) coupled-cluster method within the singles and doubles approximation (DLPNO-CCSD) to the calculation of ⁵⁷Fe isomer shifts and quadrupole splittings in a small training set of iron complexes consisting of large molecular ligands and iron atoms in varying charge, spin, and oxidation states. The electron densities and electric field gradients needed for these calculations were obtained within the recently implemented analytic derivative scheme. A method for the direct treatment of scalar relativistic effects in the calculation of effective electron densities is described by using the first-order Douglas-Kroll-Hess Hamiltonian and a Gaussian charge distribution model for the nucleus. The performance of DLPNO-CCSD is compared with four modern-day density functionals, namely, RPBE, TPSS, B3LYP, and B2PLYP, as well as with the second-order Møller-Plesset perturbation theory. An excellent correlation between the calculated electron densities and the experimental isomer shifts is attained with the DLPNO-CCSD method. The correlation constant a obtained from the slope of the linear correlation plot is found to be ≈-0.31 a.u.³ mm s^-1, which agrees very well with the experimental calibration constant α = -0.31 ± 0.04 a.u.³ mm s^-1. This value of a is obtained consistently using both nonrelativistic and scalar relativistic DLPNO-CCSD electron densities. While the B3LYP and B2PLYP functionals achieve equally good correlation between theory and experiment, the correlation constant a is found to deviate from the experimental value. Similar trends are observed also for quadrupole splittings. The value of the nuclear quadrupole moment for ⁵⁷Fe is estimated to be 0.15 b at the DLPNO-CCSD level. This is consistent with previous results and is here supported by a higher level of theory. The DLPNO-CCSD results are found to be insensitive to the intrinsic approximations in the method, in particular the PNO occupation number truncation error, while the results obtained with density functional theory (DFT) are found to depend on the choice of the functional. In a statistical sense, i.e., on the basis of the linear regression analysis, however, the accuracies of the DFT and DLPNO-CCSD results can be considered comparable.

The Spectroscopy of Nitrogenases

Nitrogenases are responsible for biological nitrogen fixation, a crucial step in the biogeochemical nitrogen cycle. These enzymes utilize a two-component protein system and a series of iron-sulfur clusters to perform this reaction, culminating at the FeMco active site (M = Mo, V, Fe), which is capable of binding and reducing N2 to 2NH3. In this review, we summarize how different spectroscopic approaches have shed light on various aspects of these enzymes, including their structure, mechanism, alternative reactivity, and maturation. Synthetic model chemistry and theory have also played significant roles in developing our present understanding of these systems and are discussed in the context of their contributions to interpreting the nature of nitrogenases. Despite years of significant progress, there is still much to be learned from these enzymes through spectroscopic means, and we highlight where further spectroscopic investigations are needed.

Epoxidation Catalyzed by the Nonheme Iron(II)- and 2-Oxoglutarate-Dependent Oxygenase, AsqJ: Mechanistic Elucidation of Oxygen Atom Transfer by a Ferryl Intermediate

Mechanisms of enzymatic epoxidation via oxygen atom transfer (OAT) to an olefin moiety is mainly derived from the studies on thiolate-heme containing epoxidases, such as cytochrome P450 epoxidases. The molecular basis of epoxidation catalyzed by nonheme-iron enzymes is much less explored. Herein, we present a detailed study on epoxidation catalyzed by the nonheme iron(II)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenase, AsqJ. The native substrate and analogues with different para substituents ranging from electron-donating groups (e.g., methoxy) to electron-withdrawing groups (e.g., trifluoromethyl) were used to probe the mechanism. The results derived from transient-state enzyme kinetics, Mössbauer spectroscopy, reaction product analysis, X-ray crystallography, density functional theory calculations, and molecular dynamic simulations collectively revealed the following mechanistic insights: (1) The rapid O2 addition to the AsqJ Fe(II) center occurs with the iron-bound 2OG adopting an online-binding mode in which the C1 carboxylate group of 2OG is trans to the proximal histidine (His134) of the 2-His-1-carboxylate facial triad, instead of assuming the offline-binding mode with the C1 carboxylate group trans to the distal histidine (His211); (2) The decay rate constant of the ferryl intermediate is not strongly affected by the nature of the para substituents of the substrate during the OAT step, a reactivity behavior that is drastically different from nonheme Fe(IV)-oxo synthetic model complexes; (3) The OAT step most likely proceeds through a stepwise process with the initial formation of a C(benzylic)-O bond to generate an Fe-alkoxide species, which is observed in the AsqJ crystal structure. The subsequent C3-O bond formation completes the epoxide installation.

The influences of carbon donor ligands on biomimetic multi-iron complexes for N2 reduction

The active site clusters of nitrogenase enzymes possess the only examples of carbides in biology. These are the only biological FeS clusters that are capable of reducing N₂ to NH₄⁺, implicating the central carbon and its interaction with Fe as important in the mechanism of N₂ reduction. This biological question motivates study of the influence of carbon donors on the electronic structure and reactivity of unsaturated, high-spin iron centers. Here, we present functional and structural models that test the impacts of carbon donors and sulfide donors in simpler iron compounds. We report the first example of a diiron complex that is bridged by an alkylidene and a sulfide, which serves as a high-fidelity structural and spectroscopic model of a two-iron portion of the active-site cluster (FeMoco) in the resting state of Mo-nitrogenase. The model complexes have antiferromagnetically coupled pairs of high-spin iron centers, and sulfur K-edge X-ray absorption spectroscopy shows comparable covalency of the sulfide for C and S bridged species. The sulfur-bridged compound does not interact with N₂ even upon reduction, but upon removal of the sulfide it becomes capable of reducing N₂ to NH₄⁺ with the addition of protons and electrons. This provides synthetic support for sulfide extrusion in the activation of nitrogenase cofactors.

High-spin diiron alkylidenes give insight into the electronic structure and functional relevance of carbon in the FeMoco active site of nitrogenase.

Mössbauer isomer shifts and effective contact densities obtained by the exact two-component (X2C) relativistic method and its local variants

A standalone program to calculate scalar relativistic effective contact densities.

Masked Radicals: Iron Complexes of Trityl, Benzophenone, and Phenylacetylene

We report the first Fe─CPh₃ complex, and show that the long Fe─C bond can be disrupted by neutral π-acceptor ligands (benzophenone and phenylacetylene) to release the triphenylmethyl radical. The products are formally iron(I) complexes, but X-ray absorption spectroscopy coupled with density functional and multireference ab initio calculations indicates that the best description of all the complexes is iron(II). In the formally iron(I) complexes, this does not imply that the π-acceptor ligand has radical character, because the iron(II) description arises from doubly-occupied frontier molecular orbitals that are shared equitably by the iron and the π-acceptor ligand, and the unpaired electrons lie on the metal. Despite the lack of substantial radical character on the ligands, alkyne and ketone fragments can couple to form a high-spin iron(III) complex with a cyclized metalladihydrofuran core.

Nitrogenase-Relevant Reactivity of a Synthetic Iron-Sulfur-Carbon Site

Simple synthetic compounds with only S and C donors offer a ligation environment similar to the active site of nitrogenase (FeMoco) and thus demonstrate reasonable mechanisms and geometries for N₂ binding and reduction in nature. We recently reported the first example of N₂ binding at a mononuclear iron site supported by only S and C donors. In this work, we report experiments that examine the mechanism of N₂ binding in this system. The reduction of an iron(II) tris(thiolate) complex with 1 equiv of KC₈ leads to a thermally unstable intermediate, and a combination of Mössbauer, EPR, and X-ray absorption spectroscopies identifies it as a high-spin (S = 3/2) iron(I) species that maintains coordination of all three sulfur atoms. DFT calculations suggest that this iron(I) intermediate has a pseudotetrahedral geometry that resembles the S₃C iron coordination environment of the belt iron sites in the resting state of the FeMoco. Further reduction to the iron(0) oxidation level under argon causes the dissociation of one of the thiolate donors and gives an η⁶-arene species which reacts with N₂. Thus, in this system the loss of thiolate and binding of N₂ require reduction beyond the iron(I) level to the iron(0) level. Further reduction of the iron(0)-N₂ complex gives a reactive, formally iron(-I) species. Treatment of the putative iron(-I) complex with weak acids gives low yields of ammonia and hydrazine, demonstrating that these nitrogenase products can be generated from N₂ at a synthetic Fe-S-C site. Catalytic N₂ reduction is not observed, which is attributed to protonation of the supporting ligand and degradation of the complex via ligand dissociation. Identification of the challenges in this system gives insight into the design features needed for functional biomimetic complexes.

Regio- and Diastereoselective Iron-Catalyzed [4+4]-Cycloaddition of 1,3-Dienes

A family of single-component iron precatalysts for the [4+4]-cyclodimerization and intermolecular cross-[4+4]-cycloaddition of monosubstituted 1,3-dienes is described. Cyclooctadiene products were obtained with high regioselectivity, and catalyst-controlled access to either cis- or trans-diastereomers was achieved using 4-substituted diene substrates. Reactions conducted either with single-component precatalysts or with iron dihalide complexes activated in situ proved compatible with common organic functional groups and were applied on multigram scale (up to >100 g). Catalytically relevant, S = 1 iron complexes bearing 2-(imino)pyridine ligands, (^RPI)FeL₂ (^RPI = [2-(2,6-R₂-C₆H₃-N═CMe)-C₅H₄N] where R = ⁱPr or Me, L₂ = bis-olefin), were characterized by single-crystal X-ray diffraction, Mößbauer spectroscopy, magnetic measurements, and DFT calculations. The structural and spectroscopic parameters are consistent with an electronic structure description comprised of a high spin iron(I) center ( S_Fe = 3/2) engaged in antiferromagnetically coupling with a ligand radical anion ( S_PI = -1/2). Mechanistic studies conducted with these single-component precatalysts, including kinetic analyses, ¹²C/¹³C isotope effect measurements, and in situ Mößbauer spectroscopy, support a mechanism involving oxidative cyclization of two dienes that determines regio- and diastereoselectivity. Topographic steric maps derived from crystallographic data provided insights into the basis for the catalyst control through stereoselective oxidative cyclization and subsequent, stereospecific allyl-isomerization and C-C bond-forming reductive elimination.

Effects of N2 Binding Mode on Iron-Based Functionalization of Dinitrogen to Form an Iron(III) Hydrazido Complex

Distinguishing the reactivity differences between N2 complexes having different binding modes is crucial for the design of effective N2-functionalizing reactions. Here, we compare the reactions of a K-bridged, dinuclear FeNNFe complex with a monomeric Fe(N2) complex where the bimetallic core is broken up by the addition of chelating agents. The new anionic iron(0) dinitrogen complex has enhanced electron density at the distal N atoms of coordinated N2, and though the N2 is not as weakened in this monomeric compound, it is much more reactive toward silylation by (CH3)3SiI (TMSI). Double silylation of N2 gives a three-coordinate iron(III) hydrazido(2-) complex, which is finely balanced between coexisting S = 1/2 and S = 3/2 states that are characterized by crystallography, spectroscopy, and computations. These results give insight into the interdependence between binding modes, alkali dependence, reactivity, and magnetic properties within an iron system that functionalizes N2.

Reversible Ligand-Centered Reduction in Low-Coordinate Iron Formazanate Complexes

Coordination of redox-active ligands to metals is a compelling strategy for making reduced complexes more accessible. In this work, we explore the use of redox-active formazanate ligands in low-coordinate iron chemistry. Reduction of an iron(II) precursor occurs at milder potentials than analogous non-redox-active β-diketiminate complexes, and the reduced three-coordinate formazanate-iron compound is characterized in detail. Structural, spectroscopic, and computational analysis show that the formazanate ligand undergoes reversible ligand-centered reduction to form a formazanate radical dianion in the reduced species. The less negative reduction potential of the reduced low-coordinate iron formazanate complex leads to distinctive reactivity with formation of a new N-I bond that is not seen with the β-diketiminate analogue. Thus, the storage of an electron on the supporting ligand changes the redox potential and enhances certain reactivity.

Selective Conversion of CO2 into Isocyanate by Low-Coordinate Iron Complexes

Discovery of the mechanisms for selective transformations of CO2 into organic compounds is a challenge. Herein, we describe the reaction of low-coordinate Fe silylamide complexes with CO2 to give trimethylsilyl isocyanate and the corresponding Fe siloxide complex. Kinetic studies show that this is a two-stage reaction, and the presence of a single equivalent of THF influences the rates of both steps. Isolation of a thermally unstable intermediate provides mechanistic insight that explains both the effect of THF in this reaction, and the way in which the reaction achieves high selectivity for isocyanate formation.

Selective [1,4]-Hydrovinylation of 1,3-Dienes with Unactivated Olefins Enabled by Iron Diimine Catalysts

The selective, intermolecular [1,4]-hydrovinylation of conjugated dienes with unactivated α-olefins catalyzed by α-diimine iron complexes is described. Value-added "skipped" diene products were obtained with exclusive [1,4]-selectivity, and the formation of branched, ( Z)-olefin products was observed with no evidence for alkene isomerization. Mechanistic studies conducted with the well-defined, single-component iron precatalyst (^MesDI)Fe(COD) (^MesDI = [2,4,6-Me₃-C₆H₂-N═CMe]₂); COD = 1,5-cyclooctadiene) provided insights into the origin of the high selectivity. An iron diene complex was identified as the catalyst resting state, and one such isoprene complex, (^iPrDI)Fe(η⁴-C₅H₈), was isolated and characterized. A combination of single crystal X-ray diffraction, Mößbauer spectroscopy, magnetic measurements, and DFT calculations established that the complex is best described as a high-spin Fe(I) center ( S_Fe = 3/2) engaged in antiferromagnetic coupling to an α-diimine radical anion ( S_DI = -1/2), giving rise to the observed S = 1 ground state. Deuterium-labeling experiments and kinetic analyses of the catalytic reaction provided support for a pathway involving oxidative cyclization of an alkene with the diene complex to generate an iron metallacycle. The observed selectivity can be understood in terms of competing steric interactions in the transition states for oxidative cyclization and subsequent β-hydrogen elimination.

Synthesis and Mechanism of Formation of Hydride-Sulfide Complexes of Iron

Iron-sulfide complexes with hydride ligands provide an experimental precedent for spectroscopically detected hydride species on the iron-sulfur MoFe₇S₉C cofactor of nitrogenase. In this contribution, we expand upon our recent synthesis of the first iron sulfide hydride complex from an iron hydride and a sodium thiolate ( Arnet, N. A.; Dugan, T. R.; Menges, F. S.; Mercado, B. Q.; Brennessel, W. W.; Bill, E.; Johnson, M. A.; Holland, P. L., J. Am. Chem. Soc. 2015 , 137 , 13220 - 13223 ). First, we describe the isolation of an analogous iron sulfide hydride with a smaller diketiminate supporting ligand, which benefits from easier preparation of the hydride precursor and easier isolation of the product. Second, we describe mechanistic studies on the C-S bond cleavage through which the iron sulfide hydride product is formed. In a key experiment, use of cyclopropylmethanethiolate as the sulfur precursor leads to products from cyclopropane ring opening, implicating an alkyl radical as an intermediate. Combined with the results of isotopic labeling studies, the data are consistent with a mechanism in which homolytic C-S bond cleavage is followed by rebound of the alkyl radical to abstract a hydrogen atom from iron to give the observed alkane and iron-sulfide products.