Theoretical study of the mechanism of valence tautomerism in cobalt complexes

Valence tautomerism is studied in the [Co(II-HS)(sq)(2)(bpy)]/[Co(III-LS)(sq)(cat)(bpy)] mononuclear cobalt complex by using DFT methods (HS, high spin; LS, low spin; cat, catecholate; sq, semiquinone; bpy, 2,2'-bipyridine). Calculations at the B3LYP* level of theory reproduce well the energy gap between the Co(II-HS) and Co(III-LS) forms giving an energy gap of 4.4 kcal/mol, which is comparable to the experimental value of 8.9 kcal/mol. Potential energy surfaces and crossing seams of the electronic states of the doublet, quartet, and sextet spin states are calculated along minimum energy paths connecting the energy minima corresponding to the different spin states. The calculated minimum energy crossing points (MECPs) are located at 8.8 kcal/mol in the doublet/sextet surfaces, at 10.2 kcal/mol in the doublet/quartet surfaces, and at 8.4 kcal/mol in the quartet/sextet surfaces relative to the doublet ground state. Considering the energy of the three spin states and the crossing points, the one-step relaxation mechanism between the Co(II-HS) and Co(III-LS) forms is the most probable. This research shows that mapping MECPs can be a useful strategy to analyze the potential energy surfaces of systems with complex deformation modes.     

Principal active species of horseradish peroxidase, compound I: a hybrid quantum mechanical/molecular mechanical study

The active species, Compound I, of horseradish peroxidase (HRP) has been investigated by quantum mechanical/molecular mechanical (QM/MM) calculations using 10 different QM regions. In accord with experimental data, the lowest doublet and quartet states are found to be virtually degenerate, with two unpaired electrons on the FeO moiety and one localized on the porphyrin in an a(2u)-dominant orbital with a minor, but nonnegligible, a(1u) component. The proximal ligand appears to be imidazole rather than imidazolate. The hydrogen-bonding network around the FeO moiety (i.e., Arg38 and His42) has significant influence on the axial bonds and the spin density distribution in the FeO moiety. Including this network in the QM region was found to be essential for reproducing the experimental M?ssbauer parameters. The protein environment shapes most of the subtle features of Compound I of HRP.     

What external perturbations influence the electronic properties of catalase compound I?

We have performed density functional theory calculations on an active-site model of catalase compound I and studied the responses of the catalytic center to external perturbations. Thus, in the gas phase, compound I has close-lying doublet and quartet spin states with three unpaired electrons: two residing in pi(FeO) orbitals and the third on the heme. The addition of a dielectric constant to the model changes the doublet-quartet energy ordering but keeps the same electronic configuration. By contrast, the addition of an external electric field along one of the principal axes of the system can change the doublet-quartet energy splitting by as much as 6 kcal mol(-1) in favor of either the quartet or the doublet spin state. This sensitivity is much stronger than the effect obtained for iron heme models with thiolate or imidazole axial ligands. Moreover, an external electric field is able to change the electronic system from a heme-based radical [Fe=O(Por*+)OTyr-] to a tyrosinate radical [Fe=O(Por)OTyr*]. This again shows that oxo-iron heme systems are chameleonic species that are influenced by external perturbations and change their character and catalytic properties depending on the local environment.     

气相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,导致反应势能面的势垒明显降低．     

Theoretical study on the spin-state energy splittings and local spin in cationic [Re]-Cn-[Re] complexes

Total spin-state energy splittings are calculated for mono- and dications of the formula {[Re]-Cn-[Re]}z+ where [Re] = eta5-(C5Me5)Re(NO)(PPh3). Cn is an even-numbered carbon chain with n ranging from 4 to 20, and z is 1 or 2. These complexes are experimentally known, and their potential role as molecular electronic devices initiated this work. We have considered the different total spin states monocation/doublet, monocation/quartet, dication/singlet, and dication/triplet. Data obtained for two density functionals BP86 and B3LYP were compared to verify the internal consistency of the results. In both ionization states, the low-spin state is the ground state, but the spin-state splittings decrease as the chain gets longer. For the dications, the splitting reaches a nearly constant value of about 10 kJ/mol with BP86 and about 4 kJ/mol with B3LYP when there are at least 14 carbon atoms in the chain, whereas for the monocations, no constant value appears to be reached asymptotically, not even if 20 carbon atoms are in the chain. For monocations, the splittings range from 138 kJ/mol (n = 4) to 68 kJ/mol (n = 20) with BP86 and from 134 kJ/mol (n = 4) to 73 kJ/mol (n = 20) with B3LYP and are thus considerably higher than those of the dications. The spin-state splittings are qualitatively mirrored by the energy splitting between the highest-occupied molecular orbital with beta spin (HOMObeta) and the lowest-unoccupied molecular orbital with alpha spin (LUMOalpha) as obtained in the low-spin state. Furthermore, the HOMOalpha-LUMOalpha gaps decrease as the carbon chain lengthens. In addition, the local distribution of the ?z expectation value is analyzed for the monocation/doublet, the monocation/quartet, and the dication/triplet state using a modified L?wdin partitioning scheme. In the monocation/doublet and the dication/triplet state, the electron spin is distributed mainly on the metal centers and slightly delocalized onto the carbon chain. In the monocation/quartet state for chain lengths of more than 8 carbon atoms, the electron spin is mainly localized on selected atoms of the chain and not on the metal centers. In all cases, the spin delocalization onto the chain increases as the chain gets longer.     

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．     

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.     

First findings of monocrystalline aragonite inclusions in garnet from diamond-grade UHPM rocks (Kokchetav Massif, Northern Kazakhstan)

The crystal field effect and microscopic origins of the Zeeman g-factors g(//) and g(⊥) for (6)S(3d(5)) state ions at tetragonal symmetry crystal filed, taking into account the spin-spin (SS), the spin-other-orbit (SOO), and the orbit-orbit (OO) magnetic interactions besides the well-known spin-orbit (SO) magnetic interaction, have been investigated using the microscopic spin Hamiltonian theory and the complete diagonalization method (CDM). It is found that the g(//)(±1/2)≠g(//)(±5/2) and g(⊥)(±1/2)≠g(⊥)(±5/2), where the g-factors g(//)(±1/2) and g(⊥)(±1/2) express the g-factors of the ground state |M?(s)=±1/2), whereas the g-factors g(//)(±5/2) and g(⊥)(±5/2) express the g-factors of the ground state |M?(s)=±5/2). It is shown that although the SO magnetic interaction is the most important one, the contributions to the shifts of g-factors Δg(//)(=2.0023-g(//)) and Δg(⊥)(=2.0023-g(⊥)) from other three magnetic interactions including the SS, SOO, and OO magnetic interactions are appreciable and should not be omitted, especially for the shifts of g-factors Δg(//)(±5/2) and Δg(⊥)(±5/2). The individual contributions to the shifts of g-factors arising from the spin quartet states and spin doublet states have been studied. The investigations show that the Δg(//)(±1/2) and Δg(⊥)(±1/2) primarily result from the spin quartet states, whereas Δg(//)(±5/2) and Δg(⊥)(±5/2) from the spin quartet states as well as the combined effects between the spin quartet states and the spin doublet states. The contribution to the shifts of g-factors from the net spin doublet states is zero.     

Electron paramagnetic resonance of the excited states of the three-spins-1/2 ZnTPP–3-NOPy complex

EPR spectra of the excited quartet and doublet molecular states of (tetraphenylporphinato)zinc(II)covalently bounded to 3-(N-nitronyl-notroxide) pyridine stable radical are modeled in terms of the spin-Hamiltonian given by the sum of the contributions from the radical and triplet moieties, and the interaction between them. The later is represented by anisotropic point dipolar and isotropic exchange electron spin–spin interactions. It is shown that the high field (W-band) EPR spectra depend on energy separation between the electronic doublet (D) and quartet (Q) states. This dependence was utilized to estimate the upper limit of the intensity of exchange interaction between the radical and porphyrin moieties.     

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.     

U+与CO2气相反应的密度泛函理论研究

运用密度泛函理论(DFT)中的B3LYP方法,U原子用含相对论有效原子实势(ECP)校正的基组(SDD),C、O原子采用6-311+G(d)基组,对气相中U+和CO2的反应进行了理论研究.通过研究二重和四重自旋态的反应势能面(PESs),优化得到了两条反应路径的反应物、中间体、过渡态和产物的结构.用"两态反应"(TSR)分析反应机理,结果表明体系的优先选择路径为高自旋态进入和低自旋态离开反应,发生在四重态和二重态的自旋多重度的改变使得整个反应系统能以一个低能反应途径进行.     

Theoretical studies of the electronic spectrum of SiC

Ab initio based multireference singles and doubles configuration interaction calculations have been carried out to study the electronic structure and spectroscopic properties of the SiC⁺ ion. Potential energy curves and spectroscopic constants (r_e, T_e, ω_e) of 14 low-lying doublet and quartet states of the ion are studied. The spin-orbit coupling has been included to see its effects on the spectroscopic properties. Transition probabilities of some quartet–quartet transitions are computed, while the spin-forbidden transitions are very weak. Dipole moments of all low-lying states are estimated by keeping the origin at the center of mass. The vertical and adiabatic ionization energies of SiC are also reported.     

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.     

On the third virial coefficient for the alkali metal vapours

Theoretical calculations of the third virial coefficient for the pure components of all alkali metal vapours have been carried out using recently reported potential energy surfaces for the doublet and quartet states in which the interaction of three ²S ground-state atoms may evolve. The discrepancy between the theoretical and available experimental estimates for those coefficients is pointed out.     

Ab initio study of low-lying electronic states of SnCl2+

Complete active space self-consistent field (CASSCF), multireference configuration interaction (MRCI), and restricted-spin coupled-cluster singles-doubles with perturbative triples [RCCSD(T)] calculations have been carried out on low-lying doublet and quartet states of SnCl2+, employing basis sets of up to aug-cc-pV5Z quality. Effects of core correlation and off-diagonal spin-orbit interaction on computed vertical ionization energies were investigated. The best theoretical estimate of the adiabatic ionization energy (including zero-point vibrational energy correction) to the X2A1 state of SnCl2+ is 10.093+/-0.010 eV. The first photoelectron band of SnCl2 has also been simulated by employing RCCSD(T)/aug-cc-pV5Z potential energy functions and including Duschinsky rotation and anharmonicity.     

Combined quantum mechanical/molecular mechanical study on the pentacoordinated ferric and ferrous cytochrome P450cam complexes

The pentacoordinated ferric and ferrous cytochrome P450(cam) complexes have been investigated by combined quantum mechanical/molecular mechanical (QM/MM) calculations in the presence of a protein/solvent environment and by QM calculations on the isolated QM regions with use of density functional theory. The B3LYP functional has been found more reliable than the BLYP and BHLYP functionals for estimating the relative state energies. The B3LYP/CHARMM calculations with an all-electron basis set for iron give high-spin ground states for the title complexes, in agreement with experiment. The comparison of the B3LYP/CHARMM results of the entire protein system with the B3LYP calculations on the naked QM regions shows that the amount of stabilization by the protein environment is largest for the intermediate-spin states, followed by the high-spin states of the complexes. The calculation of M?ssbauer parameters in the presence of the enzyme environment confirms the double occupation of the d(xz) orbital in the quintet spin state of the ferrous complex, consistent with the computed QM/MM energies in the enzyme environment, while the d(x)2(-)(y)2 orbital is doubly occupied in the gas-phase quintet state.     

Theoretical study of the electronic states of Nb4, Nb5 clusters and their anions (Nb4-,Nb5-)

Geometries and energy separations of the various low-lying electronic states of Nb(n) and Nb(n) (-) (n=4,5) clusters with various structural arrangements have been investigated. The complete active space multiconfiguration self-consistent field method followed by multireference singles and doubles configuration interaction (MRSDCI) calculations that included up to 52x10(6) configuration spin functions have been used to compute several electronic states of these clusters. The ground states of both Nb(4) ((1)A('), pyramidal) and Nb(4) (-) ((2)B(3g), rhombus) are low-spin states at the MRSDCI level. The ground state of Nb(5) cluster is a doublet with a distorted trigonal bipyramid (DTB) structure. The anionic cluster of Nb(5) has two competitive ground states with singlet and triplet multiplicities (DTB). The low-lying electronic states of these clusters have been found to be distorted due to Jahn-Teller effect. On the basis of the energy separations of our computed electronic states of Nb(4) and Nb(5), we have assigned the observed photoelectron spectrum of Nb(n) (-) (n=4,5) clusters. We have also compared our MRSDCI results with density functional calculations. The electron affinity, ionization potential, dissociation and atomization energies of Nb(4) and Nb(5) have been calculated and the results have been found to be in excellent agreement with the experiment.     

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.     

Electronic properties of vanadium-doped TiO2

The electronic properties of vanadium-doped rutile TiO(2) are investigated theoretically with a Hartree-Fock/DFT hybrid approach. The most common oxidation states (V(2+), V(3+), V(4+), and V(5+)) in different spin states are investigated and their relative stability is calculated. The most stable spin states are quartet, quintet, doublet, and singlet for V(2+), V(3+), V(4+), and V(5+) doping, respectively. By comparing the formation energy with respect to the parent oxides and gas-phase oxygen (ΔE), we conclude that V(4+) (ΔE=145.3 kJ mol(-1)) is the most likely oxidation state for vanadium doping with the possibility of V(5+) doping (ΔE=283.5 kJ mol(-1)). The energetic and electronic properties are converged with dopant concentrations in the range of 0.9 to 3.2%, which is within the experimentally accessible range. The investigation of electronic properties shows that V(4+) doping creates both occupied and unoccupied vanadium states in the band gap and V(5+) doping creates unoccupied states at the bottom of the conduction band. In both cases there is a significant reduction of the band gap by 0.65 to 0.75 eV compared to that of undoped rutile TiO(2).     

Geometries and stabilities of the carbon clusters with the rhodium impurity: a computational investigation

A density functional study of the RhCn(n = 1-6) clusters with different spin states has been carried out systematically by using the B3LYP/Lan2DZ method. The equilibrium geometries associated with total energies and natural populations of RhCn (n = 1-6) clusters are calculated and presented. Stabilities and electronic properties are discussed in detail. The relative stabilities in term of the calculated fragmentation energies show that the lowest-energy RhCn clusters with rhodium atom being located at terminal of carbon chain are the linear geometries and the ground states of the RhCn clusters alternate between doublet (for n-odd members) and quartet (for n-even members) states. Furthermore, the calculated fragmentation energies of the RhCn show strong even-odd alternations: the RhCn clusters with an odd number of carbon atoms are more stable than those with an even number ones. In addition, we comment on the charge transfer and chemical bonding properties within the clusters.