Cover Image

Spin-inversion dynamics in O₂ binding to a model heme complex, which consisted of Fe(II)-porphyrin and imidazole, were studied using nonadiabatic wave packet dynamics calculations. We considered three active nuclear degrees of freedom in the dynamics, including the motions along the Fe–O distance, Fe–O–O angle, and Fe out-of-plane distance. Spin-free potential energy surfaces for the singlet, triplet, quintet, and septet states were developed using density functional theory calculations, and spin–orbit coupling elements were obtained from CASSCF-level electronic structure calculations. The spin-inversion mainly occurred between the singlet state and one of the triplet states due to large spin–orbit couplings and the contributions of other states were extremely small. The present quantum dynamics calculations suggested that the narrow crossing region model plays a dominant role in the O₂ binding dynamics. In addition, the one-dimensional Landau–Zener model underestimated the nonadiabatic transition probability.     

Electronic structure and spin bonds in the [(NH3)5Cr-O-Cr(NH3)5]4+ complex

The electronic structure of the lowest singlet, triplet, quintet, and septet states of the binuclear complex [Cr(NH₃)₅]₂O⁴⁺ are calculated by the semiempirical MC SCF method in the INDO approximation. Schemes of spin coupling of Cr-O-Cr in these states are analyzed. Translated fromZhumal Struktumoi Khimii, Vol. 38, No. 3, pp. 457–464, May–June, 1997.     

Mo活化NH3的自旋禁阻反应机理

采用密度泛函理论的UB3LYP方法,计算研究了气相中Mo活化NH3的反应机理。为了理解由Mo活化NH3过程中自旋翻转行为,对自旋态分别为七、五、三和单等4个反应势能面进行了计算研究,其结果表明,Mo活化NH3的过程是通过各自旋态势能面交叉产生的典型的自旋禁阻反应,最低能量交叉点(MECPs)附近的系间窜越导致2步H转移和脱H2反应能垒降低。此外运用自然键(NBO)轨道理论分析了反应中较为重要的几个物种的成键特性。通过计算在最低能量交叉点(MECPs)附近不同自旋态之间的自旋-轨道耦合常数,再运用Landau-Zener跃迁几率公式估算了MECPs处系间窜越几率。所确定的最低能量反应路径为:7Mo+NH3→7IM1→7/5MECP1→5TS12→5IM2→5/3MECP2→3TS23→3IM3→3TS34→3IM4→3HMoN+H2。     

What factors influence the ratio of C-H hydroxylation versus C=C epoxidation by a nonheme cytochrome P450 biomimetic?

Density functional calculations on a nonheme biomimetic (Fe=O(TMCS+) have been performed and its catalytic properties versus propene investigated. Our studies show that this catalyst is able to chemoselectively hydroxylate C=H bonds even in the presence of C=C double bonds. This phenomenon has been analyzed and found to occur due to Pauli repusions between protons on the TMCS ligand with protons attached to the approaching substrate. The geometries of the rate determining transition states indicate that the steric hindrance is larger in the epoxidation transition states than in the hydroxylation ones with much shorter distances; hence the hydroxylation pathway is favored over the epoxidation. Although, the reactant experiences close lying triplet and quintet spin states, the dominant reaction mechanism takes place on the quintet spin state surface; i.e., Fe=O(TMCS)+ reacts via single-state reactivity. Our calculations show that this spin state selectivity is the result of geometric orientation of the transition state structures, whereby the triplet ones are destabilized by electrostatic repulsions between the substrate and the ligand while the quintet spin transition states are aligned along the ideal axis. The reactivity patterns and geometries are compared with oxoiron species of dioxygenase and monoxygenase enzymes. Thus, Fe=O(TMCS)+ shows some similarities with P450 enzyme reactivity: it chemoselectively hydroxylates C=H bonds even in the presence of a C=C double bond and therefore is an acceptable P450 biomimetic. However, the absolute barriers of substrate oxidation by Fe=O(TMCS)+ are higher than the ones obtained with heme enzymes, but the chemoselectivity is lesser affected by external perturbations such as hydrogen bonding of a methanol molecule toward the thiolate sulfur or a dielectric constant. This is the first oxoiron complex whereby we calculated a chemoselective hydroxylation over epoxidation in the gas phase.     

EPR Studies on Branched High‐Spin Arylnitrenes

The UV (λ>305 nm) photolysis of triazide 3 in 2‐methyl‐tetrahydrofuran glass at 7 K selectively produces triplet mononitrene 4 (g=2.003, D_T=0.92 cm^?1, E_T=0 cm^?1), quintet dinitrene 6 (g=2.003, D_Q=0.204 cm^?1, E_Q=0.035 cm^?1), and septet trinitrene 8 (g=2.003, D_S=?0.0904 cm^?1, E_S=?0.0102 cm^?1). After 45 min of irradiation, the major products are dinitrene 6 and trinitrene 8 in a ratio of ～1:2, respectively. These nitrenes are formed as mixtures of rotational isomers each of which has slightly different magnetic parameters D and E. The best agreement between the line‐shape spectral simulations and the experimental electron paramagnetic resonance (EPR) spectrum is obtained with the line‐broadening parameters Γ(E_Q)=180 MHz for dinitrene 6 and Γ(E_S)=330 MHz for trinitrene 8 . According to these line‐broadening parameters, the variations of the angles Θ in rotational isomers of 6 and 8 are expected to be about ±1 and ±3°, respectively. Theoretical estimations of the magnetic parameters obtained from PBE/DZ(COSMO)//UB3LYP/6‐311+G(d,p) calculations overestimate the E and D values by 1 and 8 %, respectively. Despite the large distances between the nitrene units and the extended π systems, the zero field splitting (zfs) parameters D are found to be close to those in quintet dinitrenes and septet trinitrenes, where the nitrene centers are attached to the same aryl ring. The large D values of branched septet nitrenes are due to strong negative one‐center spin–spin interactions in combination with weak positive two‐center spin–spin interactions, as predicted by theoretical considerations.     

Ultrafast Deactivation Mechanism of the Excited Singlet in the Light‐Induced Spin Crossover of [Fe(2,2′‐bipyridine)3]2+

The mechanism of the light‐induced spin crossover of the [Fe(bpy)₃]²⁺ complex (bpy=2,2′‐bipyridine) has been studied by combining accurate electronic‐structure calculations and time‐dependent approaches to calculate intersystem‐crossing rates. We investigate how the initially excited metal‐to‐ligand charge transfer (MLCT) singlet state deactivates to the final metastable high‐spin state. Although ultrafast X‐ray free‐electron spectroscopy has established that the total timescale of this process is on the order of a few tenths of a picosecond, the details of the mechanisms still remain unclear. We determine all the intermediate electronic states along the pathway from low spin to high spin and give estimates for the deactivation times of the different stages. The calculations result in a total deactivation time on the same order of magnitude as the experimentally determined rate and indicate that the complex can reach the final high‐spin state by means of different deactivation channels. The optically populated excited singlet state rapidly decays to a triplet state with an Fe d⁶(${{\rm t}{{5\hfill \atop {\rm 2g}\hfill}}}$ ${{\rm e}{{1\hfill \atop {\rm g}\hfill}}}$ ) configuration either directly or by means of a triplet MLCT state. This triplet ligand‐field state could in principle decay directly to the final quintet state, but a much faster channel is provided by internal conversion to a lower‐lying triplet state and subsequent intersystem crossing to the high‐spin state. The deactivation rate to the low‐spin ground state is much smaller, which is in line with the large quantum yield reported for the process.     

The powder X-band electron paramagnetic resonance spectroscopy of septet pyridyl-2,4,6-trinitrene

The first X-band EPR spectrum containing only non-overlapping signals of septet pyridyl-2,4,6-trinitrene and triplet pyridylnitrenes is reported. This spectrum was recorded after photolysis of 2,4,6-triazidopyridine in solid argon at 5 K. The zero-field splitting (ZFS) parameters of this trinitrene as well as of intermediate triplet mononitrenes and quintet dinitrenes formed at early stages of the photolysis were determined using the combination of modern computer line-shape spectral simulations and density functional theory (DFT) calculations. It was found that septet pyridyl-2,4,6-trinitrene has the record negative parameter D_S = −0.1031 cm⁻¹ among all known to date septet pyridyl-2,4,6-trinitrenes and may be of interest as a model multi-qubit spin system for investigations of quantum computation processing.     

UNO DMRG CAS CI calculations of binuclear manganese complex Mn(IV)2O2(NHCHCO2)4: Scope and applicability of Heisenberg model

Both direct exchange and super-exchange interactions cooperate to realize inter-spin magnetic interaction in binuclear manganese complex Mn(IV)₂O₂(NHCHCO₂)₄ with a di-μ-oxo path. We revisited this spin system using DMRG CAS methods and CAS selection procedures. Our results indicate that our previous “dynamically extended spin polarization” (DE-SP) procedure for organic polyradicals and so forth does not work well. Thus, we have examined another selection procedure, the “dynamically extended super-exchange” (DE-SE) procedure. DMRG CASCI [18,18] by UB3LYP(HS)-UNO(DE-SE) can realize antiferromagnetic J values similar to experimental ones (−87 cm⁻¹). In addition, all J values between all spin states (HS[septet],IS[quintet],IS[triplet],LS[singlet])were also shown to be correct under sufficiently large M values. © 2018 Wiley Periodicals, Inc.     

Global ab initio potential energy surfaces for the O2(3Σg-)+O2(3Σg-) interaction

Completely ab initio global potential energy surfaces (PESs) for the singlet and triplet spin multiplicities of rigid O(2)((3)Σ(g)(-))+O(2)((3)Σ(g)(-)) are reported for the first time. They have been obtained by combining an accurate restricted coupled cluster theory with singles, doubles, and perturbative triple excitations [RCCSD(T)] quintet potential [Bartolomei et al., J. Chem. Phys. 128, 214304 (2008)] with complete active space second order perturbation theory (CASPT2) or, alternatively, multireference configuration interaction (MRCI) calculations of the singlet-quintet and triplet-quintet splittings. Spherical harmonic expansions, containing a large number of terms due to the high anisotropy of the interaction, have been built from the ab initio data. The radial coefficients of these expansions are matched at long range distances with analytical functions based on recent ab initio calculations of the electric properties of the monomers [M. Bartolomei, E. Carmona-Novillo, M. I. Hernández, J. Campos-Martínez, and R. Hernández-Lamoneda, J. Comput. Chem. (2010) (in press)]. The singlet and triplet PESs obtained from either RCCSD(T)-CASPT2 or RCCSD(T)-MRCI calculations are quite similar, although quantitative differences appear in specific terms of the expansion. CASPT2 calculations are the ones giving rise to larger splittings and more attractive interactions, particularly in the region of the absolute minima (in the rectangular D(2h) geometry). The new singlet, triplet, and quintet PESs are tested against second virial coefficient B(T) data and, their spherically averaged components, against integral cross sections measured with rotationally hot effusive beams. Both types of multiconfigurational approaches provide quite similar results, which, in turn, are in good agreement with the measurements. It is found that discrepancies with the experiments could be removed if the PESs were slightly more attractive. In this regard, the most attractive RCCSD(T)-CASPT2 PESs perform slightly better than the RCCSD(T)-MRCI counterpart.     

Study of the Light‐Induced Spin Crossover Process of the [FeII(bpy)3]2+ Complex

Ab initio calculations have been performed on [Fe^II(bpy)₃]²⁺ (bpy=bipyridine) to establish the variation of the energy of the electronic states relevant to light‐induced excited‐state spin trapping as a function of the Fe? ligand distance. Light‐induced spin crossover takes place after excitation into the singlet metal‐to‐ligand charge‐transfer (MLCT) band. We found that the corresponding electronic states have their energy minimum in the same region as the low‐spin (LS) state and that the energy dependence of the triplet MLCT states are nearly identical to the ¹MLCT states. The high‐spin (HS) state is found to cross the MLCT band near the equilibrium geometry of the MLCT states. These findings give additional support to the hypothesis of a fast singlet–triplet interconversion in the MLCT manifold, followed by a ³MLCT–HS (⁵T₂) conversion accompanied by an elongation of the Fe? N distance.     

On the catalytic role of Re+ in the oxygen transport activation of N2O by CO

The mechanism of the cyclic reaction N₂O(¹∑) + CO(¹∑⁺) → N₂(¹∑ _g ⁺ ) + CO₂(¹∑ _g ⁺ ) catalyzed by Re⁺ has been investigated on quintet and septet potential energy surfaces (PES). The reactions were studied by the B3LYP density functional method and the CCSD(T) theory. The calculated results of different PES show that the reaction proceeds in a two-step manner and spin crossing between different PES occurs. The involving crossing points (CPs) between the quintet and septet PES have been discussed by means of the intrinsic reaction coordinate approach. And the O-atom affinities testified that Re⁺ can capture O from N₂O and transfer O atom to CO in the two spin state, which are thermodynamically allowed. Furthermore, the spin–orbit coupling (SOC) is calculated between electronic states of different multiplicities at the CPs. For CP1 and CP2, the computed SOC constants are 8.34 and 10.09 cm^?1, respectively, obtained by using one-electron spin–orbit Hamiltonian in GAMESS. Therefore, the intersystem crossing at CP1 and CP2 occurs with a little probability because of the small SOC involved.     

Reaction mechanism of NO with hydrolysates of NAMI-A: an MD simulation by combining the QM/MM(ABEEM) with the MD-FEP method

Nitrosylation reaction mechanisms of the hydrolysates of NAMI-A and hydrolysis reactions of ruthenium nitrosyl complexes were investigated in the triplet state and the singlet state. Activation free energies were calculated by combining the QM/MM(ABEEM) method with free energy perturbation theory, and the explicit solvent environment was simulated by an ABEEMσπ polarizable force field. Our results demonstrate that nitrosylation reactions of the hydrolysates of NAMI-A occur in both the triplet and the singlet states. The Ru-N-O angle of the triplet ruthenium nitrosyl complexes is in the range of 132.0°–138.2°. However, all the ruthenium nitrosyl complexes at the singlet state show an almost linear Ru-N-O angle. The nitrosylation reaction happens prior to the hydrolysis reaction for the first-step hydrolysates. The activation free energies of the nitrosylation reactions show that the H₂O-NO exchange reaction of [RuCl₄(Im)(H₂O)] in the singlet spin sate is the most likely one. Comparing with the activation free energies of the hydrolysis reactions of the ruthenium nitrosyl complexes, the results indicate that the rate of the DMSO–H₂O exchange reaction of [RuCl₃(NO)(Im)(DMSO)] is faster than that of [RuCl₃(H₂O)(Im)(DMSO)] in both the triplet spin state and the singlet spin state. © 2018 Wiley Periodicals, Inc.     

Nb+离子活化甲烷脱氢反应机理密度泛函(DFT)研究

通过DFT-UB3LYP方法, 计算了五重、三重和单重自旋态下的气相Nb+离子活化甲烷脱氢反应的能量变化, 并对其直接式和插入式反应机理进行了比较, 考察了自旋翻转对反应的影响. 结果表明, 插入式脱氢较直接式有利, CH4上的H转移到Nb+上形成的中间体HNbCH+3中, 多重度由五重降为三重, 反应活化能垒显著降低; HNbCH+3可经四中心过渡态转化为(H2)NbCH+2, 最后生成三重态的NbCH+2+H2. 速控步骤为(H2)NbCH+2的脱氢. 此外, 通过对V+, Nb+, Ta+活化甲烷的比较研究了三者活化甲烷的反应活性.     

基于非血红素铁模型配合物[FeⅣ(O)(N4Py)]2+的新型N杂环卡宾配合物的理论设计

基于四价非血红素铁模型配合物[FeⅣ(O)(N4Py)]2+, 通过理论计算设计出一种新型N杂环卡宾配合物[FeⅣ(O)(N4Py)]2+. 采用密度泛函理论B3LYP方法, 计算了[FeⅣ(O)(N4Py)]2+的几何结构和电子结构, 并研究了[FeⅣ(O)(N4Py)]2+使环己烷C-H键羟基化的反应机理. 结果表明, [FeⅣ(O)(N4Py)]2+的五重态能量比基态三重态能量高约5.7 kJ/mol, 故五重态几乎不能参与反应. 赤道方向的配位基N杂环卡宾(NHC)对FeO单元的σ-贡献要大于N4Py的贡献, 而它的空间位阻效应也大于N4Py, 因此2+的稳定性强于[FeⅣ(O)(N4Py)]2+. [FeⅣ(O)(N4Py)]2+的三重态的反应能垒比[FeⅣ(O)(N4Py)]2+的三重态反应能垒高2.0 kJ/mol, 且为单态反应, 所以[FeⅣ(O)(N4Py)]2+的反应活性要高于[FeⅣ(O)(N4Py)]2+.     

草酰胺衍生物桥联双核镍配合物从头算研究

采用LanL2DZ基组 ,对N,N′ 双(2 苯甲酸根)草酰胺桥联双核镍配合物Ni(obbz)Ni(H2O)4 进行从头计算研究 ,探讨该配合物单、三重态的电子组态的稳定性、电子结构特征及电子自旋布居规律等。计算结果表明 ,该配合物分子的三重态比单重态稳定 ,因此 ,该配合物择型于三重态的电子组态。电子自旋主要布居于八面体环境的Ni(1)中心上 ,而处于四方配位环境的Ni(2)中心则没有发现电子自旋布居。同时 ,Ni(2)中心主要参与的分子轨道都处于低能区 ,表明Ni(2)中心的配位是相当稳定的 ,这与实验规律相符。     

Energetic and topological analyses of the oxidation reaction between Mo(n) (n = 1, 2) and N2O

The interaction between molybdenum, atom, and dimer, with nitrous oxide has been investigated using density functional theory. The analysis of the potential energy surfaces for both reactions has revealed that a single molybdenum atom can activate the N--O bond of N2O requiring a small activation energy. However, the presence of several intersystem crossings between three different spin states, namely, septet, quintet and triplet states, seems to be the major constraint to the Mo + N2O reaction. Contrarily, the low-lying excited states (triplet and quintet) do not participate in the reaction between the molybdenum dimer and N2O. The latter reaction fully evolves on the singlet spin surface. Three different regions have been distinguished along the pathway: formation of an adduct complex, formation of an inserted compound, and the N2 detachment. The connection between the two first regions has been characterized by the formation of a special complex in which the N--O bond is so weakened that it could be considered as a first step in the insertion process. It has been shown that the topological changes along the pathways provide a clear explanation for the geometrical changes that occur along the reaction pathway. In summary, the detachment of the N2 molecule is found to be kinetically an effective process for both reactions, owing to the high exothermicity and consequently to the high internal energy of the insertion intermediates. However, in the case of Mo atom, the reaction should be a slow process due to the presence of spin-forbidden transitions. These results fully agree with previous experimental works.     

Theory of chemical bonds in metalloenzymes XII: Electronic and spin structures of metallo–oxo and isoelectronic species and spin crossover phenomena in oxygenation reactions

The broken-symmetry (BS) and multideterminant approaches to atomic oxygen (O), molecular oxygen (O₂) and iron–oxo (Fe(IV)O) core in P450 have elucidated electronic structures of the ground triplet and excited singlet states, which indicate isoelectronic characteristics of the species. The dissociation processes of the O–O and Fe–O double bonds are also examined to clarify the radical character, namely O-atom property responsible for radical mechanism of hydroxylations of alkanes and epoxidation of alkenes. This isolobal analogy has indeed enabled us to propose possible reaction mechanisms of oxygenation reactions by the Fe(IV)O species on the basis of available theoretical and experimental results for O and O₂. Similarly, an isolobal analogy of the σ^* bond among Fe(IV)O, dioxirane, peracids, etc. indicates the common electrophilic property for the oxygenation reactions. The small energy gaps between the high- and low-spin states of the transition structures and intermediates generated in the oxygenation reactions are found to be origins for spin crossover phenomena along the reaction pathways of these reactions.     

Theoretical study of inverted sandwich type complexes of 4d transition metal elements: interesting similarities to and differences from 3d transition metal complexes

Inverted sandwich type complexes (ISTCs) of 4d metals, (μ-η(6):η(6)-C(6)H(6))[M(DDP)](2) (DDPH = 2-{(2,6-diisopropylphenyl)amino}-4-{(2,6-diisopropylphenyl)imino}pent-2-ene; M = Y, Zr, Nb, Mo, and Tc), were investigated with density functional theory (DFT) and MRMP2 methods, where a model ligand AIP (AIPH = (Z)-1-amino-3-imino-prop-1-ene) was mainly employed. When going to Nb (group V) from Y (group III) in the periodic table, the spin multiplicity of the ground state increases in the order singlet, triplet, and quintet for M = Y, Zr, and Nb, respectively, like 3d ISTCs reported recently. This is interpreted with orbital diagram and number of d electrons. However, the spin multiplicity decreases to either singlet or triplet in ISTC of Mo (group VI) and to triplet in ISTC of Tc (group VII), where MRMP2 method is employed because the DFT method is not useful here. These spin multiplicities are much lower than the septet of ISTC of Cr and the nonet of that of Mn. When going from 3d to 4d, the position providing the maximum spin multiplicity shifts to group V from group VII. These differences arise from the size of the 4d orbital. Because of the larger size of the 4d orbital, the energy splitting between two d(δ) orbitals of M(AIP) and that between the d(δ) and d(π) orbitals are larger in the 4d complex than in the 3d complex. Thus, when occupation on the d(δ) orbital starts, the low spin state becomes ground state, which occurs at group VI. Hence, the ISTC of Nb (group V) exhibits the maximum spin multiplicity.     

A direct method for locating minimum-energy crossing points (MECPs) in spin-forbidden transitions and nonadiabatic reactions

An efficient computational method for locating minimum-energy crossing points (MECPs) between potential-energy surfaces in spin-crossover transitions and nonadiabatic spin-forbidden (bio)chemical reactions is introduced. The method has been tested on the phenyl cation and the computed MECP associated with its radiationless singlet-triplet spin crossover is in good agreement with available data. However, the convergence behavior of the present method is significantly more efficient than some alternative methods which allows us to study nonadiabatic processes in larger systems such as spin crossover in metal-containing compounds. The convergence rate of the method obeys a fast logarithmic law which has been verified on the phenyl cation. As an application of this new methodology, the MECPs of the ferrous complex [Fe(ptz)(6)](BF(4))(2), which exhibits light-induced excited spin state trapping, have been computed to identify their geometric and energetic parameters during spin crossover. Our calculations, in conjunction with spin-unrestricted density-functional calculations, show that the transition from the singlet ground state to a triplet intermediate and to the quintet metastable state of [Fe(ptz)(6)](BF(4))(2) is accompanied by unusually large bond-length elongations of the axial ligands ( approximately 0.26 and 0.23 A, respectively). Our results are consistent with crystallographic data available for the metastable quintet but also predict new structural and energetic information about the triplet intermediate and at the MECPs which is currently not available from experiment.     

O2 Binding to Heme is Strongly Facilitated by Near‐Degeneracy of Electronic States

This paper reports the computed O₂ binding to heme, which for the first time explains experimental enthalpies for this process of central importance to bioinorganic chemistry. All four spin states along the relaxed Fe? O₂‐binding curves were optimized using the full heme system with dispersion, thermodynamic, and scalar‐relativistic corrections, applying several density functionals. When including all these physical terms, the experimental enthalpy of O₂ binding (?59 kJ mol^?1) is closely reproduced by TPSSh‐D3 (?66 kJ mol^?1). Dispersion changes the potential energy surfaces and leads to the correct electronic singlet and heptet states for bound and dissociated O₂. The experimental activation enthalpy of dissociation (～82 kJ mol^?1) was also accurately computed (～75 kJ mol^?1) with an actual barrier height of ～60 kJ mol^?1 plus a vibrational component of ～10 and ～5 kJ mol^?1 due to the spin‐forbidden nature of the process, explaining the experimentally observed difference of ～20 kJ mol^?1 in enthalpies of binding and activation. Most importantly, the work shows how the nearly degenerate singlet and triplet states increase crossover probability up to ～0.5 and accelerate binding by ～100 times, explaining why the spin‐forbidden binding of O₂ to heme, so fundamental to higher life forms, is fast and reversible.