On Predicting Mössbauer Parameters of Iron-Containing Molecules with Density-Functional Theory

An Open-Shell FeIV Nitrido

We report the photogeneration and characterization of an open-shell, terminal iron nitrido (EmL)Fe(N) using a sterically encumbered dipyrrin ligand environment. The Fe-N distance in the solid-state, zero-field 57Fe Mössbauer spectrum, and computational analysis are consistent with a triplet electronic ground state of the iron nitrido. Notably, the attenuation of Fe-N multiple bond character through occupying π*Fe-N enables (i) primary C(sp3)-H amination, (ii) H2 cleavage, (iii) aromatic C-C cleavage, and (iv) photocatalytic N-atom transfer reactivity. These modes of reactivity have not previously been observed in low-spin Fe(N) analogues.

The temperature dependence of Mössbauer quadrupole splitting values: a quantum chemical analysis

The two key parameters extracted from Mössbauer spectroscopy, isomer shift and quadrupole splitting, have well-known temperature dependencies. While the behavior of the values following a temperature change has long been known, its microscopic origins are less clear. For quantum chemical calculations - formally representing the situation at 0 K - significant discrepancies with the experiment can arise, especially at elevated temperatures. Herein, we introduce an ensemble-based approach for capturing the temperature dependence of the quadrupole splitting. Our method is exemplified with [Fe(TPP)(2-MeHIm)], an iron(II) high spin system. We rationalise the temperature dependence by changes in the shape of the charge distribution due to vibrational distortion. By using a normal mode fitting approach, we isolated collective nuclear movements associated with the change in the quadrupole splitting.

Stereoelectronic Tuning of Bioinspired Nonheme Iron(IV)-Oxo Species by Amide Groups in Primary and Secondary Coordination Spheres for Selective Oxygenation Reactions

Two mononuclear iron(II) complexes, [(6-amide2-BPMEN)FeII](OTf)2 (1) and [(6-amide-Me-BPMEN)FeII(OTf)](OTf) (2), supported by two BPMEN-derived (BPMEN = N1,N2-dimethyl-N1,N2-bis(pyridine-2-yl-methyl)ethane-1,2-diamine) ligands bearing one or two amide functionalities have been isolated to study their reactivity in the oxygenation of C-H and C═C bonds using isopropyl 2-iodoxybenzoate (iPr-IBX ester) as the oxidant. Both 1 and 2 contain six-coordinate high-spin iron(II) centers in the solid state and in solution. The 6-amide2-BPMEN ligand stabilizes an S = 1 iron(IV)-oxo intermediate, [(6-amide2-BPMEN)FeIV(O)]2+ (1A). The oxidant (1A) oxygenates the C-H and C═C bonds with a high selectivity. Oxidant 1A, upon treatment with 2,6-lutidine, is transformed into another oxidant [{(6-amide2-BPMEN)-(H)}FeIV(O)]+ (1B) through deprotonation of an amide group, resulting in a stronger equatorial ligand field and subsequent stabilization of the triplet ground state. In contrast, no iron-oxo species could be observed from complex 2 and [(6-Me2-BPMEN)FeII(OTf)2] (3) under similar experimental conditions. The iron(IV)-oxo oxidant 1A shows the highest A/K selectivity in cyclohexane oxidation and 3°/2° selectivity in adamantane oxidation reported for any synthetic nonheme iron(IV)-oxo complexes. Theoretical investigation reveals that the hydrogen bonding interaction between the -NH group of the noncoordinating amide group and Fe═O core smears out the equatorial charge density, reducing the triplet-quintet splitting, and thus helping complex 1A to achieve better reactivity.

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.

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.

High-Pressure Polymorphism in Silver Ferrite Delafossite, AgFeO2

The delafossites are a class of layered metal oxides that are notable for being able to exhibit optical transparency alongside an in-plane electrical conductivity, making them promising platforms for the development of transparent conductive oxides. Pressure-induced polymorphism offers a direct method for altering the electrical and optical properties in this class, and although the copper delafossites have been studied extensively under pressure, the silver delafossites remain only partially studied. We report two new high-pressure polymorphs of silver ferrite delafossite, AgFeO2, that are stabilized above ∼6 and ∼14 GPa. In situ X-ray diffraction and vibrational spectroscopy measurements are used to examine the structural changes across the two phase transitions. The high-pressure structure between 6 and 14 GPa is assigned as a monoclinic C2/c structure that is analogous to the high-pressure phase reported for AgGaO2. Nuclear resonant forward scattering reveals no change in the spin state or valence state at the Fe3+ site up to 15.3(5) GPa.

FeN4 Environments upon Reduction: A Computational Analysis of Spin States, Spectroscopic Properties, and Active Species

FeN4 motifs, found, for instance, in bioinorganic chemistry as heme-type cofactors, play a crucial role in man-made FeNC catalysts for the oxygen reduction reaction. Such single-atom catalysts are a potential alternative to platinum-based catalysts in fuel cells. Since FeNC catalysts are prepared via pyrolysis, the resulting materials are amorphous and contain side phases and impurities. Therefore, the geometric and electronic nature of the catalytically active FeN4 site remains to be clarified. To further understand the behavior of FeN4 centers in electrochemistry and their expected spectroscopic behavior upon reduction, we investigate two FeN4 environments (pyrrolic and pyridinic). These are represented by the model complexes [Fe(TPP)Cl] and [Fe(phen2N2)Cl], where TPP = tetraphenylporphyrin and phen = 1,10-phenanthroline. We predict their Mössbauer, UV-vis, and NRV spectral data using density functional theory as windows into their electronic structure differences. By varying the axial ligand, we further show how well small chemical changes in both complexes can be discerned. We find that the differences in ligand field strength in pyrrolic and pyridinic coordination result in different spin ground states, which in turn leads to distinct Mössbauer spectroscopic properties. As a result, pyrrolic nitrogen donors with a weaker ligand field are predicted to show more pronounced spectroscopic differences under in situ and operando conditions, while pyridinic nitrogen donors are expected to show less pronounced spectroscopic changes upon reduction and/or ligand loss. We therefore suggest that a weaker ligand field leads to better detectability of catalytic intermediates in in situ and operando experiments.

Spin-state versatility in FeII4L6 supramolecular cages with a pyridyl-hydrazone ligand scaffold modulated by solvents and counter anions

Discrete spin crossover (SCO) tetranuclear cages are a unique class of materials that have potential use in next-generation molecular recognition and sensing. In this work, two new edge-bridged SCO FeII4L6 (L = 2,7-bis(((E)-pyridin-2-ylmethylene)amino)benzo[lmn] [3,8]phenanthroline-1,3,6,8(2H,7H)-tetraone) supramolecular cages with different counter anions: ClO4- (2) and CF3SO3- (3) were constructed via subcomponent self-assembly to investigate both solvent and anion influences on their magnetic properties and compare them to cage 1 with a BF4- anion. Pyridyl-hydrazone bidentate ligand scaffolds were employed to replace the 'classical' imidazole/thiazolyl-imine coordination units to induce SCO behaviour in these cages. 2 and 3 were structurally characterized by single-crystal X-ray diffraction analysis and electrospray ionization time-of-flight mass spectrometry. Magnetic susceptibilities of 1-3 and 1-3·desolvated indicate that the solvents' presence is in favor of the low-spin (LS) state. While different counter anions in 1-3·desolvated affect the spin-state configurations of the four FeII metal centers. According to the 57Fe Mössbauer spectral analysis, the spin-state distributions in 1-3 at 80 K are [2 high-spin (HS)-2LS], [1HS-3LS] and [2HS-2LS], respectively and density functional theory calculations were employed to investigate the reasons. These findings provide insights to regulate the spin-state versatility of SCO FeII cage systems in the solid state.

Benchmarks for transition metal spin-state energetics: why and how to employ experimental reference data?

Accurate prediction of energy differences between alternative spin states of transition metal complexes is essential in computational (bio)inorganic chemistry-for example, in characterization of spin crossover materials and in the theoretical modeling of open-shell reaction mechanisms-but it remains one of the most compelling problems for quantum chemistry methods. A part of this challenge is to obtain reliable reference data for benchmark studies, as even the highest-level applicable methods are known to give divergent results. This Perspective discusses two possible approaches to method benchmarking for spin-state energetics: using either theoretically computed or experiment-derived reference data. With the focus on the latter approach, an extensive general review is provided for the available experimental data of spin-state energetics and their interpretations in the context of benchmark studies, targeting the possibility of back-correcting the vibrational effects and the influence of solvents or crystalline environments. With a growing amount of experience, these effects can be now not only qualitatively understood, but also quantitatively modeled, providing the way to derive nearly chemically accurate estimates of the electronic spin-state gaps to be used as benchmarks and advancing our understanding of the phenomena related to spin states in condensed phases.

Precise Fingerprint Determination of Vibrational Infrared Spectra in a Series of Co(II) Clathrochelates through Experimental and Theoretical Analyses

The energetic demands of modern society for clean energy vectors, such as H2, have caused a surge in research associated with homogeneous and immobilized electrocatalysts that may replace Pt. In particular, clathrochelates have shown excellent electrocatalytic properties for the hydrogen evolution reaction (HER). However, the actual mechanism for the HER catalyzed by these d-metal complexes remains an open debate, which may be addressed via Operando spectroelectrochemistry. The prediction of electrochemical properties via density functional theory (DFT) needs access to thermodynamic functions, which are only available after Hessian calculations. Unfortunately, there is a notable lack in the current literature regarding the precise evaluation of vibrational spectra of such complexes, given their structural complexity and the associated tangled IR spectra. In this work, we have performed a detailed theoretical and experimental analysis in a family of Co(II) clathrochelates, in order to establish univocally their IR pattern, and also the calculation methodology that is adequate for such predictions. In summary, we have observed the presence of multiple common bands shared by this clathrochelate family, using the B3LYP functional, the LANL2DZ basis, and effective core potentials (ECP) for heavy atoms. The most important issue addressed in this article was therefore related to the detailed assignment of the fingerprint associated with cobalt(II) clathrochelates, which is a challenging endeavor due to the crowded nature of their spectra.

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.

Thermally Induced Solid-Phase Quasi-Intramolecular Redox Reactions of [Hexakis(urea-O)iron(III)] Permanganate: An Easy Reaction Route to Prepare Potential (Fe,Mn)Ox Catalysts for CO2 Hydrogenation

Research on new reaction routes and precursors to prepare catalysts for CO2 hydrogenation has enormous importance. Here, we report on the preparation of the permanganate salt of the urea-coordinated iron(III), [hexakis(urea-O)iron(III)]permanganate ([Fe(urea-O)6](MnO4)3) via an affordable synthesis route and preliminarily demonstrate the catalytic activity of its (Fe,Mn)Ox thermal decomposition products in CO2 hydrogenation. [Fe(urea-O)6](MnO4)3 contains O-coordinated urea ligands in octahedral propeller-like arrangement around the Fe3+ cation. There are extended hydrogen bond interactions between the permanganate ions and the hydrogen atoms of the urea ligands. These hydrogen bonds serve as reaction centers and have unique roles in the solid-phase quasi-intramolecular redox reaction of the urea ligand and the permanganate anion below the temperature of ligand loss of the complex cation. The decomposition mechanism of the urea ligand (ammonia elimination with the formation of isocyanuric acid and biuret) has been clarified. In an inert atmosphere, the final thermal decomposition product was manganese-containing wuestite, (Fe,Mn)O, at 800 °C, whereas in ambient air, two types of bixbyite (Fe,Mn)2O3 as well as jacobsite (Fe,Mn)T-4(Fe,Mn)OC-62O4), with overall Fe to Mn stoichiometry of 1:3, were formed. These final products were obtained regardless of the different atmospheres applied during thermal treatments up to 350 °C. Disordered bixbyite formed first with inhomogeneous Fe and Mn distribution and double-size supercell and then transformed gradually into common bixbyite with regular structure (and with 1:3 Fe to Mn ratio) upon increasing the temperature and heating time. The (Fe,Mn)Ox intermediates formed under various conditions showed catalytic effect in the CO2 hydrogenation reaction with <57.6% CO2 conversions and <39.3% hydrocarbon yields. As a mild solid-phase oxidant, hexakis(urea-O)iron(III) permanganate, was found to be selective in the transformation of (un)substituted benzylic alcohols into benzaldehydes and benzonitriles.

Versatile Fe-Sn Bonding Interactions in a Metallostannylene System: Multiple Bonding and C-H Bond Activation

The metallostannylene Cp*(iPr2MeP)(H)2Fe-SnDMP (1; Cp* = η5-C5Me5; DMP = 2,6-dimesitylphenyl), formed by hydrogen migration in a putative Cp*(iPr2MeP)HFe[Sn(H)DMP] intermediate, serves as a robust platform for exploration of transition-metal main-group element bonding and reactivity. Upon one-electron oxidation, 1 expels H2 to generate the coordinatively unsaturated [Cp*(iPr2MeP)Fe═SnDMP][B(C6F5)4] (3), which possesses a highly polarized Fe-Sn multiple bond that involves interaction of the tin lone pair with iron. Evidence from EPR and 57Fe Mössbauer spectroscopy, along with DFT studies, shows that 3 is primarily an iron-based radical with charge localization at tin. Upon reduction of 3, C-H bond activation of the phosphine ligand was observed to produce Cp*HFe(κ2-(P,Sn)═Sn(DMP)CH2CHMePMeiPr) (5). Complex 5 was also accessed via thermolysis of 1, and kinetics studies of this thermolytic pathway indicate that the reductive elimination of H2 from 1 to produce a stannylyne intermediate, Cp*(iPr2MeP)Fe[SnDMP] (A), is likely rate-determining. Evidence indicates that the production of 5 proceeds through a concerted C-H bond activation. DFT investigations suggest that the transition state for this transformation involves C-H cleavage across the Fe-Sn bond and that a related transition state where C-H bond activation occurs exclusively at the tin center is disfavored, illustrating an effect of iron-tin cooperativity in this system.

A Reactive, Photogenerated High-Spin (S = 2) FeIV(O) Complex via O2 Activation

Addition of dioxygen at low temperature to the non-heme ferrous complex Fe^II(Me₃TACN)((OSi^Ph2)₂O) (1) in 2-MeTHF produces a peroxo-bridged diferric complex Fe₂^III(μ-O₂)(Me₃TACN)₂((OSi^Ph2)₂O)₂ (2), which was characterized by UV-vis, resonance Raman, and variable field Mössbauer spectroscopies. Illumination of a frozen solution of 2 in THF with white light leads to homolytic O-O bond cleavage and generation of a Fe^IV(O) complex 4 (ν(Fe=O) = 818 cm^-1; δ = 0.22 mm s^-1, ΔE_Q = 0.23 mm s^-1). Variable field Mössbauer spectroscopy measurements show that 4 is a rare example of a high-spin S = 2 Fe^IV(O) complex and the first synthetic example to be generated directly from O₂. Complex 4 is highly reactive, as expected for a high-spin ferryl, and decays rapidly in fluid solution at cryogenic temperatures. This decay process in 2-MeTHF involves C-H cleavage of the solvent. However, the controlled photolysis of 2 in situ with visible light and excess phenol substrate leads to competitive phenol oxidation, via the proposed transient generation of 4 as the active oxidant.

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.

57Fe Mössbauer isomer shift of pure iron and iron oxides at high pressure-An experimental and theoretical study

The ⁵⁷Fe isomer shift (IS) of pure iron has been measured up to 100 GPa using synchrotron Mössbauer spectroscopy in the time domain. Apart from the expected discontinuity due to the α → ε structural and spin transitions, the IS decreases monotonically with increasing pressure. The absolute shifts were reproduced without semi-empirical calibrations by periodic density functional calculations employing extensive localized basis sets with several common density functionals. However, the best numerical agreement is obtained with the B1WC hybrid functional. Extension of the calculations to 350 GPa, a pressure corresponding to the Earth's inner core, predicted the IS range of 0.00 to -0.85 mm/s, covering the span from Fe(0) to Fe(VI) compounds measured at ambient pressure. The calculations also reproduced the pressure trend from polymorphs of prototypical iron oxide minerals, FeO and Fe₂O₃. Analysis of the electronic structure shows a strong donation of electrons from oxygen to iron at high pressure. The assignment of formal oxidation to the Fe atom becomes ambiguous under this condition.

Trajectory surface-hopping photoinduced dynamics from Rydberg states of trimethylamine

We present a computational study on nonadiabatic excited-state dynamics initiated from the 3p Rydberg states of trimethylamine (TMA). We utilise a methodology based on full-dimensional (39 D) trajectory surface-hopping (TSH) simulations, in which propagation is carried out on on-the-fly density functional theory (DFT)/time-dependent DFT (TD-DFT) potentials. Both our electronic structure benchmarks to high-level ab initio methods (EOM-CCSD, CASPT2) and TSH simulations demonstrate high-accuracy of the applied CAM-B3LYP functional for the description of Rydberg excited states. Based on our excited-state simulations, we construct the following mechanistic picture: when pumped resonantly to the 3p Rydberg manifold, TMA coherently vibrates along the planarisation mode with a period of 104 fs and an exponential coherence decay time constant of 240 fs. Nonadiabatic dynamics occur on a faster (∼1 ps) and a slower (∼3 ps) timescale, along the N-C stretching mode by mixing with a dissociative σN-C* state. As a minor relaxation channel, 3p → 3s internal conversion occurs via branching at the σN-C*/3s intersection. We find that photodissociaton is hardly observable within 3 ps (1%), which is a failure of the range-separated hybrid CAM-B3LYP functional, as a consequence of its static electron correlation deficiency at long range. In contrast, pure DFT (GGA-BLYP) provides an accurate long-range description (19% dissociation yield), also supported by comparison to recent ultrafast experiments, even if the Rydberg state energies are significantly underestimated (>1 eV). Finally, we reveal the crucial role of vibrational coherence and energy transfer from the planarisation mode for N-C bond activation and resulting nonadiabatic dynamics. The present work illustrates the importance of nuclear-electronic coupling for excited-state dynamics and branching at conical intersections.

Proton-Coupled Electron-Transfer Reactivity Controls Iron versus Sulfur Oxidation in Nonheme Iron-Thiolate Complexes

Reaction of the five-coordinate Fe^II(N₄S) complexes, [Fe^II(iPr₃TACN)(abt^X)](OTf) (abt = aminobenzenethiolate, X = H, CF₃), with a one-electron oxidant and an appropriate base leads to net H atom loss, generating new Fe^III(iminobenzenethiolate) complexes that were characterized by single-crystal X-ray diffraction (XRD), as well as UV-vis, EPR, and Mössbauer spectroscopies. The spectroscopic data indicate that the iminobenzenethiolate complexes have S = 3/2 ground states. In the absence of a base, oxidation of the Fe^II(abt) complexes leads to disulfide formation instead of oxidation at the metal center. Bracketing studies with separated proton-coupled electron-transfer (PCET) reagents show that the Fe^II(aminobenzenethiolate) and Fe^III(iminobenzenethiolate) forms are readily interconvertible by H⁺/e^- transfer and provide a measure of the bond dissociation free energy (BDFE) for the coordinated N-H bond between 64 and 69 kcal mol^-1. This work shows that coordination to the iron center causes a dramatic weakening of the N-H bond and that Fe- versus S-oxidation in a nonheme iron complex can be controlled by the protonation state of an ancillary amino donor.

Accurate Prediction of Mössbauer Hyperfine Parameters in Bis-Axially Coordinated Iron(II) Phthalocyanines Using Density Functional Theory Calculations: A Story of a Single Orbital Revealed by Natural Bond Orbital Analysis

Density Functional Theory (DFT) calculations coupled with several exchange-correlation functionals were used for the prediction of Mössbauer hyperfine parameters of 36 bis-axially coordinated iron(II) phthalocyanine complexes with the general formulas PcFeL₂, PcFeL'L″, and [PcFeX₂]^2-, including four new compounds. Both gas-phase and PCM calculations using BPW91 and MN12L exchange-correlation functionals were found to accurately predict both Mössbauer quadrupole splittings and the correct trends in experimentally observed isomer shifts. In comparison, hybrid exchange-correlation functionals underestimated quadrupole splittings, while still accurately predicted isomer shifts. Out of ∼40 exchange-correlation functionals tested, only MN12L was found to correctly reproduce quadrupole splitting trends in the PcFeL₂ complexes coordinated with phosphorus-donor axial ligands (i.e., P(OnBu)₃ ≈ P(OEt)₃ < PMe₃ < P[(CH₂O)₂CH₂]-p-C₆H₄NO₂ < PEt₃ ≈ PnBu₃). Natural Bond Orbital (NBO) analysis was successfully used to explain the general trends in the observed quadrupole splitting for all compounds of interest. In particular, the general trends in the quadrupole splitting correlate well with the axial ligand dependent, NBO-predicted population of the 3d_z2 orbital of the Fe ion and are reflective of the hypothesis proposed by Ohya and co-workers ( Inorg. Chem., 1984, 23, 1303) on the adaptability of the phthalocyanine's π-system toward Fe-L_ax interactions. The first X-ray crystal structure of a PcFeL₂ complex with axial phosphine ligands is also reported.

Iron-catalysed synthesis and chemical recycling of telechelic 1,3-enchained oligocyclobutanes

Closed-loop recycling offers the opportunity to mitigate plastic waste through reversible polymer construction and deconstruction. Although examples of chemical recycling of polymers are known, few have been applied to materials derived from abundant commodity olefinic monomers, which are the building blocks of ubiquitous plastic resins. Here we describe a [2+2] cycloaddition/oligomerization of 1,3-butadiene to yield a previously unrealized telechelic microstructure of (1,n'-divinyl)oligocyclobutane. This material is thermally stable, has stereoregular segments arising from chain-end control, and exhibits high crystallinity even at low molecular weight. Exposure of the oligocyclobutane to vacuum in the presence of the pyridine(diimine) iron precatalyst used to synthesize it resulted in deoligomerization to generate pristine butadiene, demonstrating a rare example of closed-loop chemical recycling of an oligomeric material derived from a commodity hydrocarbon feedstock.

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.

Spin-crossover behavior of bis[dihydrobis(4-methylpyrazol-1-yl-borate)]-(2,2'-bipyridine)iron and analogous complexes in the bulk and in thin films: Elucidating the influence of π-π-interactions on the type of spin transition

The new complex [Fe(H₂B(4-CH₃-pz)₂)₂(bipy)] (1, pz  =  pyrazole, bipy  =  2,2'-bipyridine) is synthesized and investigated by temperature-dependent magnetic measurements, Mößbauer, electronic absorption and vibrational spectroscopy as well as single crystal x-ray diffraction. In the crystal structure of 1 the complexes are pairwise linked to dimers by intermolecular π-π interactions between their bipyridine ligands, with a shortest intradimer distance between two neighboring pyridine rings of 3.575 Å. Analysis of the crystal structures of related iron(II) bis(dihydrobis(pyrazoyl)borate) complexes reveals that most of them contain similar dimers, and that at short π-π intra-dimer distances the complexes are locked in the high-spin state whereas at long distances complete thermal spin crossover (SCO) is observed. Compound 1 with an intermediate π-π intra-dimer distance shows incomplete SCO in the bulk but complete SCO in vacuum-deposited thin films where intermolecular interactions are absent. The implications of this remarkable structure-property relationship are discussed.

Bis(imino)pyrazine-Supported Iron Complexes: Ligand-Based Redox Chemistry, Dearomatization, and Reversible C-C Bond Formation

Iron complexes supported by novel π-acidic bis(imino)pyrazine (P^PzDI) ligands can be functionalized at the nonligated nitrogen atom, and this has a marked effect on the redox properties of the resulting complexes. Dearomatization is observed in the presence of cobaltocene, which reversibly reduces the pyrazine core and not the imine functionality, as observed in the case of the pyridinediimine-ligated iron analogues. The resulting ligand-based radical is prone to dimerization through the formation of a long carbon-carbon bond, which can be subsequently cleaved under mild oxidative conditions.

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.

Activation of Dioxygen by a Mononuclear Nonheme Iron Complex: Sequential Peroxo, Oxo, and Hydroxo Intermediates

The activation of dioxygen by Fe^II(Me₃TACN)(S₂SiMe₂) (1) is reported. Reaction of 1 with O₂ at -135 °C in 2-MeTHF generates a thiolate-ligated (peroxo)diiron complex Fe^III₂(O₂)(Me₃TACN)₂(S₂SiMe₂)₂ (2) that was characterized by UV-vis (λ_max = 300, 390, 530, 723 nm), Mössbauer (δ = 0.53, |ΔE_Q| = 0.76 mm s^-1), resonance Raman (RR) (ν(O-O) = 849 cm^-1), and X-ray absorption (XAS) spectroscopies. Complex 2 is distinct from the outer-sphere oxidation product 1^ox (UV-vis (λ_max = 435, 520, 600 nm), Mössbauer (δ = 0.45, |ΔE_Q| = 3.6 mm s^-1), and EPR (S = 5/2, g = [6.38, 5.53, 1.99])), obtained by one-electron oxidation of 1. Cleavage of the peroxo O-O bond can be initiated either photochemically or thermally to produce a new species assigned as an Fe^IV(O) complex, Fe^IV(O)(Me₃TACN)(S₂SiMe₂) (3), which was identified by UV-vis (λ_max = 385, 460, 890 nm), Mössbauer (δ = 0.21, |ΔE_Q| = 1.57 mm s^-1), RR (ν(Fe^IV═O) = 735 cm^-1), and X-ray absorption spectroscopies, as well as reactivity patterns. Reaction of 3 at low temperature with H atom donors gives a new species, Fe^III(OH)(Me₃TACN)(S₂SiMe₂) (4). Complex 4 was independently synthesized from 1 by the stoichiometric addition of a one-electron oxidant and a hydroxide source. This work provides a rare example of dioxygen activation at a mononuclear nonheme iron(II) complex that produces both Fe^III-O-O-Fe^III and Fe^IV(O) species in the same reaction with O₂. It also demonstrates the feasibility of forming Fe/O₂ intermediates with strongly donating sulfur ligands while avoiding immediate sulfur oxidation.

Mononuclear, Nonheme, High-Spin {FeNO}7/8 Complexes Supported by a Sterically Encumbered N4S-Thioether Ligand

The synthesis of a new nonheme iron NO binding complex, [Fe^II(CH₃CN)(N3Py^2PhSEtCN)](BF₄)₂ (1), is reported. Complex 1, which contains two sterically encumbering phenyl substituents, exhibits a high-spin (hs) Fe^II (S = 2) ground state in contrast to the S = 0 ground state for unsubstituted [Fe^II(CH₃CN)(N3PySEtCN)(BF₄)₂. Reaction of 1 with NO_(g) in CH₃CN yields an {FeNO}⁷ (S = ³/₂) complex 2, which slowly decays at 25 °C with loss of NO^• to regenerate 1. One-electron reduction of 2 with Cr(C₆H₆)₂ at -40 °C yields the metastable, S = 1 {FeNO}⁸ complex 3. The nitrosyl moieties in thioether-ligated 2 and 3 are significantly less activated than in thiolate-ligated [Fe(NO)(N3PyS)]^+/0, a structurally analogous pair of hs {FeNO}^7/8 complexes. Calculations reveal that reduction of 2 is iron-centered, which may be a general property of hs {FeNO}^7/8 complexes.

A Mononuclear, Nonheme FeII-Piloty's Acid (PhSO2NHOH) Adduct: An Intermediate in the Production of {FeNO}7/8 Complexes from Piloty's Acid

Reaction of the mononuclear nonheme complex [Fe^II(CH₃CN)(N3PyS)]BF₄ (1) with an HNO donor, Piloty's acid (PhSO₂NHOH, P.A.), at low temperature affords a high-spin ( S = 2) Fe^II-P.A. intermediate (2), characterized by ⁵⁷Fe Mössbauer and Fe K-edge X-ray absorption (XAS) spectroscopies, with interpretation of both supported by DFT calculations. The combined methods indicate that P.A. anion binds as the N-deprotonated tautomer (PhSO₂NOH^-) to [Fe^II(N3PyS)]⁺, leading to 2. Complex 2 is the first spectroscopically characterized example, to our knowledge, of P.A. anion bound to a redox-active metal center. Warming of 2 above -60 °C yields the stable {FeNO}⁷ complex [Fe(NO)(N3PyS)]BF₄ (4), as evidenced by ¹H NMR, ATR-IR, and Mössbauer spectroscopies. Isotope labeling experiments with ¹⁵N-labeled P.A. confirm that the nitrosyl ligand in 4 derives from P.A. In contrast, addition of a second equivalent of a strong base leads to S-N cleavage and production of an {FeNO}⁸ species, the deprotonated analog of an Fe-HNO complex. This work has implications for the targeted delivery of HNO/NO^-/NO· to nonheme Fe centers in biological and synthetic applications, and suggests a new role for nonheme Fe^II complexes in the assisted degradation of HNO donor molecules.

Observation of the First Spin Crossover in an Iron(II) Complex with an S6 Coordination Environment: Tris[bis(N,N-diethylamino)carbeniumdithiocarboxylato]iron(II) Hexafluorophosphate

For the first time, the spin-crossover (SCO) phenomenon has been observed in an Fe^II -S₆ system in a tris(chelate)-type iron(II) complex with a zwitterionic sulfur donor bidentate, bis(N,N-diethylamino)carbeniumdithiocarboxylate (EtL), [Fe^II (EtL)₃ ](PF₆ )₂ (1), as synthesized by the reaction of a precursor complex [Fe^II (CH₃ CN)₆ ](PF₆ )₂ with EtL. In the solid state, the high-spin (HS) d⁶ state at ambient temperature and the low-spin (LS) d⁶ state at temperatures lower than approximately 240 K were evidenced by magnetic measurements with SQUID and Mössbauer spectra in the temperature range 4-290 K. X-ray analyses of the crystals at various temperatures disclosed that the distorted trigonal prismatic coordination environments essentially do not change; however, contraction of Fe-S distances by approximately 10 % (0.22 Å), ordering of alkyl groups in EtL and PF₆ ^- counteranions, and formation of significant intermolecular S⋅⋅⋅S interactions between adjacent molecules (average distances of 3.59 Å) take place during the transition from the HS to the LS state. A large decrease in the volume of the formula unit (78.1 ų ) might be responsible for the large activation barrier, thereby resulting in a slow phase transition upon cooling.

Structures, Spectroscopic Properties, and Dioxygen Reactivity of 5- and 6-Coordinate Nonheme Iron(II) Complexes: A Combined Enzyme/Model Study of Thiol Dioxygenases

The synthesis of four new FeII(N4S(thiolate)) complexes as models of the thiol dioxygenases are described. They are composed of derivatives of the neutral, tridentate ligand triazacyclononane (R3TACN; R = Me, iPr) and 2-aminobenzenethiolate (abtx; X = H, CF3), a non-native substrate for thiol dioxygenases. The coordination number of these complexes depends on the identity of the TACN derivative, giving 6-coordinate (6-coord) complexes for FeII(Me3TACN)(abtx)(OTf) (1: X = H; 2: X = CF3) and 5-coordinate (5-coord) complexes for [FeII(iPr3TACN)(abtx)](OTf) (3: X = H; 4: X = CF3). Complexes 1-4 were examined by UV-vis, 1H/19F NMR, and Mössbauer spectroscopies, and density functional theory (DFT) calculations were employed to support the data. Mössbauer spectroscopy reveals that the 6-coord 1-2 and 5-coord 3- 4 exhibit distinct spectra, and these data are compared with that for cysteine-bound CDO, helping to clarify the coordination environment of the cys-bound FeII active site. Reaction of 1 or 2 with O2 at -95 °C leads to S-oxygenation of the abt ligand, and in the case of 2, a rare di(sulfinato)-bridged complex, [Fe2III(μ-O)((2-NH2) p-CF3C6H3SO2)2](OTf)2 ( 5), was obtained. Parallel enzymatic studies on the CDO variant C93G were carried out with the abt substrate and show that reaction with O2 leads to disulfide formation, as opposed to S-oxygenation. The combined model and enzyme studies show that the thiol dioxygenases can operate via a 6-coord FeII center, in contrast to the accepted mechanism for nonheme iron dioxygenases, and that proper substrate chelation to Fe appears to be critical for S-oxygenation.

A Nonheme Sulfur-Ligated {FeNO}6 Complex and Comparison with Redox-Interconvertible {FeNO}7 and {FeNO}8 Analogues

A nonheme {FeNO}⁶ complex, [Fe(NO)(N3PyS)]²⁺ , was synthesized by reversible, one-electron oxidation of an {FeNO}⁷ analogue. This complex completes the first known series of sulfur-ligated {FeNO}^6-8 complexes. All three {FeNO}^6-8 complexes are readily interconverted by one-electron oxidation/reduction. A comparison of spectroscopic data (UV/Vis, NMR, IR, Mössbauer, X-ray absorption) provides a complete picture of the electronic and structural changes that occur upon {FeNO}⁶ -{FeNO}⁸ interconversion. Dissociation of NO from the new {FeNO}⁶ complex is shown to be controlled by solvent, temperature, and photolysis, which is rare for a sulfur-ligated {FeNO}⁶ species.

Sample treatment in Mössbauer spectroscopy for protein-related analyses: Nondestructive possibilities to look inside metal-containing biosystems

In this review, the unique possibilities are considered of the ⁵⁷Fe transmission (TMS) and ⁵⁷Co emission (EMS) variants of Mössbauer (nuclear γ-resonance) spectroscopy as nondestructive techniques with minimal sample preparation/treatment and a significant analytical potential, with a focus on the analysis of cation-binding sites in metalloproteins. The techniques are shown to provide unique structural and quantitative information on the coordination microenvironment, the chemical state and transformations of the Mössbauer nuclides in sophisticated metal-containing proteins, including those within complicated supramolecular structures, and in microbial cells or tissues. Recent representative examples of analyses of Fe-containing proteins by ⁵⁷Fe TMS are briefly discussed, along with the newly emerging data on using ⁵⁷Co EMS for probing the structural organisation of ⁵⁷Co-doped cation-binding sites in sophisticated biocomplexes including metalloenzymes. Finally, some rare or exotic applications of Mössbauer spectroscopy (including the synchrotron-based methodology) in protein-related studies are outlined.

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

The relative ease of Mössbauer spectroscopy and of density functional theory (DFT) calculations encourages the use of Mössbauer parameters as a validation method for calculations, and the use of calculations as a double check on crystallographic structures. A number of studies have proposed correlations between the computationally determined electron density at the iron nucleus and the observed isomer shift, but deviations from these correlations in low-valent iron β-diketiminate complexes encouraged us to determine a new correlation for these compounds. The use of B3LYP/def2-TZVP in the ORCA platform provides an excellent balance of accuracy and speed. We provide here not only this new correlation and a clear guide to its use but also a systematic analysis of the limitations of this approach. We also highlight the impact of crystallographic inaccuracies, DFT model truncation, and spin states, with intent to assist experimentalists to use Mössbauer spectroscopy and calculations together.

Reliable Estimation of Prediction Uncertainty for Physicochemical Property Models

One of the major challenges in computational science is to determine the uncertainty of a virtual measurement, that is the prediction of an observable based on calculations. As highly accurate first-principles calculations are in general unfeasible for most physical systems, one usually resorts to parameteric property models of observables, which require calibration by incorporating reference data. The resulting predictions and their uncertainties are sensitive to systematic errors such as inconsistent reference data, parametric model assumptions, or inadequate computational methods. Here, we discuss the calibration of property models in the light of bootstrapping, a sampling method that can be employed for identifying systematic errors and for reliable estimation of the prediction uncertainty. We apply bootstrapping to assess a linear property model linking the ⁵⁷Fe Mössbauer isomer shift to the contact electron density at the iron nucleus for a diverse set of 44 molecular iron compounds. The contact electron density is calculated with 12 density functionals across Jacob's ladder (PWLDA, BP86, BLYP, PW91, PBE, M06-L, TPSS, B3LYP, B3PW91, PBE0, M06, TPSSh). We provide systematic-error diagnostics and reliable, locally resolved uncertainties for isomer-shift predictions. Pure and hybrid density functionals yield average prediction uncertainties of 0.06-0.08 mm s^-1 and 0.04-0.05 mm s^-1, respectively, the latter being close to the average experimental uncertainty of 0.02 mm s^-1. Furthermore, we show that both model parameters and prediction uncertainty depend significantly on the composition and number of reference data points. Accordingly, we suggest that rankings of density functionals based on performance measures (e.g., the squared coefficient of correlation, r², or the root-mean-square error, RMSE) should not be inferred from a single data set. This study presents the first statistically rigorous calibration analysis for theoretical Mössbauer spectroscopy, which is of general applicability for physicochemical property models and not restricted to isomer-shift predictions. We provide the statistically meaningful reference data set MIS39 and a new calibration of the isomer shift based on the PBE0 functional.

Rational design of Fe catalysts for olefin aziridination through DFT-based mechanistic analysis

Experimental and DFT-based mechanistic studies are used to optimise Fe catalysts for aziridination.

Formation and structure of the ferryl [FeO] intermediate in the non-haem iron halogenase SyrB2: classical and QM/MM modelling agree

All O₂ activation roads for three substrates and three spin states in SyrB2 lead to the same [FeO] structure.

Comment on “Calibration of 57Fe Mössbauer constants by first principles” Phys. Chem. Chem. Phys., 2016, 18, 10201–10206

The proportionality constant between the iron-57 Mössbauer isomer shift and the electron probablility density at the nucleus is reevaluated from the correct experimental isomer shifts.