Spinning around in Transition-Metal Chemistry

Role of non-redox innocent ligand units in the oxidation of alcohols with H2O2 catalyzed by μ-oxido-diiron(III) bis-phenolato polypyridyl complexes

Redox non-innocent ligands hold the potential to expand the redox chemistry and activity of transition metal catalysts. The impact of the additional redox chemistry of phenol ligands in oxidation catalysis is explored here in the complex μ-oxido-diiron(III) polypyridyl (1) [(L)Fe(III)(μ-O)Fe(III)(L)](ClO4)2 (where HL is 2-(((di(pyridin-2-yl)methyl) (pyridin-2-ylmethyl) amino)methyl)phenol) and its tert-butyl substituted analog 2, in which each of the Fe(III) centers is coordinated to a phenolato moiety of the ligand. Complex 1 was shown earlier to catalyse the oxidation of benzyl alcohols to aldehydes with H2O2. In particular acid was found to accelerate the reactions by removal of a lag period before catalysis initiated. Here, we use reaction monitoring with resonance Raman, UV/vis absorption and EPR spectroscopy to show that under catalytic conditions, i.e. with excess H2O2, rapid (< 5 s) loss of the phenolato moiety occurs, resulting in the formation of an N4 ligated Fe(III) complex. This N4 coordinated complex forms a Fe(III)-OOH species, which is responsible for alcohol oxidation and over time a relatively stable oxido-bridged dinuclear Fe(III) complex forms as a resting state in the catalytic system. The main role of acid in the catalysis is shown to be to facilitate the initial coordination of H2O2 by driving the formation of mononuclear complexes from 1 and 2. The data show that although the phenolato moiety imparts interesting redox properties on complex 1, it does not contribute directly to the oxidation catalysis observed with H2O2.

Leveraging natural language processing to curate the tmCAT, tmPHOTO, tmBIO, and tmSCO datasets of functional transition metal complexes

The breadth of transition metal chemical space covered by databases such as the Cambridge Structural Database and the derived computational database tmQM is not conducive to application-specific modeling and the development of structure-property relationships. Here, we employ both supervised and unsupervised natural language processing (NLP) techniques to link experimentally synthesized compounds in the tmQM database to their respective applications. Leveraging NLP models, we curate four distinct datasets: tmCAT for catalysis, tmPHOTO for photophysical activity, tmBIO for biological relevance, and tmSCO for magnetism. Analyzing the chemical substructures within each dataset reveals common chemical motifs in each of the designated applications. We then use these common chemical structures to augment our initial datasets for each application, yielding a total of 21 631 compounds in tmCAT, 4599 in tmPHOTO, 2782 in tmBIO, and 983 in tmSCO. These datasets are expected to accelerate the more targeted computational screening and development of refined structure-property relationships with machine learning.

Automated prediction of ground state spin for transition metal complexes

Exploiting crystallographic data repositories for large-scale quantum chemical computations requires the rapid and accurate extraction of the molecular structure, charge and spin from the crystallographic information file. Here, we develop a general approach to assign the ground state spin of transition metal complexes, in complement to our previous efforts on determining metal oxidation states and bond order within the cell2mol software. Starting from a database of 31k transition metal complexes extracted from the Cambridge Structural Database with cell2mol, we construct the TM-GSspin dataset, which contains 2063 mononuclear first row transition metal complexes and their computed ground state spins. TM-GSspin is highly diverse in terms of metals, metal oxidation states, coordination geometries, and coordination sphere compositions. Based on TM-GSspin, we identify correlations between structural and electronic features of the complexes and their ground state spins to develop a rule-based spin state assignment model. Leveraging this knowledge, we construct interpretable descriptors and build a statistical model achieving 98% cross-validated accuracy in predicting the ground state spin across the board. Our approach provides a practical way to determine the ground state spin of transition metal complexes directly from crystal structures without additional computations, thus enabling the automated use of crystallographic data for large-scale computations involving transition metal complexes.

Benchmarking Periodic Density Functional Theory Calculations for Spin-State Energies in Spin-Crossover Systems

Spin energetics is one of the biggest challenges associated with energy calculations for electronic structure methods. The energy differences of the spin states in spin-crossover compounds are very small, making them one of the most difficult systems to calculate. Few methods provide accurate results for calculating these energy differences. In addition, studies have usually focused on calculating energetics of single molecules, while spin-crossover properties are usually experimentally studied in the solid phase. In this paper, we have used periodic boundary conditions employing methods based on density functional theory to calculate the high- and low-spin energy differences for a test case of 20 extended systems. Compounds with different metals and ligands have been selected, and the results indicate that a semiquantitative description of the energy differences can be obtained with the combination of geometry optimization using the PBE functional including many-body dispersion approach and the use of meta-GGA functionals, such as r2SCAN but especially KTBM24, for the energy calculation. Other hybrid functionals, such as TPSSh, give generally good results, but the calculation of the exact exchange with periodic boundary conditions involves a huge increase in computer time and computational resources. It makes the proposed nonhybrid functional approach (KTBM24//PBE+MB) a great advantage for the study of periodic systems.

Predicting spin states of iron porphyrins with DFT methods including crystal packing effects and thermodynamic corrections

Accurate computational treatment of spin states for transition metal complexes, exemplified by iron porphyrins, lies at the heart of quantum bioinorganic chemistry, but at the same time represents a great challenge for approximate density functional theory (DFT) methods, which are predominantly used. Here, the accuracy of DFT methods for spin-state splittings in iron porphyrin is assessed by probing the ability to correctly predict the ground states for six FeIII or FeII complexes experimentally characterized in solid state. For each case, molecular and periodic DFT calculations are employed to quantify the effect of porphyrin side substituents and the crystal packing effect (CPE) on the spin-state splitting. It is proposed to partition the total CPE into additive components, the direct and structural one, the importance of which is shown to significantly vary from case to case. By knowing the substituent effect, the CPE, and the Gibbs free energy thermodynamic correction from calculations, one can employ the experimental ground-state information in order to derive a quantitative constraint on the electronic energy difference for a simplified (porphin) model of the experimentally characterized metalloporphyrin. The constraints derived in such a way-in the form of single or double inequalities-are used to assess the accuracy of dispersion-corrected DFT methods for 6 spin-state splittings of [FeIII(P)(2-MeIm)2]+, [FeIII(P)(2-MeIm)]+, [FeII(P)(THF)2] and [FeII(P)] models (where P is porphin, 2-MeIm is 2-methylimidazole, THF is tetrahydrofuran). These data constitute the new benchmark set of spin states for crystalline iron porphyrins (SSCIP6). The highest accuracy is obtained in the case of double-hybrid functionals (B2PLYP-D3, DSD-PBEB95-D3), whereas hybrid functionals, especially those with reduced admixture of the exact exchange (B3LYP*-D3, TPSSh-D3), are found to considerably overstabilize the intermediate spin state, leading to incorrect ground-state prediction in FeIII porphyrins. The present approach, which can be generalized to other transition metal complexes, is not only useful in method benchmarking, but also sheds light on the interpretations of experimental data for metalloporphyrins, which are important models to understand the electronic properties of heme proteins.

Chelate Complexes of 3d Transition Metal Ions─A Challenge for Electronic-Structure Methods?

Different electronic-structure methods were assessed for their ability to predict two important properties of the industrially relevant chelating agent nitrilotriacetic acid (NTA): its selectivity with respect to six different first-row transition metal ions and the spin-state energetics of its complex with Fe(III). The investigated methods encompassed density functional theory (DFT), the random phase approximation (RPA), coupled cluster (CC) theory, and the auxiliary-field quantum Monte Carlo (AFQMC) method, as well as the complete active space self-consistent field (CASSCF) method and the respective on-top methods: second-order N-electron valence state perturbation theory (NEVPT2) and multiconfiguration pair-density functional theory (MC-PDFT). Different strategies for selecting active spaces were explored, and the density matrix renormalization group (DMRG) approach was used to solve the largest active spaces. Despite somewhat ambiguous multi-reference diagnostics, most methods gave relatively good agreement with experimental data for the chemical reactions connected to the selectivity, which only involved transition-metal complexes in their high-spin state. CC methods yielded the highest accuracy followed by range-separated DFT and AFQMC. We discussed in detail that even higher accuracies can be obtained with NEVPT2, under the prerequisite that consistent active spaces along the entire chemical reaction can be selected, which was not the case for reactions involving Fe(III). A bigger challenge for electronic-structure methods was the prediction of the spin-state energetics, which additionally involved lower spin states that exhibited larger multi-reference diagnostics. Conceptually different, typically accurate methods ranging from CC theory via DMRG-NEVPT2 in combination with large active spaces to AFQMC agreed well that the high-spin state is energetically significantly favored over the other spin states. This was in contrast to most DFT functionals and RPA which yielded a smaller stabilization and some common DFT functionals and MC-PDFT even predicting the low-spin state to be energetically most favorable.

DFT Struggles to Predict the Energy Landscape for Iron Pyridine Diimine-Catalyzed [2 + 2] Cycloaddition of Alkenes: Insights into the Problem and Alternative Solutions

In recent years, noninnocent pyridine diimine (PDI) complexes featuring first-row transition metals have emerged as prominent catalysts, demonstrating efficacy in a diverse range of vital organometallic transformations. However, the inherent complexity of the fundamental reactivity paradigm in these systems arises from the presence of a noninnocent ligand and the multispin feasibility of 3d metals. While density functional theory (DFT) has been widely used to unravel mechanistic insights, its limitations as a single-reference method can potentially misrepresent spin-state energetics, compromising our understanding of these intricate systems. In this study, we employ extensive high-level ab initio state averaged-complete active space self-consistent field/N-electron valence state perturbation theory (SA-CASSCF/NEVPT2) calculations in combination with DFT to investigate an iron-PDI-catalyzed [2 + 2] cycloaddition reaction of alkenes. The transformation proceeds through two major steps: oxidative cyclization and reductive elimination. Contrary to the predictions of DFT calculations, which suggest two-state reactivity in the reaction and identify reductive elimination as the turnover-limiting step, SA-CASSCF/NEVPT2-corrected results unequivocally establish a single-state reactivity scenario with oxidative cyclization as the turnover-limiting step. SA-CASSCF/NEVPT2-based insights into electronic ground states and electron distribution elucidate the intriguing interactions between the PDI ligand and the iron center, revealing the highly multiconfigurational nature of these species and providing a precise depiction of metal-ligand cooperativity throughout the transformation. A comparative assessment of several widely recognized DFT functionals against SA-CASSCF/NEVPT2-corrected data indicates that single-point energy calculations using the modern density functional MN15 on TPSSh geometries offer the most reliable density functional methodology, in scenarios where SA-CASSCF/NEVPT2 computational cost is a consideration.

Questing for homoleptic mononuclear manganese complexes with monodentate O-donor ligands

Compounds containing Mn-O bonds are of utmost importance in biological systems and catalytic processes. Nevertheless, mononuclear manganese complexes containing all O-donor ligands are still rare. Taking advantage of the low tendency of the pentafluoroorthotellurate ligand (teflate, OTeF5) to bridge metal centers, we have synthesized two homoleptic manganese complexes with monomeric structures and an all O-donor coordination sphere. The tetrahedrally distorted MnII anion, [Mn(OTeF5)4]2-, can be described as a high spin d5 complex (S = 5/2), as found experimentally (magnetic susceptibility measurements and EPR spectroscopy) and using theoretical calculations (DFT and CASSCF/NEVPT2). The high spin d4 electronic configuration (S = 2) of the MnIII anion, [Mn(OTeF5)5]2-, was also determined experimentally and theoretically, and a square pyramidal geometry was found to be the most stable one for this complex. Finally, the bonding situation in both complexes was investigated by means of the Interacting Quantum Atoms (IQA) methodology and compared to that of hypothetical mononuclear fluoromanganates. Within each pair of [MnXn]2- (n = 4, 5) species (X = OTeF5, F), the Mn-X interaction is found to be comparable, therefore proving that the similar electronic properties of the teflate and the fluoride are also responsible for the stabilization of these unique species.

An Update in Computational Methods for Environmental Monitoring: Theoretical Evaluation of the Molecular and Electronic Structures of Natural Pigment-Metal Complexes

Metals are beneficial to life, but the presence of these elements in excessive amounts can harm both organisms and the environment; therefore, detecting the presence of metals is essential. Currently, metal detection methods employ powerful instrumental techniques that require a lot of time and money. Hence, the development of efficient and effective metal indicators is essential. Several synthetic metal detectors have been made, but due to their risk of harm, the use of natural pigments is considered a potential alternative. Experiments are needed for their development, but they are expensive and time-consuming. This review explores various computational methods and approaches that can be used to investigate metal-pigment interactions because choosing the right methods and approaches will affect the reliability of the results. The results show that quantum mechanical methods (ab initio, density functional theory, and semiempirical approaches) and molecular dynamics simulations have been used. Among the available methods, the density functional theory approach with the B3LYP functional and the LANL2DZ ECP and basis set is the most promising combination due to its good accuracy and cost-effectiveness. Various experimental studies were also in good agreement with the results of computational methods. However, deeper analysis still needs to be carried out to find the best combination of functions and basis sets.

Theoretical Photoelectron Spectroscopy of Quadruple-Bonded Dimolybdenum(II,II) and Ditungsten(II,II) Paddlewheel Complexes: Performance of Common Density Functional Theory Methods

We have revisited the gas-phase photoelectron spectra of quadruple-bonded dimolybdenum(II,II) and ditungsten(II,II) paddlewheel complexes with modern density functional theory methods and obtained valuable calibration of four well-known exchange-correlation functionals, namely, BP86, OLYP, B3LYP*, and B3LYP. All four functionals were found to perform comparably, with discrepancies between calculated and experimental ionization potentials ranging from <0.1 to ∼0.5 eV, with the lowest errors observed for the classic pure functional BP86. All four functionals were found to reproduce differences in ionization potentials (IPs) between analogous Mo2 and W2 complexes, as well as large, experimentally observed ligand field effects on the IPs, with near-quantitative accuracy. The calculations help us interpret a number of differences between analogous Mo2 and W2 complexes through the lens of relativistic effects. Thus, relativity results in not only significantly lower IPs for the W2 complexes but also smaller HOMO-LUMO gaps and different triplet states relative to their Mo2 counterparts.

Peculiar Differences between Two Copper Complexes Containing Similar Redox-Active Ligands: Density Functional and Multiconfigurational Calculations

Transition metal complexes featuring redox-active ligands often exhibit multiple redox states, influenced by the interplay between the metal center and the ligand. This study delves into the electronic structures of two mononuclear complexes of copper with two similar redox-active urea azine ligands. The ligands differ by the replacement of an NCH3 moiety by an S atom in the ligand backbone. Experimental analysis yields pronounced electronic structural disparities between these complexes, observable in both the solution and solid phases. Conventional quantum chemical methods, such as density functional theory using different functionals (B3LYP, TPSSh, and CAM-B3LYP), remain inadequate to rationalize the observed spectroscopic anomalies. However, a multiconfigurational approach elucidates the disparate behaviors of these complexes. Multireference perturbation theory, based on complete active space self-consistent field computations, identifies Cu(I) in the case of the complex with the NCH3 containing ligands and a state with substantial Cu(II) contributions in the case of the complex with the S atom containing ligands. In contrast, DFT indicates Cu(I) in both scenarios.

Predicting the FCI Energy of Large Systems to Chemical Accuracy from Restricted Active Space Density Matrix Renormalization Group Calculations

We theoretically derive and validate with large scale simulations a remarkably accurate power law scaling of errors for the restricted active space density matrix renormalization group (DMRG-RAS) method [J. Phys. Chem. A 126, 9709] in electronic structure calculations. This yields a new extrapolation method, DMRG-RAS-X, which reaches chemical accuracy for strongly correlated systems such as the chromium dimer, dicarbon up to a large cc-pVQZ basis and even a large chemical complex such as the FeMoco with significantly lower computational demands than those of previous methods. The method is free of empirical parameters, performed robustly and reliably in all examples we tested, and has the potential to become a vital alternative method for electronic structure calculations in quantum chemistry and more generally for the computation of strong correlations in nuclear and condensed matter physics.

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.

Spin state and magnetic coupling in polynuclear Ni(II) complexes from density functional theory: is there an optimal amount of Fock exchange?

Reliable prediction of the ground-state spin and magnetic coupling constants in transition-metal complexes is a well-known challenge for density functional theory (DFT). One popular strategy for addressing this long-standing issue involves the modification of the fraction of Fock exchange in a hybrid functional. Here we explore the viability of this approach using three polynuclear metal-organic complexes based on a Ni4O4 cubane motif, having different ground state spin values (S = 0, 2, 4) owing to the use of different ligands. We systematically search for an optimum fraction of Fock exchange, across various global, range-separated, and double hybrid functionals. We find that for all functionals tested, at best there only exists a very narrow range of Fock exchange fractions which results in a correct prediction of the ground-state spin for all three complexes. The useful range is functional dependent, but general trends can be identified. Typically, at least two similar systems must be used in order to determine both an upper and lower limit of the optimal range. This is likely owing to conflicting demands of minimizing delocalization errors, which typically requires a higher percentage of Fock exchange, and addressing static correlation, which typically requires a lower one. Furthermore, we find that within the optimal range of Fock exchange, the sign and relative magnitude of Ni-Ni magnetic coupling constants are reasonably well reproduced, but there is still room for quantitative improvement in the prediction. Thus, the prediction of spin state and magnetic coupling in polynuclear complexes remains an ongoing challenge for DFT.

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.

Generation of Ru(III)-hypochlorite with resemblance to the heme-dependent haloperoxidase enzyme

The reaction of [(Me/BnTPEN)RuII(NCCH3)]2+ (BnTPEN = N1-benzyl-N1,N2,N2-tris(pyridine-2-ylmethyl)ethane-1,2-diamine and MeTPEN = N1-methyl-N1,N2,N2-tris(pyridine-2-ylmethyl)ethane-1,2-diamine) with mCPBA in the presence of chloride ions in CH3CN : H2O generated a novel (Me/BnTPEN)RuIII-OCl species at room temperature. This hypochlorite adduct could also be obtained by the direct reaction of NaOCl and HClO4 with (L)RuII complexes. The current study mimics the synthesis of a metal hypochlorite adduct in a similar fashion as in the heme-dependent haloperoxidase enzyme. As an electrophilic oxidant, the ruthenium hypochlorite adduct catalyzes hydrogen atom abstraction reactions of phenols and their derivatives.

Electronic Structure of the Complete Series of Gas-Phase Manganese Acetylacetonates by X-ray Absorption Spectroscopy

Metal centers in transition metal-ligand complexes occur in a variety of oxidation states causing their redox activity and therefore making them relevant for applications in physics and chemistry. The electronic state of these complexes can be studied by X-ray absorption spectroscopy, which is, however, due to the complex spectral signature not always straightforward. Here, we study the electronic structure of gas-phase cationic manganese acetylacetonate complexes Mn(acac)1-3+ using X-ray absorption spectroscopy at the metal center and ligand constituents. The spectra are well reproduced by multiconfigurational wave function theory, time-dependent density functional theory as well as parameterized crystal field and charge transfer multiplet simulations. This enables us to get detailed insights into the electronic structure of ground-state Mn(acac)1-3+ and extract empirical parameters such as crystal field strength and exchange coupling from X-ray excitation at both the metal and ligand sites. By comparison to X-ray absorption spectra of neutral, solvated Mn(acac)2,3 complexes, we also show that the effect of coordination on the L3 excitation energy, routinely used to identify oxidation states, can contribute about 40-50% to the observed shift, which for the current study is 1.9 eV per oxidation state.

New insights into pertinent Fe-complexes for the synthesis of iron via the instant polyol process

Chemically synthesized iron is in demand for biomedical applications due to its large saturation magnetization compared to iron oxides. The polyol process, suitable for obtaining Co and Ni particles and their alloys, is laborious in synthesizing Fe. The reaction yields iron oxides, and the reaction pathway remains unexplored. This study shows that a vicinal polyol, such as 1,2-propanediol, is suitable for obtaining Fe rather than 1,3-propanediol owing to the formation of a reducible Fe intermediate complex. X-ray absorption spectroscopy analysis reveals the ferric octahedral geometry and tetrahedral geometry in the ferrous state of the reaction intermediates in 1,2-propanediol and 1,3-propanediol, respectively. The final product obtained using a vicinal polyol is Fe with a γ-Fe2O3 shell, while the terminal polyol is favourable for Fe3O4. The distinct Fe-Fe and Fe-O bond lengths suggest the presence of a carboxylate group and a terminal alkoxide ligand in the intermediate of 1,2-propanediol. A large Fe-Fe bond distance suggests diiron complexes with bidentate carboxylate bridges. Prominent high-spin and low-spin states indicate the possibility of transition, which favors the reduction of iron ions in the reaction using 1,2-propanediol.

Spin-state energetics and magnetic anisotropy in penta-coordinated Fe(III) complexes with different axial and equatorial ligand environments

The penta-coordinated trigonal-bi-pyramidal (TBP) Fe(III) complex (PMe₂Ph)₂FeCl₃ shows a reduced magnetic anisotropy in its intermediate-spin (IS) state as compared to its methyl-analog (PMe₃)₂Fe(III)Cl₃. In this work, the ligand environment in (PMe₂Ph)₂FeCl₃ is systematically altered by replacing the axial -P with -N and -As, the equatorial -Cl with other halides, and the axial methyl group with an acetyl group. This has resulted in a series of Fe(III) TBP complexes modelled in their IS and high-spin (HS) states. Lighter ligands -N and -F stabilize the complex in the HS state, while the magnetically anisotropic IS state is stabilized by -P and -As at the axial site, and -Cl, -Br, and -I at the equatorial site. Larger magnetic anisotropies appear for complexes with nearly degenerate ground electronic states that are well separated from the higher excited states. This requirement, largely controlled by the d-orbital splitting pattern due to the changing ligand field, is achieved with a certain combination of axial and equatorial ligands, such as -P and -Br, -As and -Br, and -As and -I. In most cases, the acetyl group at the axial site enhances the magnetic anisotropy compared to its methyl counterpart. In contrast, the presence of -I at the equatorial site compromises the uniaxial type of anisotropy of the Fe(III) complex leading to an enhanced rate of quantum tunneling of magnetization.

Design and Mechanism of Rare-Earth Singlet Oxygen Sensing: An Experimental and Quantum Chemical Approach

The sensitive detection of singlet oxygen (¹O₂) is one key issue in various photochemical analyses, reactions, and processes; it is indispensable for designing catalysts for photodynamic therapies. Corresponding fluorescence-based organic ¹O₂ monitor luminophores may be equipped with rare-earth complexes with several intrinsic advantages. The design of the necessary ligands being a tedious, time-consuming effort, often involving empirical guesswork, we decided to support our experimental work with quantum chemical calculations. Hence, next to the experimental core, this paper suggests the additional use of time-dependent density functional theory (TDDFT) on suitable, free β-diketonate ligands to devise corresponding Eu³⁺ complexes as ¹O₂ probes eventually; the free ligand calculations obviously allow profoundly reduced computational efforts. Novel β-diketonate-substituted dimethyl anthracene complexes of Eu³⁺, Tb³⁺, and Gd³⁺ and their endoperoxidized descendants were thus synthesized, compared to known related complexes and analyzed with regard to their electronic characteristics; in addition, spectroscopy of a Eu³⁺ complex with ancillary epoxiphenanthroline for subsequent attachment to biological substrates featuring -NH₂ or -SH groups was included. The spectroscopic determination of the decisive lowest triplet (T₁) states of the Gd complexes could be matched by the Tamm-Dancoff approximation (TDA)/TDDFT calculations on the free ligands satisfactorily if suitable functionals were applied. Most significantly, the results suffice to describe the luminescence "switch-on" mechanism of this complex in the presence of ¹O₂.

Photodissociation Spectroscopy and Photofragment Imaging of the Fe+(Acetylene) Complex

Tunable laser photodissociation spectroscopy in the 700-400 nm region and photofragment imaging experiments are employed to investigate the Fe⁺(acetylene) ion-molecule complex. At energies above a threshold at 679 nm, continuous dissociation is detected throughout the visible wavelength region, with regions of broad structure. Comparison to the spectrum predicted by time-dependent density functional theory (TD-DFT) indicates that the complex has a quartet ground state. The dissociation threshold for Fe⁺(acetylene) at 679 nm provides the dissociation energy on the quartet potential energy surface. Correction for the atomic quartet-sextet spin state energy difference provides an adiabatic dissociation energy of 36.8 ± 0.2 kcal/mol. Photofragment imaging of the Fe⁺ photoproduct produced at 603.5 nm produces significant kinetic energy release (KER). The photon energy and the maximum value of the KER provide an upper limit on the dissociation energy of D₀ ≤ 34.6 ± 3.2 kcal/mol. The dissociation energies determined from the spectroscopy and photofragment imaging experiments agree nicely with the value determined previously by collision-induced dissociation (38.0 ± 2.6 kcal/mol). However, both values are significantly lower than those produced by computational chemistry at the DFT level using different functionals recommended for transition-metal chemistry.

Rational design of iron catalysts for C-X bond activation

We have quantum chemically studied the iron-mediated CX bond activation (X = H, Cl, CH₃ ) by d⁸ -FeL₄ complexes using relativistic density functional theory at ZORA-OPBE/TZ2P. We find that by either modulating the electronic effects of a generic iron-catalyst by a set of ligands, that is, CO, BF, PH₃ , BN(CH₃ )₂ , or by manipulating structural effects through the introduction of bidentate ligands, that is, PH₂ (CH₂ )_n PH₂ with n = 6-1, one can significantly decrease the reaction barrier for the CX bond activation. The combination of both tuning handles causes a decrease of the CH activation barrier from 10.4 to 4.6 kcal mol^-1 . Our activation strain and Kohn-Sham molecular orbital analyses reveal that the electronic tuning works via optimizing the catalyst-substrate interaction by introducing a strong second backdonation interaction (i.e., "ligand-assisted" interaction), while the mechanism for structural tuning is mainly caused by the reduction of the required activation strain because of the pre-distortion of the catalyst. In all, we present design principles for iron-based catalysts that mimic the favorable behavior of their well-known palladium analogs in the bond-activation step of cross-coupling reactions.

Reductive Coupling of Nitric Oxide by Cu(I): Stepwise Formation of Mono- and Dinitrosyl Species En Route to a Cupric Hyponitrite Intermediate

Transition-metal-mediated reductive coupling of nitric oxide (NO(g)) to nitrous oxide (N2O(g)) has significance across the fields of industrial chemistry, biochemistry, medicine, and environmental health. Herein, we elucidate a density functional theory (DFT)-supplemented mechanism of NO(g) reductive coupling at a copper-ion center, [(tmpa)CuI(MeCN)]+ (1) {tmpa = tris(2-pyridylmethyl)amine}. At -110 °C in EtOH (<-90 °C in MeOH), exposing 1 to NO(g) leads to a new binuclear hyponitrite intermediate [{(tmpa)CuII}2(μ-N2O22-)]2+ (2), exhibiting temperature-dependent irreversible isomerization to the previously characterized κ2-O,O'-trans-[(tmpa)2Cu2II(μ-N2O22-)]2+ (OOXray) complex. Complementary stopped-flow kinetic analysis of the reaction in MeOH reveals an initial mononitrosyl species [(tmpa)Cu(NO)]+ (1-(NO)) that binds a second NO molecule, forming a dinitrosyl species [(tmpa)CuII(NO)2] (1-(NO)2). The decay of 1-(NO)2 requires an available starting complex 1 to form a dicopper-dinitrosyl species hypothesized to be [{(tmpa)Cu}2(μ-NO)2]2+ (D) bearing a diamond-core motif, en route to the formation of hyponitrite intermediate 2. In contrast, exposing 1 to NO(g) in 2-MeTHF/THF (v/v 4:1) at <-80 °C leads to the newly observed transient metastable dinitrosyl species [(tmpa)CuII(NO)2] (1-(NO)2) prior to its disproportionation-mediated transformation to the nitrite product [(tmpa)CuII(NO2)]+. Our study furnishes a near-complete profile of NO(g) activation at a reduced Cu site with tripodal tetradentate ligation in two distinctly different solvents, aided by detailed spectroscopic characterization of metastable intermediates, including resonance Raman characterization of the new dinitrosyl and hyponitrite species detected.

Coordination and Spin States in Fe+(C2H2)n Complexes Studied with Selected-Ion Infrared Spectroscopy

Fe⁺(acetylene)_n ion-molecule complexes are produced in a supersonic molecular beam with pulsed laser vaporization. These ions are mass selected and studied with infrared photodissociation spectroscopy in the C-H stretching region, complemented by computational chemistry calculations. All C-H stretch vibrations are shifted to frequencies lower than the vibrations of isolated acetylene because of the charge transfer that occurs between the metal ion and the molecules. Complexes in the size range of n = 1-4 are found to have structures with individual acetylene molecules bound to the core metal ion via cation-π interactions. The coordination is completed with four ligands in a structure close to a distorted tetrahedron. Larger complexes in the range of n = 5-8 have external acetylene molecules solvating this n = 4 core ion via CH-π bonding to inner-shell ligands. DFT computations predict that quartet spin states are more stable for all complex sizes, but infrared spectra for quartet and doublet spin states are quite similar, precluding definitive determination of the spin states. There is no evidence for any of these complexes having acetylenes coupled into reacted structures. This is consistent with computed thermochemistry, which finds significant activation barriers to such reactions.

Chemoton 2.0: Autonomous Exploration of Chemical Reaction Networks

Fueled by advances in hardware and algorithm design, large-scale automated explorations of chemical reaction space have become possible. Here, we present our approach to an open-source, extensible framework for explorations of chemical reaction mechanisms based on the first-principles of quantum mechanics. It is intended to facilitate reaction network explorations for diverse chemical problems with a wide range of goals such as mechanism elucidation, reaction path optimization, retrosynthetic path validation, reagent design, and microkinetic modeling. The stringent first-principles basis of all algorithms in our framework is key for the general applicability that avoids any restrictions to specific chemical systems. Such an agile framework requires multiple specialized software components of which we present three modules in this work. The key module, Chemoton, drives the exploration of reaction networks. For the exploration itself, we introduce two new algorithms for elementary-step searches that are based on Newton trajectories. The performance of these algorithms is assessed for a variety of reactions characterized by a broad chemical diversity in terms of bonding patterns and chemical elements. Chemoton successfully recovers the vast majority of these. We provide the resulting data, including large numbers of reactions that were not included in our reference set, to be used as a starting point for further explorations and for future reference.

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.

Three Novel Thiazole-Arm Containing 1,3,4-Oxadiazole-Based [HS-HS] Fe(II) Dinuclear Complexes

Novel synthesis of 2,5-bis[(1,3-thiazol-2-ylmethyl)amino]methyl-1,3,4-oxadiazole (LThiazole) is reported, along with the preparation of three new dinuclear Fe(II) complexes with different counterions [FeII2(μ2-LThiazole)2](BF4)4·2CH3CN (1), [FeII2(μ2-LThiazole)2](ClO4)4 (2) and [FeII2(μ2-LThiazole)2](CF3SO3)4·2CH3CN (3). The obtained complexes were characterized by single-crystal X-ray crystallography, SQUID magnetometry and IR spectroscopy. The structure of the crystalline material was determined at 120 K for 1 and 3. The magnetic properties of all three complexes (1–3) were measured between 2–300 K and clearly show that all three complexes stay in the high-spin state over the measured temperature range.

CASPT2 molecular geometries of Fe(II) spin-crossover complexes

Using fully internally contracted (FIC)-CASPT2 analytical gradients, geometry optimizations of spin-crossover complexes are reported. This approach is tested on a series of Fe(II) complexes with different sizes, ranging from 13 to 61 atoms. A combination of active space and basis set choices are employed to investigate their role in determining reliable molecular geometries. The reported strategy demonstrates that a wave function-based level of theory can be used to optimize the geometries of metal complexes in reasonable times and enables one to treat the molecular geometry and electronic structure of the complexes using the same level of theory. For a series of smaller Fe(II) SCO complexes, strong field ligands in the LS state result in geometries with the largest differences between DFT and CASPT2; however, good agreement overall is observed between DFT and CASPT2. For the larger complexes, moderate sized basis sets yield geometries that compare well with DFT and available experimental data. We recommend using the (10e,12o) active space since convergence to a minimum structure was more efficient than with truncated active spaces despite having similar Fe-ligand bond distances.

Antiaromaticity-Promoted Radical Anion stability in α-vinyl Heterocyclics

As an electron-rich species, radical anions have a wide range of applications in organic synthesis. In addition, aromaticity is an essential concept in chemistry that has attracted considerable attention from...

Fate of Sc-Ion Interaction With Water: A Computational Study to Address Splitting Water Versus Solvating Sc Ion

An exhaustive study of Sc-ion interaction with water molecules in all its possible oxidation and spin states has been carried out to delineate the relative propensity of Sc ions toward solvation and water splitting. Potential energy surface analysis of the Sc-ion reaction with water molecules, topological analysis of bonds, and the effect of sequential solvation up to 6 water molecules have been examined. Calculated values showed good agreement with the available experimental results. Close-shell systems such as singlet mono- and tricationic Sc ions prefer to split the water molecules. In contrast, the open-shell systems such as triplet mono- and doublet dicationic Sc ions prefer to get solvated than split the water molecule. Topological analysis of electron density predicted the Sc^+/2+–water bond as a noncovalent bond while Sc³⁺–OH₂, Sc²⁺–OH, and Sc⁺–H bonds as partially covalent in nature. Energy decomposition analysis revealed that Sc ion–water interactions are driven by electrostatic energy followed by polarization energy. The current study reveals that transition metal catalysis can be one of the most effective tools to employ in water splitting, by properly tuning the electrons, spin, and ligands around the catalytic center.

Tetrahedral Cyclopentadienylmetal Carbonyl Clusters of Manganese and Chromium: A Theoretical Study

Tetranuclear Cp₄M₄(CO)₄ clusters have been synthesized for iron and vanadium but not for the intermediate first-row transition metals manganese and chromium. All of the low-energy structures of these "missing" Cp₄M₄(CO)₄ (M = Mn, Cr) species are shown by density functional theory to consist of a central M₄ tetrahedron with each of the four faces capped by a μ₃-CO group. The individual low-energy structures differ in their spin states and in their formal metal-metal bond orders along the six edges of their central M₄ tetrahedra. The two low-energy Cp₄Mn₄(μ₃-CO)₄ structures are a triplet structure with all Mn-Mn single bonds and a singlet structure with one Mn≡Mn triple bond and five Mn-Mn single bonds along the six tetrahedral edges. Related low-energy Cp₄Cr₄(μ₃-CO)₄ structures include a quintet structure with all Cr-Cr single bonds and a singlet structure with two Cr≡Cr triple bonds and four Cr-Cr single bonds. However, the potential energy surface of the Cp₄Cr₄(CO)₄ system is complicated by three other structures of comparable energies including two triplet structures and one quintet structure with various combinations of single, double, and triple chromium-chromium bonds in the central Cr₄ tetrahedron.

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.

Computational Discovery of Transition-metal Complexes: From High-throughput Screening to Machine Learning

Transition-metal complexes are attractive targets for the design of catalysts and functional materials. The behavior of the metal-organic bond, while very tunable for achieving target properties, is challenging to predict and necessitates searching a wide and complex space to identify needles in haystacks for target applications. This review will focus on the techniques that make high-throughput search of transition-metal chemical space feasible for the discovery of complexes with desirable properties. The review will cover the development, promise, and limitations of "traditional" computational chemistry (i.e., force field, semiempirical, and density functional theory methods) as it pertains to data generation for inorganic molecular discovery. The review will also discuss the opportunities and limitations in leveraging experimental data sources. We will focus on how advances in statistical modeling, artificial intelligence, multiobjective optimization, and automation accelerate discovery of lead compounds and design rules. The overall objective of this review is to showcase how bringing together advances from diverse areas of computational chemistry and computer science have enabled the rapid uncovering of structure-property relationships in transition-metal chemistry. We aim to highlight how unique considerations in motifs of metal-organic bonding (e.g., variable spin and oxidation state, and bonding strength/nature) set them and their discovery apart from more commonly considered organic molecules. We will also highlight how uncertainty and relative data scarcity in transition-metal chemistry motivate specific developments in machine learning representations, model training, and in computational chemistry. Finally, we will conclude with an outlook of areas of opportunity for the accelerated discovery of transition-metal complexes.

Controlling Radical-Type Single-Electron Elementary Steps in Catalysis with Redox-Active Ligands and Substrates

Advances in (spectroscopic) characterization of the unusual electronic structures of open-shell cobalt complexes bearing redox-active ligands, combined with detailed mapping of their reactivity, have uncovered several new catalytic radical-type protocols that make efficient use of the synergistic properties of redox-active ligands, redox-active substrates, and the metal to which they coordinate. In this perspective, we discuss the tools available to study, induce, and control catalytic radical-type reactions with redox-active ligands and/or substrates, contemplating recent developments in the field, including some noteworthy tools, methods, and reactions developed in our own group. The main topics covered are (i) tools to characterize redox-active ligands; (ii) novel synthetic applications of catalytic reactions that make use of redox-active carbene and nitrene substrates at open-shell cobalt-porphyrins; (iii) development of catalytic reactions that take advantage of purely ligand- and substrate-based redox processes, coupled to cobalt-centered spin-changing events in a synergistic manner; and (iv) utilization of redox-active ligands to influence the spin state of the metal. Redox-active ligands have emerged as useful tools to generate and control reactive metal-coordinated radicals, which give access to new synthetic methodologies and intricate (electronic) structures, some of which are yet to be exposed.

Mechanistic Insights into the Oxygen Atom Transfer Reactions by Nonheme Manganese Complex: A Computational Case Study on the Comparative Oxidative Ability of Manganese-Hydroperoxo vs High-Valent MnIV═O and MnIV-OH Intermediates

Understanding the comparative oxidative abilities of high-valent metal-oxo/hydroxo/hydroperoxo species holds the key to robust biomimic catalysts that perform desired organic transformations with very high selectivity and efficiency. The comparative oxidative abilities of popular high-valent iron-oxo and manganese-oxo species are often counterintuitive, for example, oxygen atom transfer (OAT) reaction by [(Me₂EBC)Mn^IV-OOH]³⁺, [(Me₂EBC)Mn^IV-OH]³⁺, and [(Me₂EBC)Mn^IV═O]²⁺ (Me₂EBC = 4,11-dimethyl-1,4,8,11-tetraazabicyclo[6.6.2]hexadecane) shows extremely high reactivity for Mn^IV-OOH species and no reactivity for Mn^IV-OH and Mn^IV═O species toward alkyl/aromatic sulfides. Using a combination of density functional theory (DFT) and ab initio domain-based local pair natural orbital coupled-cluster with single, double, and perturbative triples excitation (DLPNO-CCSD(T)) and complete-active space self-consistent field/N-electron valence perturbation theory second order (CASSCF/NEVPT2) calculations, here, we have explored the electronic structures and sulfoxidation mechanism of these species. Our calculations unveil that Mn^IV-OOH reacts through distal oxygen atom with the substrate via electron transfer (ET) mechanism with a very small kinetic barrier (16.5 kJ/mol), placing this species at the top among the best-known catalysts for such transformations. The Mn^IV-OH and Mn^IV═O species have a much larger barrier. The mechanism has also been found to switch from ET in the former to concerted in the latter, rendering both unreactive under the tested experimental conditions. Intrinsic differences in the electronic structures, such as the presence and absence of the multiconfigurational character coupled with the steric effects, are responsible for such variations observed. This comparative oxidative ability that runs contrary to the popular iron-oxo/hydroperoxo reactivity will have larger mechanistic implications in understanding the reactivity of biomimic catalysts and the underlying mechanisms in PSII.

Is the Electrophilicity of the Metal Nitrene the Sole Predictor of Metal-Mediated Nitrene Transfer to Olefins? Secondary Contributing Factors as Revealed by a Library of High-Spin Co(II) Reagents

Recent research has highlighted the key role played by the electron affinity of the active metal-nitrene/imido oxidant as the driving force in nitrene additions to olefins to afford valuable aziridines. The present work showcases a library of Co(II) reagents that, unlike the previously examined Mn(II) and Fe(II) analogues, demonstrate reactivity trends in olefin aziridinations that cannot be solely explained by the electron affinity criterion. A family of Co(II) catalysts (17 members) has been synthesized with the assistance of a trisphenylamido-amine scaffold decorated by various alkyl, aryl, and acyl groups attached to the equatorial amidos. Single-crystal X-ray diffraction analysis, cyclic voltammetry and EPR data reveal that the high-spin Co(II) sites (S = 3/2) feature a minimal [N3N] coordination and span a range of 1.4 V in redox potentials. Surprisingly, the Co(II)-mediated aziridination of styrene demonstrates reactivity patterns that deviate from those anticipated by the relevant electrophilicities of the putative metal nitrenes. The representative L4Co catalyst (-COCMe3 arm) is operating faster than the L8Co analogue (-COCF3 arm), in spite of diminished metal-nitrene electrophilicity. Mechanistic data (Hammett plots, KIE, stereocontrol studies) reveal that although both reagents follow a two-step reactivity path (turnover-limiting metal-nitrene addition to the C b atom of styrene, followed by product-determining ring-closure), the L4Co catalyst is associated with lower energy barriers in both steps. DFT calculations indicate that the putative [L4Co]NTs and [L8Co]NTs species are electronically distinct, inasmuch as the former exhibits a single-electron oxidized ligand arm. In addition, DFT calculations suggest that including London dispersion corrections for L4Co (due to the polarizability of the tert-Bu substituent) can provide significant stabilization of the turnover-limiting transition state. This study highlights how small ligand modifications can generate stereoelectronic variants that in certain cases are even capable of overriding the preponderance of the metal-nitrene electrophilicity as a driving force.

On the Use of Normalized Metrics for Density Sensitivity Analysis in DFT

In the past years, there has been a discussion about how the errors in density functional theory might be related to errors in the self-consistent densities obtained from different density functional approximations. This, in turn, brings up the discussion about the different ways in which we can measure such errors and develop metrics that assess the sensitivity of calculated energies to changes in the density. It is important to realize that there cannot be a unique metric in order to look at this density sensitivity, simultaneously needing size-extensive and size-intensive metrics. In this study, we report two metrics that are widely applicable to any density functional approximation. We also show how they can be used to classify different chemical systems of interest with respect to their sensitivity to small variations in the density.

A Pseudotetrahedral Terminal Oxoiron(IV) Complex: Mechanistic Promiscuity in C-H Bond Oxidation Reactions

S=2 oxoiron(IV) species act as reactive intermediates in the catalytic cycle of nonheme iron oxygenases. The few available synthetic S=2 FeIV =O complexes known to date are often limited to trigonal bipyramidal and very rarely to octahedral geometries. Herein we describe the generation and characterization of an S=2 pseudotetrahedral FeIV =O complex 2 supported by the sterically demanding 1,4,7-tri-tert-butyl-1,4,7-triazacyclononane ligand. Complex 2 is a very potent oxidant in hydrogen atom abstraction (HAA) reactions with large non-classical deuterium kinetic isotope effects, suggesting hydrogen tunneling contributions. For sterically encumbered substrates, direct HAA is impeded and an alternative oxidative asynchronous proton-coupled electron transfer mechanism prevails, which is unique within the nonheme oxoiron community. The high reactivity and the similar spectroscopic parameters make 2 one of the best electronic and functional models for a biological oxoiron(IV) intermediate of taurine dioxygenase (TauD-J).

Zeolites at the Molecular Level: What Can Be Learned from Molecular Modeling

This review puts the development of molecular modeling methods in the context of their applications to zeolitic active sites. We attempt to highlight the utmost necessity of close cooperation between theory and experiment, resulting both in advances in computational methods and in progress in experimental techniques.

Tuning spin-crossover transition temperatures in non-symmetrical homoleptic meridional/facial [Fe(didentate)3]2+ complexes: what for and who cares about it?

The [FeN₆] chromophores found in [Fe(didentate)₃]²⁺ complexes, where didentate is a non-symmetrical 2-(6-membered-heterocyclic ring)-benzimidazole ligand (Lk), exist as mixtures of two geometrical mer (C₁-symmetry) and fac (C₃-symmetry) isomers. Specific alkyl-substituted six-membered heterocyclic rings connected to the benzimidazole unit (pyridines in ligands L1-L3, pyrazines in L4-L5 and pyrimidines in L6-L7) control the ligand field strength and the electron delocalization so that [Fe^II(Lk)₃]²⁺ display tunable thermally-induced spin transitions in solution. Thermodynamic, spectroscopic (UV-Vis, NMR) and magnetic studies in solution demonstrate that [Fe(L6)₃]²⁺ (L6 = 1-methyl-2-(pyrimidin-2-yl)-1H-benzo[d]imidazole) exhibits a close to room temperature spin transition (T_1/2 = 273(3) K) combined with a high stability formation constant ( in acetonitrile), which makes this complex suitable for the potential modulation of lanthanide-based luminescence in polymetallic helicates. A novel method is proposed for assigning specific thermodynamic spin crossover parameters to fac-[Fe(L6)₃]²⁺ and mer-[Fe(L6)₃]²⁺ isomers in solution. The observed difference relies mainly on the entropic content ΔS-ΔS = 11(1) J mol^-1 K^-1, which favors the spin transition for the meridional isomer. Intermolecular interactions occurring in the crystalline state largely overcome minor thermodynamic trends operating in diluted solutions and a single configurational isomer is usually observed in the solid state. Among the thirteen solved crystal structures 1-13 containing the [M(Lk)₃]²⁺ cations (M = Fe, Ni, Zn, Lk = L6-L7), pure meridional isomers are observed six times, pure facial isomers also six times and a mixture (44% mer and 56% fac) is detected only once. Solid-state magnetic data recorded for the Fe^II complexes show the operation of slightly cooperative spin transitions in 7 (fac-[Fe(L6)₃]²⁺) and 12 (mer-[Fe(L7)₃]²⁺). For the meridional isomer in 6, a two-step spin state transition curve, associated with two phase transitions, is detected.

Widening the Window of Spin-Crossover Temperatures in Bis(formazanate)iron(II) Complexes via Steric and Noncovalent Interactions

Bis(formazanate)iron(II) complexes undergo a thermally induced S = 0 to S = 2 spin transition in solution. Here we present a study of how steric effects and π-stacking interactions between the triarylformazanate ligands affect the spin-crossover behavior, in addition to electronic substituent effects. Moreover, the effect of increasing the denticity of the formazanate ligands is explored by including additional OMe donors in the ligand (7). In total, six new compounds (2–7) have been synthesized and characterized, both in solution and in the solid state, via spectroscopic, magnetic, and structural analyses. The series spans a broad range of spin-crossover temperatures (T_1/2) for the LS ⇌ HS equilibrium in solution, with the exception of compound 6 which remains high-spin (S = 2) down to 210 K. In the solid state, 6 was shown to exist in two distinct forms: a tetrahedral high-spin complex (6a, S = 2) and a rare square-planar structure with an intermediate-spin state (6b, S = 1). SQUID measurements, ⁵⁷Fe Mössbauer spectroscopy, and differential scanning calorimetry indicate that in the solid state the square-planar form 6b undergoes an incomplete spin-change-coupled isomerization to tetrahedral 6a. The complex that contains additional OMe donors (7) results in a six-coordinate (NNO)₂Fe coordination geometry, which shifts the spin-crossover to significantly higher temperatures (T_1/2 = 444 K). The available experimental and computational data for 7 suggest that the Fe···OMe interaction is retained upon spin-crossover. Despite the difference in coordination environment, the weak OMe donors do not significantly alter the electronic structure or ligand-field splitting, and the occurrence of spin-crossover (similar to the compounds lacking the OMe groups) originates from a large degree of metal–ligand π-covalency.

A series of Fe(II) complexes with formazanate ligands are reported, and ligand substituent effects on structure and spin-crossover properties are examined. These ligand modifications allow isolation of compounds with tetrahedral geometries in both low- and high-spin ground states as well as an intermediate-spin square-planar complex. Steric properties, π-stacking interactions, and additional donor substituents lead to a wide range of spin-crossover temperatures (T_1/2) in this class of compounds.

Spin-state energetics of metallocenes: How do best wave function and density functional theory results compare with the experimental data?

We benchmark the accuracy of quantum-chemical methods, including wave function theory methods [coupled cluster theory at the CCSD(T) level, multiconfigurational perturbation-theory (CASPT2, NEVPT2) and internally contracted multireference configuration interaction (MRCI)] and 30 density functional theory (DFT) approximations, in reproducing the spin-state splittings of metallocenes. The reference values of the electronic energy differences are derived from the experimental spin-crossover enthalpy for manganocene and the spectral data of singlet-triplet transitions for ruthenocene, ferrocene, and cobaltocenium. For ferrocene and cobaltocenium we revise the previous experimental interpretations regarding the lowest triplet energy; our argument is based on the comparison with the lowest singlet excitation energy and herein reported, carefully determined absorption spectrum of ferrocene. When deriving vertical energies from the experimental band maxima, we go beyond the routine vertical energy approximation by introducing vibronic corrections based on simulated vibrational envelopes. The benchmarking result confirms the high accuracy of the CCSD(T) method (in particular, for UCCSD(T) based on Hartree-Fock orbitals we find for our dataset: maximum error 0.12 eV, weighted mean absolute error 0.07 eV, weighted mean signed error 0.01 eV). The high accuracy of the single-reference method is corroborated by the analysis of a multiconfigurational character of the complete active space wave function for the triplet state of ferrocene. On the DFT side, our results confirm the non-universality problem with approximate functionals. The present study is an important step toward establishing an extensive and representative benchmark set of experiment-derived spin-state energetics for transition metal complexes.

Closed Shell Iron(IV) Oxo Complex with an Fe-O Triple Bond: Computational Design, Synthesis, and Reactivity

Iron(IV)-oxo intermediates in nature contain two unpaired electrons in the Fe-O antibonding orbitals, which are thought to contribute to their high reactivity. To challenge this hypothesis, we designed and synthesized closed-shell singlet iron(IV) oxo complex [(quinisox)Fe(O)]+ (1+ ; quinisox-H=(N-(2-(2-isoxazoline-3-yl)phenyl)quinoline-8-carboxamide). We identified the quinisox ligand by DFT computational screening out of over 450 candidates. After the ligand synthesis, we detected 1+ in the gas phase and confirmed its spin state by visible and infrared photodissociation spectroscopy (IRPD). The Fe-O stretching frequency in 1+ is 960.5 cm-1 , consistent with an Fe-O triple bond, which was also confirmed by multireference calculations. The unprecedented bond strength is accompanied by high gas-phase reactivity of 1+ in oxygen atom transfer (OAT) and in proton-coupled electron transfer reactions. This challenges the current view of the spin-state driven reactivity of the Fe-O complexes.

Improvement of d-d interactions in density functional tight binding for transition metal ions with a ligand field model: assessment of a DFTB3+U model on nickel coordination compounds

To improve the description of interactions among the localized d, f electrons in transition metals, we have introduced a ligand-field motivated contribution into the Density Functional Tight Binding (DFTB) model. Referred to as DFTB3+U, the approach treats the d, f electron repulsions with rotationally invariant orbital-orbital interactions and a Hartree-Fock model; this represents a major conceptual improvement over the original DFTB3 approach, which treats the d, f-shell interactions in a highly averaged fashion without orbital level of description. The DFTB3+U approach is tested using a series of nickel compounds that feature Ni(ii) and Ni(iii) oxidation states. By using parameters developed with the original DFTB3 Hamiltonian and empirical +U parameters (F0/2/4 Slater integrals), we observe that the DFTB3+U model indeed provides substantial improvements over the original DFTB3 model for a number of properties of the nickel compounds, including the population and spin polarization of the d-shell, nature of the frontier orbitals, ligand field splitting and the energy different between low and high spin states at OPBE optimized structures. This proof-of-concept study suggests that with self-consistent parameterization of the electronic and +U parameters, the DFTB3+U model can develop into a promising model that can be used to efficiently study reactive events involving transition metals ion condensed phase systems. The methodology can be integrated with other approximate QM methods as well, such as the extended tight binding (xTB) approach.