Ab Initio Calculations for Spin-Gaps of Non-Heme Iron Complexes

Force-Activated Spin-Crossover in Fe2+ and Co2+ Transition Metal Mechanophores

Transition metal mechanophores exhibiting force-activated spin-crossover are attractive design targets, yet large-scale discovery of them has not been pursued due in large part to the time-consuming nature of trial-and-error experiments. Instead, we leverage density functional theory (DFT) and external force explicitly included (EFEI) modeling to study a set of 395 feasible Fe2+ and Co2+ mechanophore candidates with tridentate ligands that we curate from the Cambridge Structural Database. Among nitrogen-coordinating low-spin complexes, we observe the prevalence of spin crossover at moderate force, and we identify 155 Fe2+ and Co2+ spin-crossover mechanophores and derive their threshold force for low-spin to high-spin transition (FSCO). The calculations reveal strong correlations of FSCO with spin-splitting energies and coordination bond lengths, facilitating rapid prediction of FSCO using force-free DFT calculations. Then, among all Fe2+ and Co2+ spin-crossover mechanophores, we further identity 11 mechanophores that combine labile spin-crossover and good mechanical robustness that are thus predicted to be the most versatile for force-probing applications. We discover two classes of mer-symmetric complexes comprising specific heteroaromatic rings within extended π-conjugation that give rise to Fe2+ mechanophores with these characteristics. We expect the set of spin-crossover mechanophores, the design principles, and the computational approach to be useful in guiding the high-throughput discovery of transition metal mechanophores with diverse functionalities and broad applications, including mechanically activated catalysis.

Unraveling the Geometry-Driven C═C Epoxidation and C-H Hydroxylation Reactivity of Tetra-Coordinated Nonheme Iron(IV)-Oxo Complexes

The electronic structure and reactivity of tetra-coordinated nonheme iron(IV)-oxo complexes have remained unexplored for years. The recent synthesis of a closed-shell iron(IV)-oxo complex [(quinisox)FeIV(O)]+ (1) has set up a platform to understand how such complexes compare with the celebrated open-shell iron-oxo chemistry. Herein, using density functional theory and ab initio calculations, we present an in-depth electronic structure investigation of the C═C epoxidation [oxygen atom transfer (OAT)] and C-H hydroxylation [hydrogen atom transfer (HAT)] reactivity of 1. Using a solvent-coordinated geometry of 1 (1') and other potential tetra-coordinated iron(IV)-oxo complexes bearing rigid ligands (2 and 3), we established the geometric origin of spin-state energetics and reactivity of 1. Complex 1 featuring a strong Fe-O bond exhibits OAT and HAT reactivity in its quintet state. The lowest quintet OAT pathway has a lower barrier by ∼4 kcal/mol than the quintet HAT pathway, corroborating the experimentally observed gas-phase OAT reactivity preference. A conventional HAT reactivity preference for 2 and a comparable OAT and HAT reactivity for 3 are observed. This further supports the geometry-driven reactivity preference for 1. Noncovalent interaction analyses reveal a pronounced π-π interaction between the substrate and ligand in the OAT transition state, rationalizing the origin of the observed reactivity preference for 1.

The spin-forbidden transition in iron(IV)-oxo catalysts relevant to two-state reactivity

Quintet oxoiron(IV) intermediates are often invoked in nonheme iron enzymes capable of performing selective oxidation, while most well-characterized synthetic model oxoiron(IV) complexes have a triplet ground state. These differing spin states lead to the proposal of a two-state reactivity model, where the complexes cross from the triplet to an excited quintet state. However, the energy of this quintet state has never been measured experimentally. Here, magnetic circular dichroism is used to assign the singlet and triplet excited states in a series of triplet oxoiron(IV) complexes. These transition energies are used to determine the energies of the quintet state via constrained fitting of 2p3d resonant inelastic x-ray scattering. This allowed for a direct correlation between the quintet energies and substrate C─H oxidation rates.

Correlating Structure with Spectroscopy in Ascorbate Peroxidase Compound II

Structural and spectroscopic investigations of compound II in ascorbate peroxidase (APX) have yielded conflicting conclusions regarding the protonation state of the crucial Fe(IV) intermediate. Neutron diffraction and crystallographic data support an iron(IV)-hydroxo formulation, whereas Mössbauer, X-ray absorption (XAS), and nuclear resonance vibrational spectroscopy (NRVS) studies appear consistent with an iron(IV)-oxo species. Here we examine APX with spectroscopy-oriented QM/MM calculations and extensive exploration of the conformational space for both possible formulations of compound II. We establish that irrespective of variations in the orientation of a vicinal arginine residue and potential reorganization of proximal water molecules and hydrogen bonding, the Fe-O distances for the oxo and hydroxo forms consistently fall within distinct, narrow, and nonoverlapping ranges. The accuracy of geometric parameters is validated by coupled-cluster calculations with the domain-based local pair natural orbital approach, DLPNO-CCSD(T). QM/MM calculations of spectroscopic properties are conducted for all structural variants, encompassing Mössbauer, optical, X-ray absorption, and X-ray emission spectroscopies and NRVS. All spectroscopic observations can be assigned uniquely to an Fe(IV)═O form. A terminal hydroxy group cannot be reconciled with the spectroscopic data. Under no conditions can the Fe(IV)═O distance be sufficiently elongated to approach the crystallographically reported Fe-O distance. The latter is consistent only with a hydroxo species, either Fe(IV) or Fe(III). Our findings strongly support the Fe(IV)═O formulation of APX-II and highlight unresolved discrepancies in the nature of samples used across different experimental studies.

Rebound or Cage Escape? The Role of the Rebound Barrier for the Reactivity of Non-Heme High-Valent FeIV =O Species

Owing to their high reactivity and selectivity, variations in the spin ground state and a range of possible pathways, high-valent FeIV =O species are popular models with potential bioinspired applications. An interesting example of a structure-reactivity pattern is the detailed study with five nonheme amine-pyridine pentadentate ligand FeIV =O species, including N4py: [(L1 )FeIV =O]2+ (1), bntpen: [(L2 )FeIV =O]2+ (2), py2 tacn: [(L3 )FeIV =O]2+ (3), and two isomeric bispidine derivatives: [(L4 )FeIV =O]2+ (4) and [(L5 )FeIV =O]2+ (5). In this set, the order of increasing reactivity in the hydroxylation of cyclohexane differs from that with cyclohexadiene as substrate. A comprehensive DFT, ab initio CASSCF/NEVPT2 and DLPNO-CCSD(T) study is presented to untangle the observed patterns. These are well reproduced when both activation barriers for the C-H abstraction and the OH rebound are taken into account. An MO, NBO and deformation energy analysis reveals the importance of π(pyr) → π*xz (FeIII -OH) electron donation for weakening the FeIII -OH bond and thus reducing the rebound barrier. This requires that pyridine rings are oriented perpendicularly to the FeIII -OH bond and this is a subtle but crucial point in ligand design for non-heme iron alkane hydroxylation.

Benchmarking Density Functional Theory Methods for Metalloenzyme Reactions: The Introduction of the MME55 Set

We present a new benchmark set of metalloenzyme model reaction energies and barrier heights that we call MME55. The set contains 10 different enzymes, representing eight transition metals, both open and closed shell systems, and system sizes of up to 116 atoms. We use four DLPNO-CCSD(T)-based approaches to calculate reference values against which we then benchmark the performance of a range of density functional approximations with and without dispersion corrections. Dispersion corrections improve the results across the board, and triple-ζ basis sets provide the best balance of efficiency and accuracy. Jacob's ladder is reproduced for the whole set based on averaged mean absolute (percent) deviations, with the double hybrids SOS0-PBE0-2-D3(BJ) and revDOD-PBEP86-D4 standing out as the most accurate methods for the MME55 set. The range-separated hybrids ωB97M-V and ωB97X-V also perform well here and can be recommended as a reliable compromise between accuracy and efficiency; they have already been shown to be robust across many other types of chemical problems, as well. Despite the popularity of B3LYP in computational enzymology, it is not a strong performer on our benchmark set, and we discourage its use for enzyme energetics.

Local Oxidation States in {FeNO}6-8 Porphyrins: Insights from DMRG/CASSCF-CASPT2 Calculations

A first DMRG/CASSCF-CASPT2 study of a series of paradigmatic {FeNO}6, {FeNO}7, and {FeNO}8 heme-nitrosyl complexes has led to substantial new insight as well as uncovered key shortcomings of the DFT approach. By virtue of its balanced treatment of static and dynamic correlation, the calculations have provided some of the most authoritative information available to date on the energetics of low- versus high-spin states of different classes of heme-nitrosyl complexes. Thus, the calculations indicate low doublet-quartet gaps of 1-4 kcal/mol for {FeNO}7 complexes and high singlet-triplet gaps of ≳20 kcal/mol for both {FeNO}6 and {FeNO}8 complexes. In contrast, DFT calculations yield widely divergent spin state gaps as a function of the exchange-correlation functional. DMRG-CASSCF calculations also help calibrate DFT spin densities for {FeNO}7 complexes, pointing to those obtained from classic pure functionals as the most accurate. The general picture appears to be that nearly all the spin density of Fe[P](NO) is localized on the Fe, while the axial ligand imidazole (ImH) in Fe[P](NO)(ImH) pushes a part of the spin density onto the NO moiety. An analysis of the DMRG-CASSCF wave function in terms of localized orbitals and of the resulting configuration state functions in terms of resonance forms with varying NO(π*) occupancies has allowed us to address the longstanding question of local oxidation states in heme-nitrosyl complexes. The analysis indicates NO(neutral) resonance forms [i.e., Fe(II)-NO0 and Fe(III)-NO0] as the major contributors to both {FeNO}6 and {FeNO}7 complexes. This finding is at variance with the common formulation of {FeNO}6 hemes as Fe(II)-NO+ species but is consonant with an Fe L-edge XAS analysis by Solomon and co-workers. For the {FeNO}8 complex {Fe[P](NO)}-, our analysis suggests a resonance hybrid description: Fe(I)-NO0 ↔ Fe(II)-NO-, in agreement with earlier DFT studies. Vibrational analyses of the compounds studied indicate an imperfect but fair correlation between the NO stretching frequency and NO(π*) occupancy, highlighting the usefulness of vibrational data as a preliminary indicator of the NO oxidation state.

The Nature of the Chemical Bonds of High-Valent Transition-Metal Oxo (M=O) and Peroxo (MOO) Compounds: A Historical Perspective of the Metal Oxyl-Radical Character by the Classical to Quantum Computations

This review article describes a historical perspective of elucidation of the nature of the chemical bonds of the high-valent transition metal oxo (M=O) and peroxo (M-O-O) compounds in chemistry and biology. The basic concepts and theoretical backgrounds of the broken-symmetry (BS) method are revisited to explain orbital symmetry conservation and orbital symmetry breaking for the theoretical characterization of four different mechanisms of chemical reactions. Beyond BS methods using the natural orbitals (UNO) of the BS solutions, such as UNO CI (CC), are also revisited for the elucidation of the scope and applicability of the BS methods. Several chemical indices have been derived as the conceptual bridges between the BS and beyond BS methods. The BS molecular orbital models have been employed to explain the metal oxyl-radical character of the M=O and M-O-O bonds, which respond to their radical reactivity. The isolobal and isospin analogy between carbonyl oxide R2C-O-O and metal peroxide LFe-O-O has been applied to understand and explain the chameleonic chemical reactivity of these compounds. The isolobal and isospin analogy among Fe=O, O=O, and O have also provided the triplet atomic oxygen (3O) model for non-heme Fe(IV)=O species with strong radical reactivity. The chameleonic reactivity of the compounds I (Cpd I) and II (Cpd II) is also explained by this analogy. The early proposals obtained by these theoretical models have been examined based on recent computational results by hybrid DFT (UHDFT), DLPNO CCSD(T0), CASPT2, and UNO CI (CC) methods and quantum computing (QC).

Metastable Charged Dimers in Organometallic Species: A Look into Hydrogen Bonding between Metallocene Derivatives

Multiply charged complexes bound by noncovalent interactions have been previously described in the literature, although they were mostly focused on organic and main group inorganic systems. In this work, we show that similar complexes can also be found for organometallic systems containing transition metals and deepen in the reasons behind the existence of these species. We have studied the structures, binding energies, and dissociation profiles in the gas phase of a series of charged hydrogen-bonded dimers of metallocene (Ru, Co, Rh, and Mn) derivatives isoelectronic with the ferrocene dimer. Our results indicate that the carboxylic acid-containing dimers are more strongly bonded and present larger barriers to dissociation than the amide ones and that the cationic complexes tend to be more stable than the anionic ones. Additionally, we describe for the first time the symmetric proton transfer that can occur while in the metastable phase. Finally, we use a density-based energy decomposition analysis to shine light on the nature of the interaction between the dimers.

Unraveling the Spin-State Energetics of FeN4 Complexes with Ab Initio Methods

A systematic analysis was conducted to explore the spin-state energetics of a series of 19 FeN4 complexes. The performance of a large number of multireference methods was assessed, highlighting the significant challenges associated with accurately describing the spin-state energetics of FeN4 complexes. Most multireference methods were found to be susceptible to errors originating from the reference CASSCF wavefunction, leading to an overstabilization of high-spin states. Nonetheless, a few multireference methods, namely, CASPT2/CC, DSRG-MRPT3, and LDSRG(2), demonstrated promising performance compared to the benchmark CCSD(T) method. Furthermore, our study revealed that FeN4 complexes having a quintet ground state are exceedingly rare. Accordingly, only one specific model (Fe(L2)) and one synthesized complex (Fe(OTBP)) have the quintet ground state among the studied complexes. This scarcity of quintet FeN4 complexes highlights the unique nature of these systems and raises intriguing questions regarding the factors influencing spin states, such as the size of the macrocycle cavity, the introduction of substituents, or the induction of out-of-plane deformation.

MP2-Based Correction Scheme to Approach the Limit of a Complete Pair Natural Orbitals Space in DLPNO-CCSD(T) Calculations

The domain-based local pair natural orbital (PNO) coupled-cluster DLPNO-CCSD(T) method has been proven to provide accurate single-point energies at a fraction of the cost of canonical CCSD(T) calculations. However, the desired "chemical accuracy" can only be obtained with a large PNO space and extended basis set. We present a simple yet accurate and efficient correction scheme based on a perturbative approach. Here, in addition to DLPNO-CCSD(T) energy, one calculates DLPNO-MP2 correlation energy with the same settings as in the preceding coupled-cluster calculation. In the next step, the canonical MP2 correlation energy is obtained in the same orbital basis. This can be efficiently performed for essentially all molecule sizes accessible with the DLPNO-CCSD(T) method. By taking the difference between the canonical MP2 and DLPNO-MP2 energies, we obtain a correction term that can be added to the DLPNO-CCSD(T) correlation energy. This way, one can obtain the total correlation energy close to the limit of the complete PNO space (cPNO). The presented approach allows us to significantly increase the accuracy of the DLPNO-CCSD(T) method for both closed- and open-shell systems. The latter are known to be especially challenging for locally correlated methods. Unlike the previously developed PNO extrapolation procedure by Altun, Neese, and Bistoni ( J. Chem. Theory Comput. 2020, 16, 6142-6149), this strategy allows us to get the DLPNO-CCSD(T) correlation energy at the cPNO limit in a cost-efficient way, resulting in a minimal overall increase in calculation time as compared to the uncorrected method.

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within "chemical accuracy" (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob's ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 - N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Advances in Parallel Heat Bath Configuration Interaction

Heat-bath configuration interaction (HCI) is a deterministic method that approaches the full CI limit at greatly reduced computational cost. In this work, computational improvements to the HCI algorithm are introduced targeting speed, parallel efficiency, and memory requirements. The new implementation introduces a hash function to distribute determinants and takes advantage of MPI and OpenMP for parallelism allowing for a (22e,168o) active space to be studied, which explicitly includes 2.39 × 10⁷ variational determinants and 8.95 × 10¹⁰ perturbative determinants. Benchmarks show up to 86% parallel efficiency of the perturbative step on 32 nodes (4096 cores) and a total efficiency of 74%. The new HCI implementation is benchmarked for accuracy against prior results and applied to study the triplet-quintet gap in the challenging [FeO(NH₃)₅]²⁺ complex.

Comparative oxidative ability of mononuclear and dinuclear high-valent iron-oxo species towards the activation of methane: does the axial/bridge atom modulate the reactivity?

Over the years, mononuclear Fe^IVO species have been extensively studied, but the presence of dinuclear Fe^IVO species in soluble methane monooxygenase (sMMO) has inspired the development of biomimic models that could activate inert substrates such as methane. There are some successful attempts; particularly the [(Por)(m-CBA) Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- species has been reported to activate methane and yield decent catalytic turnover numbers and therefore regarded as the closest to the sMMO enzyme functional model, as no mononuclear Fe^IVO analogues could achieve this feat. In this work, we have studied a series of mono and dinuclear models using DFT and ab initio DLPNO-CCSD(T) calculations to probe the importance of nuclearity in enhancing the reactivity. We have probed the catalytic activities of four complexes: [(HO)Fe^IV(O)(Por)]^- (1), [(HO)Fe^IV(O)(Por˙⁺)] (2), μ-oxo dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-O)Fe^IV(O) (Por˙⁺)]^- (3) and N-bridged dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- (4) towards the activation of methane. Additionally, calculations were performed on the mononuclear models [(X)Fe^IV(O)(Por˙⁺)]ⁿ {X = N 4a (n = -2), NH 4b (n = -1) and NH₂4c (n = 0)} to understand the role of nuclearity in the reactivity. DFT calculations performed on species 1-4 suggest an interesting variation among them, with species 1-3 possessing an intermediate spin (S = 1) as a ground state and species 4 possessing a high-spin (S = 2) as a ground state. Furthermore, the two Fe^IV centres in species 3 and 4 are antiferromagnetically coupled, yielding a singlet state with a distinct difference in their electronic structure. On the other hand, species 2 exhibits a ferromagnetic coupling between the Fe^IV and the Por˙⁺ moiety. Our calculations suggest that the higher barriers for the C-H bond activation of methane and the rebound step for species 1 and 3 are very high in energy, rendering them unreactive towards methane, while species 2 and 4 have lower barriers, suggesting their reactivity towards methane. Studies on the system reveal that model 4a has multiple FeN bonds facilitating greater reactivity, whereas the other two models have longer Fe-N bonds and less radical character with steeper barriers. Strong electronic cooperativity is found to be facilitated by the bridging nitride atom, and this cooperativity is suppressed by substituents such as oxygen, rendering them inactive. Thus, our study unravels that apart from enhancing the nuclearity, bridging atoms that facilitate strong cooperation between the metals are required to activate very inert substrates such as methane, and our results are broadly in agreement with earlier experimental findings.

Reconciling Local Coupled Cluster with Multireference Approaches for Transition Metal Spin-State Energetics

Spin-state energetics of transition metal complexes remain one of the most challenging targets for electronic structure methods. Among single-reference wave function approaches, local correlation approximations to coupled cluster theory, most notably the domain-based local pair natural orbital (DLPNO) approach, hold the promise of bringing the accuracy of coupled cluster theory with single, double, and perturbative triple excitations, CCSD(T), to molecular systems of realistic size with acceptable computational cost. However, recent studies on spin-state energetics of iron-containing systems raised doubts about the ability of the DLPNO approach to adequately and systematically approximate energetics obtained by the reference-quality complete active space second-order perturbation theory with coupled-cluster semicore correlation, CASPT2/CC. Here, we revisit this problem using a diverse set of iron complexes and examine several aspects of the application of the DLPNO approach. We show that DLPNO-CCSD(T) can accurately reproduce both CASPT2/CC and canonical CCSD(T) results if two basic principles are followed. These include the consistent use of the improved iterative (T₁) versus the semicanonical perturbative triple corrections and, most importantly, a simple two-point extrapolation to the PNO space limit. The latter practically eliminates errors arising from the default truncation of electron-pair correlation spaces and should be viewed as standard practice in applications of the method to transition metal spin-state energetics. Our results show that reference-quality results can be readily achieved with DLPNO-CCSD(T) if these principles are followed. This is important also in view of the applicability of the method to larger single-reference systems and multinuclear clusters, whose treatment of dynamic correlation would be challenging for multireference-based approaches.

Comprehensive Theoretical View of the [Cu2O2] Side-on-Peroxo-/Bis-μ-Oxo Equilibria

Coupled binuclear copper (CBC) sites are employed by many metalloenzymes to catalyze a broad set of biochemical transformations. Typically, the CBC catalytic sites are activated by the O2 molecule to form various [Cu2O2] reactive species. This has also inspired synthesis and development of various biomimetic inorganic complexes featuring the CBC core. From theoretical perspective, the [Cu2O2] reactivity often hinges on the side-on-peroxo-dicopper(II) (P) vs. bis-μ-oxo-dicopper(III) (O) isomerism - an equilibrium that has become almost iconic in theoretical bioinorganic chemistry. Herein, we present a comprehensive calibration and evaluation of the performance of various composite computational protocols available in contemporary computational chemistry, involving coupled-cluster and multi-reference (relativistic) wave function methods, popular density functionals and solvation models. Starting with the well-studied reference [Cu2O2(NH3)6]2+ system, we compared the performance of electronic structure methods and discussed the relativistic effects. This allowed us to select several 'calibrated' DFT functionals that can be conveniently employed to study ten experimentally well-characterized [Cu2O2] inorganic systems. We mostly correctly predicted the lowest-energy structures (P vs. O) of the studied systems. In addition, we present calibration of the used electronic structure methods for prediction of the spectroscopic features of the [Cu2O2] core.

The many-body electronic interactions of Fe(II)–porphyrin

Fe(II)–porphyrin complexes exhibit a diverse range of electronic interactions between the metal and macrocycle. Herein, the incremental full configuration interaction method is applied to the entire space of valence orbitals of a Fe(II)–porphyrin model using a modest basis set. A novel visualization framework is proposed to analyze individual many-body contributions to the correlation energy, providing detailed maps of this complex’s highly correlated electronic structure. This technique is used to parse the numerous interactions of two low-lying triplet states (3A2g and 3Eg) and to show that strong metal d–d and macrocycle π–π orbital interactions preferentially stabilize the 3A2g state. d–π interactions, on the other hand, preferentially stabilize the 3Eg state and primarily appear when correlating six electrons at a time. Ultimately, the Fe(II)–porphyrin model’s full set of 88 valence electrons are correlated in 275 orbitals, showing the interactions up to the 4-body level, which covers the great majority of correlations in this system.

Hydrosilylation of Aldehydes and Ketones Catalyzed by a 2-Iminopyrrolyl Alkyl-Manganese(II) Complex

A well-defined and very active single-component manganese(II) catalyst system for the hydrosilylation of aldehydes and ketones is presented. First, the reaction of 5-(2,4,6-iPr₃C₆H₂)-2-[N-(2,6-iPr₂C₆H₃)formimino]pyrrolyl potassium (KL) and [MnCl₂(Py)₂] afforded the binuclear 2-iminopyrrolyl manganese(II) pyridine chloride complex [Mn₂{κ²N,N'-5-(2,4,6-iPr₃C₆H₂)-NC₄H₂-2-C(H)═N(2,6-iPr₂C₆H₃)}₂(Py)₂(μ-Cl)₂] 1. Subsequently, the alkylation reaction of complex 1 with LiCH₂SiMe₃ afforded the respective (trimethylsilyl)methyl-Mn(II) complex [Mn{κ²N,N'-5-(2,4,6-iPr₃C₆H₂)-NC₄H₂-2-C(H)═N(2,6-iPr₂C₆H₃)}(Py)CH₂SiMe₃] 2 in a good yield. Complexes 1 and 2 were characterized by elemental analysis, ¹H NMR spectroscopy, Evans' method, FTIR spectroscopy, and single-crystal X-ray diffraction. While the crystal structure of complex 1 has been identified as a binuclear entity, in which the Mn(II) centers present pentacoordinate coordination spheres, that of complex 2 corresponds to a monomer with a distorted tetrahedral coordination geometry. Complex 2 proved to be a very active precatalyst for the atom-economic hydrosilylation of several aldehydes and ketones under very mild conditions, with a maximum turnover frequency of 95 min^-1, via a silyl-Mn(II) mechanistic route, as asserted by a combination of experimental and theoretical efforts, the respective silanes were cleanly converted to the respective alcoholic products in high yields.

Redox Isomerism in the S3 State of the Oxygen-Evolving Complex Resolved by Coupled Cluster Theory

The electronic and geometric structures of the water-oxidizing complex of photosystem II in the steps of the catalytic cycle that precede dioxygen evolution remain hotly debated. Recent structural and spectroscopic investigations support contradictory redox formulations for the active-site Mn4 CaOx cofactor in the final metastable S3 state. These range from the widely accepted MnIV 4 oxo-hydroxo model, which presumes that O-O bond formation occurs in the ultimate transient intermediate (S4 ) of the catalytic cycle, to a MnIII 2 MnIV 2 peroxo model representative of the contrasting "early-onset" O-O bond formation hypothesis. Density functional theory energetics of suggested S3 redox isomers are inconclusive because of extreme functional dependence. Here, we use the power of the domain-based local pair natural orbital approach to coupled cluster theory, DLPNO-CCSD(T), to present the first correlated wave function theory calculations of relative stabilities for distinct redox-isomeric forms of the S3 state. Our results enabled us to evaluate conflicting models for the S3 state of the oxygen-evolving complex (OEC) and to quantify the accuracy of lower-level theoretical approaches. Our assessment of the relevance of distinct redox-isomeric forms for the mechanism of biological water oxidation strongly disfavors the scenario of early-onset O-O formation advanced by literal interpretations of certain crystallographic models. This work serves as a case study in the application of modern coupled cluster implementations to redox isomerism problems in oligonuclear transition metal systems.

Deciphering the origin of million-fold reactivity observed for the open core diiron [HO-FeIII-O-FeIV[double bond, length as m-dash]O]2+ species towards C-H bond activation: role of spin-states, spin-coupling, and spin-cooperation

High-valent metal-oxo species have been characterised as key intermediates in both heme and non-heme enzymes that are found to perform efficient aliphatic hydroxylation, epoxidation, halogenation, and dehydrogenation reactions. Several biomimetic model complexes have been synthesised over the years to mimic both the structure and function of metalloenzymes. The diamond-core [Fe2(μ-O)2] is one of the celebrated models in this context as this has been proposed as the catalytically active species in soluble methane monooxygenase enzymes (sMMO), which perform the challenging chemical conversion of methane to methanol at ease. In this context, a report of open core [HO(L)FeIII-O-FeIV(O)(L)]2+ (1) gains attention as this activates C-H bonds a million-fold faster compared to the diamond-core structure and has the dual catalytic ability to perform hydroxylation as well as desaturation with organic substrates. In this study, we have employed density functional methods to probe the origin of the very high reactivity observed for this complex and also to shed light on how this complex performs efficient hydroxylation and desaturation of alkanes. By modelling fifteen possible spin-states for 1 that could potentially participate in the reaction mechanism, our calculations reveal a doublet ground state for 1 arising from antiferromagnetic coupling between the quartet FeIV centre and the sextet FeIII centre, which regulates the reactivity of this species. The unusual stabilisation of the high-spin ground state for FeIV[double bond, length as m-dash]O is due to the strong overlap of with the orbital, reducing the antibonding interactions via spin-cooperation. The electronic structure features computed for 1 are consistent with experiments offering confidence in the methodology chosen. Further, we have probed various mechanistic pathways for the C-H bond activation as well as -OH rebound/desaturation of alkanes. An extremely small barrier height computed for the first hydrogen atom abstraction by the terminal FeIV[double bond, length as m-dash]O unit was found to be responsible for the million-fold activation observed in the experiments. The barrier height computed for -OH rebound by the FeIII-OH unit is also smaller suggesting a facile hydroxylation of organic substrates by 1. A strong spin-cooperation between the two iron centres also reduces the barrier for second hydrogen atom abstraction, thus making the desaturation pathway competitive. Both the spin-state as well as spin-coupling between the two metal centres play a crucial role in dictating the reactivity for species 1. By exploring various mechanistic pathways, our study unveils the fact that the bridged μ-oxo group is a poor electrophile for both C-H activation as well for -OH rebound. As more and more evidence is gathered in recent years for the open core geometry of sMMO enzymes, the idea of enhancing the reactivity via an open-core motif has far-reaching consequences.

Seeing Is Believing: Experimental Spin States from Machine Learning Model Structure Predictions

Determination of ground-state spins of open-shell transition-metal complexes is critical to understanding catalytic and materials properties but also challenging with approximate electronic structure methods. As an alternative approach, we demonstrate how structure alone can be used to guide assignment of ground-state spin from experimentally determined crystal structures of transition-metal complexes. We first identify the limits of distance-based heuristics from distributions of metal–ligand bond lengths of over 2000 unique mononuclear Fe(II)/Fe(III) transition-metal complexes. To overcome these limits, we employ artificial neural networks (ANNs) to predict spin-state-dependent metal–ligand bond lengths and classify experimental ground-state spins based on agreement of experimental structures with the ANN predictions. Although the ANN is trained on hybrid density functional theory data, we exploit the method-insensitivity of geometric properties to enable assignment of ground states for the majority (ca. 80–90%) of structures. We demonstrate the utility of the ANN by data-mining the literature for spin-crossover (SCO) complexes, which have experimentally observed temperature-dependent geometric structure changes, by correctly assigning almost all (>95%) spin states in the 46 Fe(II) SCO complex set. This approach represents a promising complement to more conventional energy-based spin-state assignment from electronic structure theory at the low cost of a machine learning model.

Unravelling the mechanism of cobalt-catalysed remote C-H nitration of 8-aminoquinolinamides and expansion of substrate scope towards 1-naphthylpicolinamide

Previously, an unexpected Co-catalysed remote C-H nitration of 8-aminoquinolinamide compounds was developed. This report provided a novel reactivity for Co which was assumed to proceed through the mechanistic pathway already known for analogous Cu-catalysed remote couplings of the same substrates. In order to shed light into this intriguing, and previously unobserved reactivity for Co, a thorough computational study has now been performed, which has allowed for a full understanding of the operative mechanism. This study demonstrates that the Co-catalysed remote coupling does not occur through the previously proposed Single Electron Transfer (SET) mechanism, but actually operates through a high-spin induced remote radical coupling mechanism, through a key intermediate with significant proportion of spin density at the 5- and 7-positions of the aminoquinoline ring. Additionally, new experimental data provides expansion of the synthetic utility of the original nitration procedure towards 1-naphthylpicolinamide which unexpectedly appears to operate via a subtly different mechanism despite having a similar chelate environment.

Large-scale comparison of 3d and 4d transition metal complexes illuminates the reduced effect of exchange on second-row spin-state energetics

The origin of distinct 3d vs. 4d transition metal complex sensitivity to exchange is explored over a large data set.

Accurate computed spin-state energetics for Co(iii) complexes: implications for modelling homogeneous catalysis

DLPNO-CCSD(T) calculations provide accurate spin state energetics for a range of Co(iii) complexes and so represent a promising approach to modelling homogeneous catalysis based on Co(iii) species.

Appraising spin-state energetics in transition metal complexes using double-hybrid models: accountability of SOS0-PBESCAN0-2(a) as a promising paradigm

Through a comprehensive survey, reliable double-hybrid models have been validated and proposed for spin-state energetics in transition metal complexes.

Spin state and reactivity of iron(iv)oxido complexes with tetradentate bispidine ligands

The iron(iv)oxido complex [(bispidine)Fe^IVO(Cl)]⁺is shown by experiment and high-level DLPNO-CCSD(T) quantum-chemical calculations to be an extremely short-lived and very reactive intermediate-spin (S= 1) species.

Energetics of non-heme iron reactivity: can ab initio calculations provide the right answer?

We use a variety of computational methods to characterize and compare the hydrogen atom transfer (HAT) and epoxidation reaction pathways for oxidation of cyclohexene by an iron(iv)-oxo complex.