Theoretical study of low-spin, high-spin, and intermediate-spin states of [Fe(III)(pap)2](+) (pap = N-2-pyridylmethylidene-2-hydroxyphenylaminato). Mechanism of light-induced excited spin state trapping

The mechanism of light-induced excited spin state trapping (LIESST) of [FeIII(pap)2]+ (pap = N-2-pyridylmethylidene-2-hydroxyphenylaminato) was discussed on the basis of potential energy surfaces (PESs) of several important spin states, where the PESs were evaluated with the DFT(B3LYP) method. The PES of the quartet spin state crosses those of the doublet and sextet spin states around its minimum. This means that the spin transition occurs from the quartet spin state to either the doublet spin state or the sextet spin state around the PES minimum of the quartet spin state. The PES minimum of the sextet spin state is slightly less stable than that of the doublet spin state by 0.18 eV (4.2 kcal/mol). This small energy difference is favorable for the LIESST. The doublet-sextet spin crossover point is 0.41 eV (9.6 kcal/mol) above the PES minimum of the sextet spin state. Because of this considerably large activation barrier, the thermal spin transition and the tunneling process do not occur easily. In the doublet spin state, the ligand to ligand charge transfer (LLCT) transition is calculated to be 2.16 eV with the TD-DFT(B3LYP) method, in which the pi orbital of the phenoxy moiety and the pi* orbital of the imine moiety in the pap ligand participate. This transition energy is moderately smaller than the visible light of 550 nm used experimentally. In the sextet spin state, the ligand to metal charge transfer (LMCT) transition is calculated to be at 2.36 eV, which is moderately higher than the visible light (550 nm). These results indicate that the irradiation of the visible light induces the LIESST to generate the sextet spin state but the reverse-LIESST is also somewhat induced by the visible light, indicating that the complete spin conversion from the doublet spin state to the sextet one does not occur, as reported experimentally.     

Electronic properties of the low‐lying spin states of dimethylnitrosamine coordinated to Fe(III) heme models: An ab initio study

The unusual O‐coordination mode of nitrosamines to Fe(III) heme models has been observed in the bis(dimethylnitrosamine)(meso‐tetraphenylporphyrinate)iron(III) cation. For the first time, this latter as well as the simpler bis(dimethylnitrosamine)(porphinate)iron(III) heme model cations have been studied through ab initio methods. The sextet, quartet, and doublet spin states of both cations have been studied through single‐point calculations based on the experimental (X‐ray) geometry. Their energies, charges, and spin densities have been analyzed. The obtained results (at the UHF/cc‐pVDZ and ROHF/cc‐pVDZ levels) indicate that the peripheral benzene rings are of secondary importance for the coordination of dimethylnitrosamine to the Fe(III) porphyrin core. The obtained energy ordering is sextet < quartet < doublet, at all computational levels. The UHF, ROHF, and UMP2 results indicate an excess of alpha spin density around the Fe atom, a low covalency for the Fe? O bond and a substantial charge transfer to the Fe atom. Our best estimates [obtained at ROMP2 level with the mixed cc‐pVDZ/cc‐pVTZ‐DK(Fe) basis set] for the energy differences (in eV) between the three spin states considered are 0.929 for the sextet‐quartet gap and 0.812 for the quartet‐doublet gap, which indicate that the spin crossover (at room temperature) is very unlikely. These results represent the substantial decrease in the uncorrelated values. The implications of spin contaminations at the UHF and UMP2 levels for subsequent geometry optimizations to be performed in the smaller cation have also been discussed. © 2013 Wiley Periodicals, Inc.     

Multi-state epoxidation of ethene by cytochrome P450: a quantum chemical study.

The epoxidation of ethene by a model for Compound I of cytochrome P450, studied by the use of density functional B3LYP calculations, involves two-state reactivity (TSR) with multiple electromer species, hence "multi-state epoxidation". The reaction is found to proceed in stepwise and effectively concerted manners. Several reactive states are involved; the reactant is an (oxo)iron(IV) porphyrin cation radical complex with two closely lying spin states (quartet and doublet), both of which react with ethene to form intermediate complexes with a covalent C-O bond and a carbon-centered radical (radical intermediates). The radical intermediates exist in two electromers that differ in the oxidation state of iron; Por(+)(*)Fe(III)OCH(2)CH(2)(*) and PorFe(IV)OCH(2)CH(2)(*) (Por = porphyrin). These radical intermediates exist in both the doublet- and quartet spin states. The quartet spin intermediates have substantial barriers for transformation to the quartet spin PorFe(III)-epoxide complex (2.3 kcal mol(-)(1) for PorFe(IV)OCH(2)CH(2)(*) and 7.2 kcal mol(-)(1) for Por(+)(*)Fe(III)OCH(2)CH(2)(*)). In contrast, the doublet spin radicals collapse to the corresponding PorFe(III)-epoxide complex with virtually no barriers. Consequently, the lifetimes of the radical intermediates are much longer on the quartet- than on the doublet spin surface. The loss of isomeric identity in the epoxide and rearrangements to other products arise therefore mostly, if not only, from the quartet process, while the doublet state epoxidation is effectively concerted (Scheme 7). Experimental trends are discussed in the light of the computed mechanistic scheme, and a comparison is made with closely related mechanistic schemes deduced from experiment.     

The Kinetic Isotope Effect as a Probe of Spin Crossover in the CH Activation of Methane by the FeO+ Cation

Two‐state reactivity (TSR) is often used to explain the reaction of transition‐metal–oxo reagents in the bare form or in the complex form. The evidence of the TSR model typically comes from quantum‐mechanical calculations for energy profiles with a spin crossover in the rate‐limiting step. To prove the TSR concept, kinetic profiles for C? H activation by the FeO⁺ cation were explored. A direct dynamics approach was used to generate potential energy surfaces of the sextet and quartet H‐transfers and rate constants and kinetic isotope effects (KIEs) were calculated using variational transition‐state theory including multidimensional tunneling. The minimum energy crossing point with very large spin–orbit coupling matrix element was very close to the intrinsic reaction paths of both sextet and quartet H‐transfers. Excellent agreement with experiments were obtained when the sextet reactant and quartet transition state were used with a spin crossover, which strongly support the TSR model.     

How molecular oxygen binds to bis[trifluoroacetylacetonato(-1)]cobalt(II)--ab initio and density functional theory studies

Cobalt(II) diketonate complexes, such as bis[trifluoroacetylacetonato(-1)]cobalt(II) [Co(tfa)(2)], catalyze the aerobic oxidation of alkenols into functionalized tetrahydrofurans. To gain insight into activation of triplet dioxygen by Co(tfa)(2) in a protic solvent, as used in oxidation catalysis, the electronic structure of aquabis[trifluoroacetylacetonato(-1)]cobalt(II)--Co(tfa)(2)(H(2)O)--and the derived dioxygen adduct were characterized using ab initio (CASSCF, NEVPT2) and density functional theory (BP86, TPSSh, B3LYP) methods. The ground state of Co(tfa)(2)(H(2)O) is a high-spin, quartet state. As dioxygen approaches the cobalt atom, the quartet state couples with a triplet dioxygen molecule and forms a sextet, a quartet, and a doublet spin state with the high-spin state being the lowest in energy. At the equilibrium Co-O(2) distance of 1.9 ?, Co(tfa)(2)(H(2)O)(O(2)) has a doublet superoxo Co(III) ground state with the unpaired electron residing on the oxygen moiety, in a nearly unchanged O(2)π* orbital.     

The Kinetic Isotope Effect as a Probe of Spin Crossover in the C?H Activation of Methane by the FeO+ Cation

Two‐state reactivity (TSR) is often used to explain the reaction of transition‐metal–oxo reagents in the bare form or in the complex form. The evidence of the TSR model typically comes from quantum‐mechanical calculations for energy profiles with a spin crossover in the rate‐limiting step. To prove the TSR concept, kinetic profiles for C H activation by the FeO⁺ cation were explored. A direct dynamics approach was used to generate potential energy surfaces of the sextet and quartet H‐transfers and rate constants and kinetic isotope effects (KIEs) were calculated using variational transition‐state theory including multidimensional tunneling. The minimum energy crossing point with very large spin–orbit coupling matrix element was very close to the intrinsic reaction paths of both sextet and quartet H‐transfers. Excellent agreement with experiments were obtained when the sextet reactant and quartet transition state were used with a spin crossover, which strongly support the TSR model.     

High-level ab initio studies of unimolecular dissociation of the ground-state N3 radical

A comprehensive study of the unimolecular dissociation of the N(3) radical on the ground doublet and excited quartet potential energy surfaces has been carried out with multireference single and double excitation configuration interaction and second-order multireference perturbation methods. Two forms of the N(3) radical have been located in the linear and cyclic region of the lowest doublet potential energy surface with an isomerization barrier of 62.2 kcal/mol above the linear N(3). Three equivalent C(2v) minima of cyclic N(3) are connected by low barrier, meaning the molecule is free to undergo pseudorotation. The cyclic N(3) is metastable with respect to ground state products, N((4)S)+N(2), and dissociation must occur via intersystem crossing to a quartet potential energy surface. Minima on the seams of crossing between the doublet and quartet potential surfaces are found to lie substantially higher in energy than the cyclic N(3) minima. This strongly suggests that cyclic N(3) possesses a long collision-free lifetime even if formed with substantial internal excitation.     

Effect of axial ligand on the electronic configuration, spin states, and reactivity of iron oxophlorin

Iron-oxophlorin is an intermediate in heme degradation, and the nature of the axial ligand can alter the spin, electron distribution, and reactivity of the metal and the oxophlorin ring. The structure and reactivity of iron-oxophlorin in the presence of imidazole, pyridine, and t-butyl isocyanide as axial ligands was investigated using the B3LYP and OPBE methods with the 6-31+G* and 6-311+G** basis sets. OPBE/6-311+G** has shown that the doublet state of [(Py)(2)Fe(III)(PO)] (where pyridines are in perpendicular planes and PO is the oxophlorin trianion) is 3.45 and 5.27 kcal/mol more stable than the quartet and sextet states, respectively. The ground-state electronic configuration of the aforementioned complex is π(xz)(2) π(yz)(2) a(2u)(2) d(xy)(1) at low temperatures and changes to π(xz)(2) π(yz)(2) d(xy)(2) a(2u)(1) at high temperatures. This latter electronic configuration is consistently seen for the [(t-BuNC)(2)Fe(II)(PO(?))] complex (where PO(?) is the oxophlorin dianion radical). The complex [(Im)(2)Fe(III)(PO)] adopted the d(xy)(2) (π(xz) π(yz))(3) ground state and has low-lying quartet excited state which is readily populated when the temperature is increased.     

DFT Study on the Os^＋-catalytic Reduction Reaction of N2O with H2 in the Gas Phase

The Os+-catalytic reduction of N2O by H2 in gas phase has been theoretically investigated with B3LYP method.The reaction mechanisms on the sextet and quartet surfaces were found to be similar.The calculated sextet potential energy profiles show that the two reactions involved in the catalytic cycle,Os+ + N2O → OsO+ + N2 and OsO+ + H2 → Os+ + H2O,have barriers of 28.3 and 123.3 kJ/mol,respectively.In contrast,the reactions on the quartet surfaces are energetically much more favorable.These results rationalize the experimentally observed low catalytic reactivity of sextet(ground-state) Os+.Further,the crossing between the sextet and quartet surfaces are also suggested and qualitatively discussed.     

Two states and two more in the mechanisms of hydroxylation and epoxidation by cytochrome P450

Past studies have shown that oxidation reactions by P450 Compound I (Cpd I) can be described by two competing quartet and doublet spin states, which possess three unpaired electrons, hence tri-radicals. One electron excitation from the delta orbital to sigma* xy generates two states that possess five unpaired electrons, so-called penta-radicals, in sextet and quartet situations, and which were shown by theory to lie only approximately 12-14 kcal/mol higher in energy than the tri-radical ground states (ref 7). The present study focuses on the C-H hydroxylation and C=C epoxidation of propene by these penta-radical states. It is shown that the initial energy differences, between the penta-radical and tri-radical states, diminish along the reaction pathway, due to the favorable and cumulative exchange stabilization of the more open-shell species. Furthermore, theory suggests that hydrogen bonding to the thiolate ligand, and general polarity of the environment, reduce these gaps further, thereby making the penta-radical states accessible to ground-state reactivity. The interconversion between the tri-radical and penta-radical states along the reaction coordinate will depend on the dynamics of spin-flips and energy barriers between the states. Especially interesting should be the region of the reaction intermediates; for both epoxidation and hydroxylation, this region is typified by a dense manifold of spin states and electromeric states (that differ by the oxidation state of iron), such that the total reactivity would be expected to reflect the interplay of these states, giving rise to multistate reactivity.     

DFT Study on the Os+-catalytic Reduction Reaction of N2O with He in the Gas Phase

The Os+-catalytic reduction of N<,2>O by H<,2> in gas phase has been theoretically investigated with B3LYP method.The reaction mechanisms on the sextet and quartet surfaces were found to be similar.The calculated sextet potential energy profiles show that the two reactions involved in the catalytic cycle, Os<'+>+ N<,2>O→OsO<'+>+N<,2> and OsO<'+>+H2→Os<'+>+H<,2>O, have barriers of 28.3 and 123.3 kJ/mol, respectively.In contrast, the reactions on the quartet surfaces are energetically much more favorable.These results rationalize the experimentally observed low catalytic reactivity of sextet (ground-state) Os<'+>.Further, the crossing between the sextet and quartet Surfaces are also suggested and qualitatively discussed.     

Substitution effects on olefin epoxidation catalyzed by Oxoiron(IV) porphyrin π-cation radical complexes: A dft study

The effects of peripheral fluorine atoms on epoxidation reactions of ethylene by oxoiron(IV) porphyrin cation radical complex in the quartet and sextet spin multiplicities are systematically investigated using the DFT method. The overall reaction routes are determined using a model system of ethylene and Fe(IV)OCl-porphyrin with substituted fluorine atoms. By obtaining the energy diagrams and electron- and spin-density difference contour maps of the transition states and intermediate compounds, we confirm that the electron-withdrawing by peripheral fluorine atoms enhances the reactivity as the number of fluorine atoms increases, as is observed experimentally. The intersystem crossing between the quartet and sextet spin multiplicities is discussed by means of the intrinsic reaction coordinate method. We conclude that the rate-determining step is located at the first transition state (TS1) for the activation of CC and FeO bonds, and the ground electronic state changes from quartet to sextet around the TS1. © 2019 Wiley Periodicals, Inc.     

气相VO（∑＋）与CH3OH（1^A′）分子自旋禁阻反应C—H，O—H键活化机理的理论研究

采用密度泛函理论（DFT）B3LYP与cCsD方法研究了二重态和四重态势能面自旋禁阻反应VO（∑’）活化cH30H（1^A′）分子c—H,0—H键的微观机理．通过自旋一轨道耦合的计算讨论了势能面交叉点和可能的自旋翻转过程．在MEcP处,四重态和二重态问的旋轨耦合常数为131．14cm^-1．自旋多重度发生改变,从四重态系间穿越到二重态势能面形成中间体2^IM1,导致反应势能面的势垒明显降低．     

Ti^＋离子和C2H4分子自旋禁阻反应中C—H键活化机理的理论研究

用密度泛函B3LYP方法,研究了二重态和四重态势能面自旋禁阻反应Ti^＋（^4F,3d^24s^1）＋C2H4→TiC2H2^＋（^2A2）＋H2的微观机理．通过自旋．轨道耦合的计算讨论了势能面交叉点和可能的自旋翻转过程．中间体IM1-^4B2处,四重态和二重态间的旋-轨耦合值为59．3cm^-1．自旋多重度必将发生变化,从四重态系间穿越到二重态势能面形成共价型复合物IM1-^2A1,同时导致四重态势能面的势垒明显降低．到插入中间体IM2后,二重态势能面上有两条不同的反应路径,即分步和协同路径,后者在二重态势能面上得到放热产物TiC2H2^＋（^2A2）＋H2具有较低的活化势垒,4．52kcal／mol,其主反应路径为：Ti^＋＋C2H4→^4IC→IM1—^4B2→4.2ISC→IM1—^2A1→[^2TSins]→IM2-^2A”→[^2TSMCTS]→IM5→TiC2H2^＋（^2A2）＋H2．     

Cobalt-centered ten-vertex germanium clusters: the pentagonal prism as an alternative to polyhedra predicted by the Wade-Mingos rules

One of the most exciting recent (2009) discoveries in metal cluster chemistry is the pentagonal prismatic Co@Ge(10)(3-) ion, found in [K(2,2,2-crypt)](4)[Co@Ge(10)][Co(1,5-C(8)H(12))(2)]·toluene and characterized structurally by X-ray diffraction. The complete absence of triangular faces in the pentagonal prismatic structure of Co@Ge(10)(3-) contradicts expectations from the well-established Wade-Mingos rules, which predict polyhedral structures having mainly or entirely triangular faces. A theoretical study on Co@Ge(10)(z) systems (z = -5 to +1) predicts a singlet D(5h) pentagonal prismatic global minimum for the trianion Co@Ge(10)(3-) in accord with this experimental result. Redox reactions on this pentagonal prismatic Co@Ge(10)(3-) trianion generate low-energy pentagonal prismatic structures for Co@Ge(10)(z) where z = 0, -1, -2, -4, and -5 having quartet, triplet, doublet, doublet, and triplet spin states, respectively. Similar theoretical methods predict a singlet C(3v) polyhedral structure for the monoanion Co@Ge(10)(-), similar to previous theoretical predictions on the isoelectronic neutral Ni@Ge(10) and the structure realized experimentally in the isoelectronic Ni@In(10)(10-) found in the K(10)In(10)Ni intermetallic. Redox reactions on this C(3v) polyhedral Co@Ge(10)(-) monoanion generate low energy C(3v) polyhedral structures for Co@Ge(10)(z) where z = 0, -2, -3, and -4 having doublet, doublet, triplet, and quartet spin states, respectively.     

Theoretical study of the reaction of Ni with OCS

The reactivity of Ni⁺ with OCS on both doublet and quartet potential energy surfaces (PES) has been investigated at the B3LYP/6-311+G(d) level. The object of this investigation was the elucidation of the reaction mechanism. The calculated results indicated that both the CS and CO bond activations proceed via an insertion–elimination mechanism. Intersystem crossing between the doublet and quartet surfaces may occur along both the CS and CO bond activation branches. The ground states of NiS⁺ and NiO⁺ were found to be quartets, whereas NiCO⁺ and NiCS⁺ have doublet ground states. The CS bond activation is energetically much more favorable than the CO bond activation. All theoretical results are in line with early experiments.     

Theoretical Investigation for Two-state Reactivity of CO_2 Hydrogenation Catalyzed by Ru in Gas Phase

Gas-phase CO_2 catalyzed activation hydrogenation by Ru atoms was studied with density functional theory. Based on the structure optimization of different potential energy surfaces,there are two crossing points between singlet and triplet potential energy surfaces and there is a crossing point between quintet and triplet potential energy surfaces in the whole catalytic cycle. Spin transition probabilities in the vicinity of the intersections have been calculated by the Landau-Zener model theory. There are three minimum energy crossing points(MECPs) with strong spin-orbital coupling effect and higher spin transition probability,and all spin inversion occurred in s orbital and different d orbitals of ruthenium,indicating this is a typical two-state reactivity(TSR) reaction. Finally,the lowest energy reaction path is ensured.     

OsO~+氧化活化氢分子气相反应机理的密度泛函理论计算

为了探寻OsO+与H2气相反应的机理,用密度泛函理论方法 UB3LYP,全优化了该反应的加成(氧化加成和[2+2]环加成)-消除、氢抽提-反弹,以及氧端插入等四种可能路径中所有可能的反应物、中间体、过渡态和产物在六重态、四重态和二重态等三个自旋态下的几何结构,计算了各种机理反应的势能面.结果表明,标题反应为自旋禁阻反应,反应起始自四重态,最终产物为六重态基态,整个反应放热21.0 kJ mol-1.因反应络合物相对于入口通道有太正Gibbs函数,氧端插入机理是高能的过程.其他三种机理都具有多(或二)态反应性(MSR或TSR).其中,两种加成-消除机理的最低能量路径都可能经由四重态-二重态-四重态-六重态的三次自旋翻转,抽提-反弹机理的最低能量路径可能经历由四重态-六重态的自旋翻转.抽提-反弹机理由势能面一路攀升的吸热氢抽提过程和几乎无能垒的强放热的反弹过程组成,所以按该机理反应在常温常压下难以发生.两种加成-消去机理的决速步(第二个H的迁移步)相同,虽然其位垒稍高,为156.9 kJ mol-1,但与其进程中前面的强放热步骤耦合,常温常压下该反应是可以发生的.其中,协同环加成步的位垒仅28.7 kJ mol-1,比第一个H的还原消去步的位垒低113.7 kJ mol-1,所以竞争的结果是,常温常压下[2+2]环加成-消去机理比氧化加成-消去机理在动力学上更有利.     

Dissociative addition of water to a main group 13 metal cluster: computational study of the reaction of gallium dimer with H2O

We have investigated the lowest triplet and singlet potential energy surfaces (PESs) for the reaction of Ga(2) dimer with water. Under thermal conditions, we predict formation of the triplet ground state addition complex Ga(2)···OH(2)((3)B(1)) involving Ga···O···Ga bridge interaction. At the coupled cluster CCSD(T)/AE (CCSD(T)/ECP) computational levels, Ga(2)···OH(2)((3)B(1)) is bound by 5.5 (5.7) kcal/mol with respect to the ground state reactants Ga(2)((3)Π(u)) + H(2)O. Identification of the addition complex is in agreement with the experimental evidence from matrix isolation infrared (IR) spectroscopy reported recently by Macrae and Downs. The located minimum energy crossing points (MECPs) between the triplet and singlet energy surfaces on the entrance channel of Ga(2) + H(2)O are not expected to be energetically accessible under the matrix conditions, consistent with the lack of occurrence of Ga(2) insertion into the O-H bond under such conditions. The computed energies and harmonic and anharmonic vibrational frequencies for the triplet and singlet Ga(2)(H)(OH) insertion isomers indicate the singlet double-bridged Ga(μ-H)(μ-OH)Ga isomer to be the most stable and support the experimental IR identification of this species. The energy barrier for elimination of H(2) from the second most stable singlet HGa(μ-OH)Ga insertion isomer found to be 13.9 (12.9) kcal/mol is also consistent with the available experimental data.     

Electron spin multiplicities of transition-metal aromatic radicals and ions: M[C6(CH3)6] and M(+)[(C6(CH3)6] (M = Ti, V, and Co)

Determination of electron spin multiplicities of transition-metal radicals and ions challenges both experimentalists and theoreticians. In this work, we report preferred electron spin states of M[C(6)(CH(3))(6)] and M(+)(C(6)(CH(3))(6)], where M = Ti, V, and Co. The neutral radicals were formed in a supersonic metal cluster beam source, and their masses were measured with time-of-flight mass spectrometry. Precise ionization energies of the radicals and metal-ligand stretching frequencies of the ions were measured by pulsed field ionization zero electron kinetic energy spectroscopy. C-H stretching frequencies of the methyl group in the radicals were obtained by infrared-ultraviolet two-photon ionization. Electron spin multiplicities of the radicals and ions were investigated by combining the spectroscopic measurements, density functional theory, and Franck-Condon factor calculations. The preferred spin states are quintet, sextet, and quartet for the neutral Ti, V, and Co radicals, respectively; for the corresponding singly charged cations, they are quartet, quintet, and triplet. In these high-spin states, the aromatic ring remains nearly planar. This finding contrasts to the previous study of Sc(hmbz), for which low-spin states are favored, and the aromatic ring is severely bent.