Two-State Reactivity in Organometallic Gas-Phase Ion Chemistry

The spin polarization strategy regulates heterogeneous catalytic activity performance: from fundamentals to applications

In recent years, there has been significant attention towards the development of catalysts that exhibit superior performance and environmentally friendly attributes. This surge in interest is driven by the growing demands for energy utilization and storage as well as environmental preservation. Spin polarization plays a crucial role in catalyst design, comprehension of catalytic mechanisms, and reaction control, offering novel insights for the design of highly efficient catalysts. However, there are still some significant research gaps in the current study of spin catalysis. Therefore, it is urgent to understand how spin polarization impacts catalytic reactions to develop superior performance catalysts. Herein, we present a comprehensive summary of the application of spin polarization in catalysis. Firstly, we summarize the fundamental mechanism of spin polarization in catalytic reactions from two aspects of kinetics and thermodynamics. Additionally, we review the regulation mechanism of spin polarization in various catalytic applications and several approaches to modulate spin polarization. Moreover, we discuss the future development of spin polarization in catalysis and propose several potential avenues for further progress. We aim to improve current catalytic systems through implementing a novel and distinctive spin engineering strategy.

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.

How Do Differences in Electronic Structure Affect the Use of Vanadium Intermediates as Mimics in Nonheme Iron Hydroxylases?

We study active-site models of nonheme iron hydroxylases and their vanadium-based mimics using density functional theory to determine if vanadyl is a faithful structural mimic. We identify crucial structural and energetic differences between ferryl and vanadyl isomers owing to the differences in their ground electronic states, i.e., high spin (HS) for Fe and low spin (LS) for V. For the succinate cofactor bound to the ferryl intermediate, we predict facile interconversion between monodentate and bidentate coordination isomers for ferryl species but difficult rearrangement for vanadyl mimics. We study isomerization of the oxo intermediate between axial and equatorial positions and find the ferryl potential energy surface to be characterized by a large barrier of ca. 10 kcal/mol that is completely absent for the vanadyl mimic. This analysis reveals even starker contrasts between Fe and V in hydroxylases than those observed for this metal substitution in nonheme halogenases. Analysis of the relative bond strengths of coordinating carboxylate ligands for Fe and V reveals that all of the ligands show stronger binding to V than Fe owing to the LS ground state of V in contrast to the HS ground state of Fe, highlighting the limitations of vanadyl mimics of native nonheme iron hydroxylases.

Investigation of novel radical scavenging mechanisms in the alkaline environment: Green, sustainable and environmentally friendly antioxidative agent(s)

Pharmaceutical and industrial utilization of synthetic chemicals has an immerse impact on the environment. In that sense, novel chemicals with potential for industrial application should be investigated for their behaviour in reactions with hydroxyl radical, simulating AOPs (Advanced Oxidation Processes). AOPs are known for being highly effective in wastewater management and natural water remediation. In this paper, exhaustive research on the radical scavenging activity of a newly synthesized coumarin derivative (4HCBH), as a representative of the series of coumarin-benzohydrazides with high antioxidative potential was conducted. This study took into consideration the pH value range significant for practically all living organisms (pH = 7.0-8.5). According to the experimentally obtained results, the 4HCBH showed an increase in radical scavenging activity, following the slight increase in pH values, which suggested that the formation of anionic form of 4HCBH is responsible for its antiradical activity. Further investigations led to the postulation of a novel mechanistic approach called Sequential Proton Loss Electron Transfer - Radical-Radical Coupling (SPLET-RRC), in which, by a series of steps, a new, stable compound was formed. Furthermore, it was demonstrated that the product generated through SPLET-RRC showed lower toxicity than the parent molecule.

Ion-molecule reactions of mass-selected ions

Gas-phase reactions of mass-selected ions with neutrals covers a very broad area of fundamental and applied mass spectrometry (MS). Oftentimes, ion-molecule reactions (IMR) can serve as a viable alternative to collision-induced dissociation and other ion dissociation techniques when using tandem MS. This review focuses on the literature pertaining applications of IMR since 2013. During the past decade considerable efforts have been made in analytical applications of IMR, including advances in one of the major techniques for characterization of unsaturated fatty acids and lipids, ozone-induced dissociation, and the development of a new technique for sequencing of large ions, hydrogen atom attachment/abstraction dissociation. Many advances have also been made in identifying gas-phase chemistry specific to a functional group in organic and biological compounds, which are useful in structure elucidation of analytes and differentiation of isomers/isobars. With "soft" ionization techniques like electrospray ionization having become mainstream for quite some time now, the efforts in the area of metal ion catalysis have firmly moved into exploring chemistry of ligated metal complexes in their "natural" oxidation states allowing to model individual steps of mechanisms in homogeneous catalysis, especially in combination with high-level DFT calculations. Finally, IMR continue to contribute to the body of knowledge in the area of chemistry of interstellar processes.

Observation of inserted oxocarbonyl species in the tantalum cation-mediated activation of carbon dioxide dictated by two-state reactivity

Reductive activation of carbon dioxide (CO2) has drawn increasing attention as an effective and convenient method to unlock this stable molecule, especially via transition metal-catalyzed reactions. Taking the [TaC4O8]+ ion-molecule complex formed in the laser ablation source as a representative, the reactivity of the tantalum metal cation towards CO2 molecules is explored using infrared photodissociation spectroscopy combined with quantum chemical calculations. The strong absorption in the carbonyl stretching region provides solid evidence for the insertion reactions into CO bonds by the tantalum cation. Two inserted oxocarbonyl products are identified based on the great agreement between the experimental results and simulated infrared spectra of energetically low-lying structures in the singlet and triplet states. The pivotal role of two-state reactivity in driving CO2 activation among three different spin states is rationalized by potential energy surface analysis. Our conclusion provides valuable insight into the intrinsic mechanisms of CO2 activation by the tantalum metal cation, highlighting the affinity of tantalum for CO bond insertion in addition to typical "end-on" binding configurations.

Quantitative Aspects of Gas-Phase Metal Ion Chemistry: Conservation of Spin, Participation of f Orbitals, and C-H Activation and C-C Coupling

In this Featured Article, I reflect on over 40 years of guided ion beam tandem mass spectrometry (GIBMS) studies involving atomic metal cations and their clusters throughout the periodic table. Studies that have considered the role of spin conservation (or lack thereof) are a primary focus with a quantitative assessment of the effects examined. A need for state-specific studies of heavier elements is noted, as is a more quantitative assessment of spin-orbit interactions in reactivity. Because GIBMS experiments explicitly evaluate the kinetic energy dependence of reactions over a wide range, several interesting and unusual observations are highlighted. More detailed studies of such unusual reaction events would be welcome. Activation of C-H bonds and ensuing C-C coupling events are reviewed, with future work encouraged. Finally, studies of lanthanides and actinides are examined with an eye on understanding the role of f orbitals in the chemistry, both as participants (or not) in the bonding and as sources/sinks of electron density. This area seems to be ripe for more quantitative experiments.

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).

Two state reactivity (TSR) and hydrogen tunneling reaction kinetics measured in the Co+ mediated decomposition of CH3CHO

The electronic structure of a transition metal atom allows it to act as a catalytic active site by providing lower energy alternative pathways in chemical transformations. We have identified and kinetically characterized three such pathways in the title reaction. One is an adiabatic pathway that occurs on a single potential energy surface described within the Born-Oppenheimer approximation. A second pathway opens microseconds into the reaction as a portion of the reacting population competitively transitions from triplet to singlet multiplicity to circumvent energetic barriers on the triplet surface. These pathways are single- and two-state reactive (SSR and TSR) where the Co+ cation mediates an oxidative addition/reductive elimination sequence of the CH3CHO molecule. The third observed reaction pathway is the aldehyde hydrogen tunneling through an Eyring barrier to form high-spin products. First-order rate constants for the adiabatic and nonadiabatic energy lowered pathways, and the hydrogen tunneling pathway, are each measured using the single photon initiated dissociative rearrangement reaction (SPIDRR) experimental technique. We believe that this is the first experimental study where such disparate dynamic features (SSR, TSR, and H-tunneling) are disentangled in a system's chemistry, attributing specific rate constant values to each effect and quantifying the various competitions. Moreover, multi-reference CASSCF/CASPT2 calculations indicate that structures with covalent Co-H bonds are present exclusively along the excited singlet surface. This phenomenon significantly reduces these structures' energy relative to their triplet counterparts, thus enabling the surface crossing and spin inversion that cause the observed two-state reactivity.

A general method for locating stationary points on the mixed-spin surface of spin-forbidden reaction with multiple spin states

Some chemical reactions proceed on multiple potential energy surfaces and are often accompanied by a change in spin multiplicity, being called spin-forbidden reactions, where the spin-orbit coupling (SOC) effects play a crucial role. In order to efficiently investigate spin-forbidden reactions with two spin states, Yang et al. [Phys. Chem. Chem. Phys. 20, 4129-4136 (2018)] proposed a two-state spin-mixing (TSSM) model, where the SOC effects between the two spin states are simulated by a geometry-independent constant. Inspired by the TSSM model, we suggest a multiple-state spin-mixing (MSSM) model in this paper for the general case with any number of spin states, and its analytic first and second derivatives have been developed for locating stationary points on the mixed-spin potential energy surface and estimating thermochemical energies. To demonstrate the performance of the MSSM model, some spin-forbidden reactions involving 5d transition elements are calculated using the density functional theory (DFT), and the results are compared with the two-component relativistic ones. It is found that MSSM DFT and two-component DFT calculations may provide very similar stationary-point information on the lowest mixed-spin/spinor energy surface, including structures, vibrational frequencies, and zero-point energies. For the reactions containing saturated 5d elements, the reaction energies by MSSM DFT and two-component DFT agree very well within 3 kcal/mol. As for the two reactions OsO+ + CH4 → OOs(CH2)+ + H2 and W + CH4 → WCH2 + H2 involving unsaturated 5d elements, MSSM DFT may also yield good reaction energies of similar accuracy but with some counterexamples. Nevertheless, the energies may be remarkably improved by a posteriori single point energy calculations using two-component DFT at the MSSM DFT optimized geometries, and the maximum error of about 1 kcal/mol is almost independent of the SOC constant used. The MSSM method as well as the developed computer program provides an effective utility for studying spin-forbidden reactions.

Reactivities of iron(IV)-oxido compounds with pentadentate bispidine ligands

The Fe^IVO complexes of bispidines (3,7-diazabicyclo[3.3.1]nonane derivatives) are known to be highly reactive oxidants - with the tetradentate bispidine, the so far most reactive ferryl complex has been reported and two isomeric pentadentate ligands also lead to very reactive high-valent oxidants. With a series of 4 new bispidine derivatives we now try to address the question why the bispidine scaffold in general leads to very reactive oxidants and how this can be tuned by ligand modifications. The study is based on a full structural, spectroscopic and electrochemical analysis of the iron(II) precursors, spectroscopic data of the iron(IV)-oxido complexes, a kinetic analysis of the stoichiometric oxidation of thioanisole by five different bispidine‑iron(IV)-oxido complexes and on product analyses of reactions by the five ferryl oxidants with thioanisole, β-methylstyrene and cis-stilbene as substrates.

Direct Coupling of Methane and Carbon Dioxide on Tantalum Cluster Cations

Understanding molecular-scale reaction mechanisms is crucial for the design of modern catalysts with industrial prospect. Through joint experimental and computational studies, we investigate the direct coupling reaction of CH₄ and CO₂ , two abundant greenhouse gases, mediated by Ta_1,4 ⁺ ions to form larger oxygenated hydrocarbons. Coherent with proposed elementary steps, we expose products of CH₄ dehydrogenation [Ta_1,4 CH₂ ]⁺ to CO₂ in a ring electrode ion trap. Product analysis and reaction kinetics indicate a predisposition of the tetramers for C-O coupling with a conversion to products of CH₂ O, whereas atomic cations enable C-C coupling yielding CH₂ CO. Selected experimental findings are supported by thermodynamic computations, connecting structure, electronic properties, and catalyst function. Moreover, the study of bare Ta_1,4 ⁺ compounds indicates that methane dehydrogenation is a significant initial step in the direct coupling reaction, enabling new, yet unknown reaction pathways.

Mechanistic insights into aerobic oxidative cleavage of glycol catalyzed by an Anderson-type polyoxometalate [IMo6O24]5

A computational investigation of the aerobic oxidative C-C bond cleavage reaction of glycol catalyzed by an Anderson-type heteropolyanion HPA [IMo₆O₂₄]^5- in the presence of acetonitrile as solvent has been performed at the WB97XD/6-31G(d,p)/lanl2dz level. Two reaction pathways have been identified. The catalytic cycle of each pathway consists of three steps: oxidation cleavage of a glycol molecule by the HPA, oxidation of the HPA by one dioxygen molecule, and, finally, oxidation of a second glycol and regeneration of the catalyst. These reaction pathways have been thoroughly investigated in terms of energetic, natural bond orbital (NBO), natural charges, and geometrical parameters. It is found that (i) even though the top oxygen atoms of the Anderson heteropolyanion are not the most negatively charged ones, they are more likely to react with the diol hydroxyl groups, (ii) a direct relationship between the presence of the iodine ion I(VII) and the studied oxidation reaction could not be identified, and (iii) in terms of energy, the transfer of the two hydrogen atoms is the most energetic step.

Spin-State Splittings in 3d Transition-Metal Complexes Revisited: Toward a Reliable Theory Benchmark

A new composite method for the calculation of spin-crossover energies in 3d transition-metal complexes based on multireference methods is presented. The method reduces to MRCISD+Q at the complete-basis-set (CBS) level for atomic ions, for which it gives excitation energies with a mean absolute error of only ca. 0.01 eV. For molecular complexes, the CASPT2+δMRCI composite approach corresponds to a CASPT2/CBS calculation augmented by a high-level MRCISD+Q-CASPT2 correction with a smaller ligand basis set. For a set of [Fe(He)6]n+ test complexes, the approach reproduces full MRCISD+Q/CBS results to within better than 0.04 eV, without depending on any arbitrary IPEA shifts. The high-quality CASPT2+δMRCI method has then been applied to a series of 3d transition-metal hexaqua complexes in aqueous solution, augmented by an elaborate 3D-RISM-SCF solvent treatment of the underlying structures. It provides unprecedented agreement with experiment for the lowest-lying vertical spin-flip excitation energies, except for the Fe³⁺ system. Closer examination of the latter case provides strong evidence that the observed lowest-energy excitation at 1.56 eV, which has been used frequently for evaluating quantum-chemical methods, does not arise from the iron(III) hexaqua complex in solution, but from its singly deprotonated counterpart, [Fe(H2O)5OH]2+.

Spin-State Splittings in 3d Transition-Metal Complexes Revisited: Benchmarking Approximate Methods for Adiabatic Spin-State Energy Differences in Fe(II) Complexes

The CASPT2+δMRCI composite approach reported in a companion paper has been extended and used to provide high-quality reference data for a series of adiabatic spin gaps (defined as ΔE = E_quintet - E_singlet) of [Fe^IIL₆]²⁺ complexes (L = CNH, CO, NCH, NH₃, H₂O), either at nonrelativistic level or including scalar relativistic effects. These highly accurate data have been used to evaluate the performance of various more approximate methods. Coupled-cluster theory with singles, doubles, and perturbative triples, CCSD(T), is found to agree well with the new reference data for Werner-type complexes but exhibits larger underestimates by up to 70 kJ/mol for the π-acceptor ligands, due to appreciable static correlation in the low-spin states of these systems. Widely used domain-based local CCSD(T) calculations, DLPNO-CCSD(T), are shown to depend very sensitively on the cutoff values used to construct the localized domains, and standard values are not sufficient. A large number of density functional approximations have been evaluated against the new reference data. The B2PLYP double hybrid gives the smallest deviations, but several functionals from different rungs of the usual ladder hierarchy give mean absolute deviations below 20 kJ/mol. This includes the B97-D semilocal functional, the PBE0* global hybrid with 15% exact-exchange admixture, as well as the local hybrids LH07s-SVWN and LH07t-SVWN. Several further functionals achieve mean absolute errors below 30 kJ/mol (M06L-D4, SSB-D, B97-1-D4, LC-ωPBE-D4, LH12ct-SsirPW92-D4, LH12ct-SsifPW92-D4, LH14t-calPBE-D4, LHJ-HFcal-D4, and several further double hybrids) and thereby also still overall outperform CCSD(T) or uncorrected CASPT2. While exact-exchange admixture is a crucial factor in favoring high-spin states, the present evaluations confirm that other aspects can be important as well. A number of the better-performing functionals underestimate the spin gaps for the π-acceptor ligands but overestimate them for L = NH₃, H₂O. In contrast to a previous suggestion, non-self-consistent density functional theory (DFT) computations on top of Hartree-Fock orbitals are not a promising path to produce accurate spin gaps in such complexes.

Facile Addition of B-H and B-B Bonds to an Iron(IV) Nitride Complex

The nitride ligand in the iron(IV) complex PhB(ⁱPr₂Im)₃Fe≡N reacts with boron hydrides to afford PhB(ⁱPr₂Im)₃FeN(B)H (B = 9-BBN (1), Bpin (2)) and with (Bpin)₂ to afford PhB(ⁱPr₂Im)₃FeN(Bpin)₂ (3). The iron(II) borylamido products have all been structurally and spectroscopically characterized, demonstrating facile insertion into B-H and B-B bonds by PhB(ⁱPr₂Im)₃Fe≡N. Density functional theory (DFT) calculations reveal that the quintet state (S = 2) is significantly lower in energy than the singlet (S = 0) and triplet (S = 1) states for all products. Stoichiometric reaction with (Bpin)₂ does not produce the mono-borylated iron imido species PhB(ⁱPr₂Im)₃FeN(Bpin). DFT calculations suggest that this is because PhB(ⁱPr₂Im)₃FeN(Bpin) is unstable toward disproportionation to the starting iron(IV) nitride and PhB(ⁱPr₂Im)₃FeN(Bpin)₂. Attempts at B-C bond insertion using phenyl- and benzyl-pinacol borane were unsuccessful, which we attribute to unfavorable kinetics.

Temperature and energy dependences of ion-molecule reactions: Studies inspired by Diethard Böhme

Diethard Böhme has had a long career covering many topics in ion-molecule reactivity. In this review, we describe the work done at the Air Force Research Laboratory (and its variously named preceding organizations) that was inspired by his studies. These fall into two main areas: nucleophilic displacement (S_N 2) and metal cation chemistry. In S_N 2 chemistry, we revisited many of the reactions Diethard pioneered and studied them in more detail. We found nonstatistical behavior, both competition and noncompetition between multiple channels. New channels were found as hydration occurred, with more solution-like behavior occurring as only a few ligands were added. Temperature-dependent studies revealed details that were not observable at room temperature. These and other highlights will be discussed. In metal cation reactions, Diethard's use of an inductively coupled ion source allowed him to systematically study the periodic table of elements with a number of simple neutrals. We have taken the most interesting of these and studied them in greater detail. In doing so, we were able to identify curve crossing rates, in a few instances information about product states, and the importance of multiple entrance channels.

Effect of Intersystem Crossings on the Kinetics of Thermal Ion-Molecule Reactions: Ti+ + O2, CO2, and N2O

A selected-ion flow tube apparatus has been used to measure rate constants and product branching fractions of ²Ti⁺ reacting with O₂, CO₂, and N₂O over the range of 200-600 K. Ti⁺ + O₂ proceeds at near the Langevin capture rate constant of 6-7 × 10^-10 cm³ s^-1 at all temperatures to yield ⁴TiO⁺ + O. Reactions initiated on doublet or quartet surfaces are formally spin-allowed; however, the 50% of reactions initiated on sextet surfaces must undergo an intersystem crossing (ISC). Statistical theory is used to calculate the energy and angular momentum dependences of the specific rate constants for the competing isomerization and dissociation channels. This acts as an internal clock on the lifetime to ISC, setting an upper limit on the order of τ_ISC < 1e^-11 s. ²Ti⁺ + CO₂ produces ⁴TiO⁺ + CO less efficiently, with a rate constant fit as 5.5 ± 1.3 × 10^-11 (T/300)^{-1.1 ± 0.2} cm³ s^-1. The reaction is formally spin-prohibited, and statistical modeling shows that ISC, not a submerged transition state, is rate-limiting, occurring with a lifetime on the order of 10^-7 s. Ti⁺ + N₂O proceeds at near the capture rate constant. In this case, both Ti⁺ON₂ and Ti⁺N₂O entrance channel complexes are formed and can interconvert over a barrier. The main product is >90% TiO⁺ + N₂, and the remainder is TiN⁺ + NO. Both channels need to undergo ISC to form ground-state products but TiO⁺ can be formed in an excited state exothermically. Therefore, kinetic information is obtained only for the TiN⁺ channel, where ISC occurs with a lifetime on the order of 10^-9 s. Statistical modeling indicates that the dipole-preferred Ti⁺ON₂ complex is formed in ∼80% of collisions, and this value is reproduced using a capture model based on the generic ion-dipole + quadrupole long-range potential.

Intermediate-spin iron(IV)-oxido species with record reactivity

The nonheme iron(IV)-oxido complex trans-N3-[(L¹)Fe^IVO(Cl)]⁺, where L¹ is a derivative of the tetradentate bispidine 2,4-di(pyridine-2-yl)-3,7-diazabicyclo[3.3.1]nonane-1-one, has an S = 1 electronic ground state and is the most reactive nonheme iron model system known so far, of a similar order of reactivity as nonheme iron enzymes (C-H abstraction of cyclohexane, -90 °C (propionitrile), t_1/2 = 3.5 s). The reaction with cyclohexane selectively leads to chlorocyclohexane, but "cage escape" at the [(L¹)Fe^III(OH)(Cl)]⁺/cyclohexyl radical intermediate lowers the productivity. Ligand field theory is used herein to analyze the d-d transitions of [(L¹)Fe^IVO(X)]ⁿ⁺ (X = Cl^-, Br^-, MeCN) in comparison with the thoroughly characterized ferryl complex of tetramethylcyclam (TMC = L²; [(L²)Fe^IVO(MeCN)]²⁺). The ligand field parameters and d-d transition energies are shown to provide important information on the triplet-quintet gap and its correlation with oxidation reactivity.

Measurement of time dependent product branching ratios indicates two-state reactivity in metal mediated chemical reactions

For several decades, the influence of Two State Reactivity (TSR) has been implicated in a host of reactions, but has lacked a stand-alone, definitive experimental kinetic signature identifying its occurrence. Here, we demonstrate that the measurement of a temporally dependent product branching ratio is indicative of spin inversion and is a kinetic signature of TSR. This is caused by products exiting different hypersurfaces with different rates and relative exothermicities. The composite measurement of product intensities with the same mass but with different multiplicities yield biexponential temporal dependences with the sampled product ratio changing in time. These measurements are made using the single photon initiated dissociative rearrangement reaction (SPIDRR) technique which identifies TSR but further determines the kinetic parameters for reaction along the original ground electronic surface in competition with spin inversion and its consequent TSR.

Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory

Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.

This feature article overviews recent work on active spaces, matrix product reference states, treatment of quasidegeneracy, hybrid theory, density-coherence functionals, machine-learned functionals, spin–orbit coupling, gradients, and dipole moments.

Experiment and Theory Clarify: Sc+ Receives One Oxygen Atom from SO2 to Form ScO+, which Proves to be a Catalyst for the Hidden Oxygen-Exchange with SO2

Using Fourier‐transform ion cyclotron resonance mass spectrometry, it was experimentally determined that Sc⁺ in the highly diluted gas phase reacts with SO₂ to form ScO⁺ and SO. By ¹⁸O labeling, ScO⁺ was shown to play the role of a catalyst when further reacting with SO₂ in a Mars‐van Krevelen‐like (MvK) oxygen exchange process, where a solid catalyst actively reacts with the substrate but emerges apparently unchanged at the end of the cycle. High‐level quantum chemical calculations confirmed that the multi‐step process to form ScO⁺ and SO is exoergic and that all intermediates and transition states in between are located energetically below the entrance level. The reaction starts from the triplet surface; although three spin‐crossing points with minimal energy have been identified by computational means, there is no evidence that a two‐state scenario is involved in the course of the reaction, by which the reactants could switch from the triplet to the singlet surface and back. Pivotal to the oxygen exchange reaction of ScO⁺ with SO₂ is the occurrence of a highly symmetric four‐membered cyclic intermediate by which two oxygen atoms become equivalent.

Precise control of the degree and regioselectivity of functionalization in nitro- and amino-functionalized di(trispyrazolylborato)iron(II) spin crossover complexes

Di(trispyrazolylborato)iron(II) ([Tp₂Fe]) complexes represent one of the most robust classes of spin-crossover complexes. Their stability renders them particularly suitable for integration in nanoscale devices, e.g. as sensors or information storage units. While prior studies of the functionalization of those derivatives have been focused on the electronic and steric effects of alkyl and -CF₃ groups in position 3, a pyrazole exchange reaction between nitropyrazole and either trispyrazolylborate or its iron complex allows the regioselective installation of nitro substituents in positions 3, 4 and 5 of the [Tp₂Fe] complexes. The degree of substitution can be varied from 1 to 4 functionalized pyrazoles per complex. The amine-functionalized analogues are accessed by reduction of the nitro analogues under hydrogen transfer conditions. With the exception of di- and tetra-3-NO₂ substituted complexes, all derivatives display spin crossover properties in the solid state, with transition temperatures ranging from 180 to 380 K and showing different degrees of abruptness but no hysteresis. The Slichter-Drickamer model was used to extract the empirical thermodynamic transition parameters, allowing a systematic investigation of the influence of stoichiometry, position, and electronic nature of the substitution on the magnetic properties of the complexes. The steric effects dominate for substitution in position 3 but the electronic effects are significant for the other positions.

Pathways of the Extremely Reactive Iron(IV)-oxido complexes with Tetradentate Bispidine Ligands

The nonheme iron(IV)-oxido complex trans-N3-[(L1 )FeIV =O(Cl)]+ , where L1 is a derivative of the tetradentate bispidine 2,4-di(pyridine-2-yl)-3,7-diazabicyclo[3.3.1]nonane-1-one, is known to have an S=1 electronic ground state and to be an extremely reactive oxidant for oxygen atom transfer (OAT) and hydrogen atom abstraction (HAA) processes. Here we show that, in spite of this ferryl oxidant having the "wrong" spin ground state, it is the most reactive nonheme iron model system known so far and of a similar order of reactivity as nonheme iron enzymes (C-H abstraction of cyclohexane, -90 °C (propionitrile), t1/2 =3.5 sec). Discussed are spectroscopic and kinetic data, supported by a DFT-based theoretical analysis, which indicate that substrate oxidation is significantly faster than self-decay processes due to an intramolecular demethylation pathway and formation of an oxido-bridged diiron(III) intermediate. It is also shown that the iron(III)-chlorido-hydroxido/cyclohexyl radical intermediate, resulting from C-H abstraction, selectively produces chlorocyclohexane in a rebound process. However, the life-time of the intermediate is so long that other reaction channels (known as cage escape) become important, and much of the C-H abstraction therefore is unproductive. In bulk reactions at ambient temperature and at longer time scales, there is formation of significant amounts of oxidation product - selectively of chlorocyclohexane - and it is shown that this originates from oxidation of the oxido-bridged diiron(III) resting state.

Effect of the Ligand Backbone on the Reactivity and Mechanistic Paradigm of Non-Heme Iron(IV)-Oxo during Olefin Epoxidation

The oxygen atom transfer (OAT) reactivity of the non-heme [Fe^IV (2PyN2Q)(O)]²⁺ (2) containing the sterically bulky quinoline-pyridine pentadentate ligand (2PyN2Q) has been thoroughly studied with different olefins. The ferryl-oxo complex 2 shows excellent OAT reactivity during epoxidations. The steric encumbrance and electronic effect of the ligand influence the mechanistic shuttle between OAT pathway I and isomerization pathway II (during the reaction stereo pure olefins), resulting in a mixture of cis-trans epoxide products. In contrast, the sterically less hindered and electronically different [Fe^IV (N4Py)(O)]²⁺ (1) provides only cis-stilbene epoxide. A Hammett study suggests the role of dominant inductive electronic along with minor resonance effect during electron transfer from olefin to 2 in the rate-limiting step. Additionally, a computational study supports the involvement of stepwise pathways during olefin epoxidation. The ferryl bend due to the bulkier ligand incorporation leads to destabilization of both d z 2 and d x 2 - y 2 orbitals, leading to a very small quintet-triplet gap and enhanced reactivity for 2 compared to 1. Thus, the present study unveils the role of steric and electronic effects of the ligand towards mechanistic modification during olefin epoxidation.

Mechanism-Based Inactivation of Cytochrome P450 Enzymes: Computational Insights

Mechanism-based inactivation (MBI) refers to the metabolic bioactivation of a xenobiotic by cytochrome P450s to a highly reactive intermediate which subsequently binds to the enzyme and leads to the quasi-irreversible or irreversible inhibition. Xenobiotics, mainly drugs with specific functional units, are the major sources of MBI. Two possible consequences of MBI by medicinal compounds are drug-drug interaction and severe toxicity that are observed and highlighted by clinical experiments. Today almost all of these latent functional groups (e.g., thiophene, furan, alkylamines, etc.) are known, and their features and mechanisms of action, owing to the vast experimental and theoretical studies, are determined. In the past decade, molecular modeling techniques, mostly density functional theory, have revealed the most feasible mechanism that a drug undergoes by P450 enzymes to generate a highly reactive intermediate. In this review, we provide a comprehensive and detailed picture of computational advances toward the elucidation of the activation mechanisms of various known groups with MBI activity. To this aim, we briefly describe the computational concepts to carry out and analyze the mechanistic investigations, and then, we summarize the studies on compounds with known inhibition activity including thiophene, furan, alkylamines, terminal acetylene, etc. This study can be reference literature for both theoretical and experimental (bio)chemists in several different fields including rational drug design, the process of toxicity prevention, and the discovery of novel inhibitors and catalysts.

Molecular oxofluorides OMFn of nickel, palladium and platinum: oxyl radicals with moderate ligand field inversion

High-valent late transition metal oxo compounds attracted attention because of their peculiar metal–oxygen bond. Their oxo ligands exhibit an electrophilic and distinct radical oxyl (O˙⁻) rather than the more common nucleophilic (O²⁻) character.

Catalytic cycle of the partial oxidation of methane to methanol over Cu-ZSM-5 revealed using DFT calculations

Density functional theory (DFT) calculations were performed to investigate the catalytic cycle of methane conversion to methanol over both [Cu2(O2)]2+ and [Cu2(μ-O)]2+ active sites in the Cu-ZSM-5 catalyst. The [Cu2(O2)]2+ site is found to be active for the partial oxidation of methane to methanol, and although it has a higher energy barrier in the methane activation step, it involves a very low energy barrier in the methanol formation step (36.3 kJ mol-1) as well as a lower methanol desorption energy (52.5 kJ mol-1). As the [Cu2(O2)]2+ active site is also thermodynamically stable, it may play an important role during methane conversion to methanol. Furthermore, the methane activation step follows the homolytic route and the heterolytic route for the [Cu2(O2)]2+ and [Cu2(μ-O)]2+ active sites, respectively, whereas the methanol formation step follows the direct radical rebound mechanism and the indirect rebound mechanism, respectively. Our calculations further indicate that the electronic properties of the reactive O atoms in the active site influence their reactivity toward methane oxidation. More specifically, the higher the spin density and the more negative the charge of the reactive O atom at the active site are, the lower the energy barrier for methane activation will be; and the more negative the charge of the hydroxyl group in the reaction intermediate during the partial oxidation of methane to methanol is, the higher energy barrier of the methanol formation step will be in the triplet state. Furthermore, we used a larger cluster model to predict the mechanism of the methane activation step and the effect of atomic charge of the O atom at the [Cu2(μ-O)]2+ and [Cu2(O2)]2+ active sites on the energy barriers of partial oxidation of methane to methanol, and the conclusions drawn employing the larger cluster model are consistent with those drawn using the smaller double-5T-ring cluster model. In addition, different from the traditional mechanism for methane activation at [Cu2(O2)]2+, which consists of two transition states, we find that the partial oxidation of methane at [Cu2(O2)]2+ can also occur via a single step by direct insertion of one of the O atoms at the active site into the C-H bond of methane.

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.

Experimental and Theoretical Studies on the Mechanism of DDQ-Mediated Oxidative Cyclization of N-Aroylhydrazones

The controversial single-electron-transfer process, frequently proposed in many 2,3-dichloro-5,6-dicyano-1,4-benzoquinone (DDQ)-mediated reactions, was investigated experimentally and theoretically using the oxidative cyclization of aroylhydrazone with DDQ. DDQ-mediated oxadiazole formation involves several processes, including cyclization to form an oxadiazole ring and N-H bond cleavage, either by proton, hydride, or hydrogen atom transfer. The detailed mechanistic study using the M06-2X density functional theory, and the 6-31+G(d,p) basis set, suggests that the pathways involving radical ion pair (RIP) intermediates, which resulted from single-electron transfer (SET), were found to be energetically nearly identical to the pathway without the SET. The substituent-dependent reactivity of oxadiazole formation was consistent with the free energy profiles of both pathways, with or without the SET. This result indicates that in addition to the electron-transfer pathway, the nucleophilic addition/elimination pathway for DDQ should be considered as a possible mechanism of the oxidative transformation reaction using DDQ.

Allenes for Versatile Iron-Catalyzed C-H Activation by Weak O-Coordination: Mechanistic Insights by Kinetics, Intermediate Isolation, and Computation

The iron-catalyzed hydroarylation of allenes was accomplished by weak phenone assistance. The C-H activation proceeded with excellent efficacy and high ortho-regioselectivity in proximity to the weakly coordinating carbonyl group for a range of substituted phenones and allenes. Detailed mechanistic studies, including the isolation of key intermediates, the structural characterization of an iron-metallacycle, and kinetic analysis, allowed the sound elucidation of a plausible catalytic working mode. This mechanistic rationale is supported by detailed computational density functional theory studies, which fully address multi-spin-state reactivity. Furthermore, in operando nuclear magnetic resonance monitoring of the catalytic reaction provided detailed insights into the mode of action of the iron-catalyzed C-H alkylation with allenes.

On the Crucial Role of Isolated Electronic States in the Thermal Reaction of ReC+ with Dihydrogen

Presented here is that isolated, long‐lived electronic states of ReC⁺ serve as the root cause for distinctly different reactivities of this diatomic ion in the thermal activation of dihydrogen. Detailed high‐level quantum chemical calculations support the experimental findings obtained in the highly diluted gas phase using FT‐ICR mass spectrometry. The origin for the existence of these long‐lived excited electronic states and the resulting implications for the varying mechanisms of dihydrogen splitting are addressed.

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.

Spin-Forbidden Reactivity of Transition Metal Oxo Species: Exploring the Potential Energy Surfaces

Spin-forbidden reactions are frequently encountered when transition metal oxo species are involved, particularly in oxygen transfer reactivity. The computational study of such reactions is challenging, because reactants and products are located on different spin potential energy surfaces (PESs). One possible approach to describe these reactions is the so-called minimum energy crossing point (MECP) between the diabatic reactants and products PESs. Alternatively, inclusion of spin-orbit coupling (SOC) effects allows to locate a saddle point on a single adiabatic PES (TS SOC). The TS SOC approach is rarely applied because of its high computational cost. Recently evidence for a TS SOC impact on significantly lowering the activation barrier in dioxygen addition to a carbene-gold(I)-hydride complex reaction (Chem. Sci. 2016, 7, 7034-7039) or even on predicting a qualitatively different reaction mechanism in mercury methylation by cobalt corrinoid (Angew. Chem. Int. Ed. 2016, 55, 11503-11506) has been put forward. Using MECP and TS SOC approaches a systematic analysis is provided here of three prototypical transition metal oxo spin-forbidden processes to investigate their implications on reactivity. Cycloaddition of ethylene to chromyl chloride (CrO₂ Cl₂ +C₂ H₄ ), iron oxide cation insertion into the hydrogen molecule (FeO⁺ +H₂ ) and H-abstraction from toluene by a Mn^V -oxo-porphyrin cation (MnOP(H₂ O)⁺ +C₆ H₅ CH₃ ) are case studies. For all these processes the MECP and TS SOC results are compared, which show that the spin-forbidden reactivity of transition metal oxo species can be safely described by a MECP approach, at least for the first-row transition metals investigated here, where the spin-orbit coupling is relatively weak. However, for the Mn-oxo reactivity, the MECP and TS SOC have been found to be crucial for a correct description of the reaction mechanism. In particular, the TS SOC approach allows to straightforwardly explore detailed features of the adiabatic potential energy surface which in principle could affect the overall reaction rate in cases where the involved diabatic PESs are tricky.