Theoretical investigations of spin-orbit coupling and kinetics in reaction W + NH3 → N≡WH3

The hydrogen-transfer reaction of W + NH(3) incorporates four possible diabatic reaction pathways along with septet, quintet, triplet, and singlet states. The intersystem crossings thus play an important role in the reaction mechanisms. In this work, ab initio and DFT methods are used to determine all possible intermediates, transition states, products, and intersystem crossing points as well as the spin-orbit couplings. The mechanism of hydrogen elimination is further revealed by the natural bond orbital analysis. From the rate constants yielded by a nonadiabatic transition state theory, we find that two intersystem crossings significantly change the reaction pathways. Finally, we suggest a feasible reaction pathway with exothermicity 72.8 kcal/mol, which is consistent with the experimental measurements.     

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.     

New insight into the gas-phase bimolecular self-reaction of the HOO radical

The singlet and triplet potential energy surfaces (PESs) for the gas-phase bimolecular self-reaction of HOO*, a key reaction in atmospheric environments, have been investigated by means of quantum-mechanical electronic structure methods (CASSCF and CASPT2). All the reaction pathways on both PESs consist of a first step involving the barrierless formation of a prereactive doubly hydrogen-bonded complex, which is a diradical species lying about 8 kcal/mol below the energy of the reactants at 0 K. The lowest energy reaction pathway on both PESs is the degenerate double hydrogen exchange between the HOO* moieties of the prereactive complex via a double proton transfer mechanism involving an energy barrier of only 1.1 kcal/mol for the singlet and 3.3 kcal/mol for the triplet at 0 K. The single H-atom transfer between the two HOO* moieties of the prereactive complex (yielding HOOH + O2) through a pathway keeping a planar arrangement of the six atoms involves a conical intersection between either two singlet or two triplet states of A' and A" symmetries. Thus, the lowest energy reaction pathway occurs via a nonplanar cisoid transition structure with an energy barrier of 5.8 kcal/mol for the triplet and 17.5 kcal/mol for the singlet at 0 K. The simple addition between the terminal oxygen atoms of the two HOO* moieties of the prereactive complex, leading to the straight chain H2O4 intermediate on the singlet PES, involves an energy barrier of 7.3 kcal/mol at 0 K. Because the decomposition of such an intermediate into HOOH + O2 entails an energy barrier of 45.2 kcal/mol at 0 K, it is concluded that the single H-atom transfer on the triplet PES is the dominant pathway leading to HOOH + O2. Finally, the strong negative temperature dependence of the rate constant observed for this reaction is attributed to the reversible formation of the prereactive complex in the entrance channel rather than to a short-lived tetraoxide intermediate.     

Spin-inversion mechanisms in O2 binding to a model heme complex revisited by density function theory calculations

Spin-inversion mechanisms in O₂ binding to a model heme complex, consisting of Fe(II)-porphyrin and imidazole, were investigated using density-functional theory calculations. First, we applied the recently proposed mixed-spin Hamiltonian method to locate spin-inversion structures between different total spin multiplicities. Nine spin-inversion structures were successfully optimized for the singlet–triplet, singlet–quintet, triplet–quintet, and quintet–septet spin-inversion processes. We found that the singlet–triplet spin-inversion points are located around the potential energy surface region at short Fe–O distances, whereas the singlet–quintet and quintet–septet spin-inversion points are located at longer Fe–O distances. This suggests that both narrow and broad crossing models play roles in O₂ binding to the Fe-porphyrin complex. To further understand spin-inversion mechanisms, we performed on-the-fly Born-Oppenheimer molecular dynamics calculations. The reaction coordinates, which are correlated to the spin-inversion dynamics between different spin multiplicities, are also discussed.     

Computational studies of the thermal fragmentation of P-arylphosphiranes: have arylphosphinidenes been generated by this method?

CASSCF, CASPT2, CCSD(T), and (U)B3LYP electronic structure calculations have been performed in order to investigate the thermal fragmentation of P-phenylphosphirane (1) to phenylphosphinidene (PhP) and ethylene. The calculations show that generation of triplet PhP via a stepwise pathway is 21 kcal mol(-1) less endothermic and has a 12 kcal mol(-1) lower barrier height than concerted fragmentation of 1 to give singlet PhP. The formation of singlet PhP via a concerted pathway is predicted to be stereospecific, whereas formation of triplet PhP is predicted to occur with complete loss of stereochemistry. However, calculations on fragmentation of anti-cis-2,3-dimethyl-P-mesitylphosphirane (cis-1Me) to triplet mesitylphosphinidene (MesP) indicate that this reaction should be more stereospecific, in agreement with the experimental results of Li and Gaspar. Nevertheless, with a predicted free energy of activation of 42 kcal mol(-1), the formation of MesP from cis-1Me is not likely to have occurred in an uncatalyzed reaction at the temperatures at which this phosphirane has been pyrolyzed.     

Spin-forbidden CO binding to iron–sulfur cluster-free hydrogenase: A density functional study

Spin-forbidden CO binding to the iron–sulfur cluster-free hydrogenase (Hmd) is studied by the DFT calculation. The result shows that the surface of the triplet causes a PHmd–CO minimum and that ^3,5MECP is the lowest energy path to PHmd–CO. It is found that this CO binding involves a low barrier of 0.931 kcal/mol because of the need to change from a bound triplet state to the Hmd quintet ground state.     

(Salen)Mn(III)-catalyzed epoxidation reaction as a multichannel process with different spin states. Electronic tuning of asymmetric catalysis: a theoretical study

The (salen)Mn(III)-catalyzed epoxidation reaction mechanism has been investigated using density functional theory (DFT). There is considerable interest in and controversy over the mechanism of this reaction. The results of experimental studies have offered some support for three different reaction mechanisms: concerted, stepwise radical, and metallooxetane mediated. In this paper, a theoretical examination of the reaction suggests a novel mechanism that describes the reaction as a multichannel process combining both concerted and stepwise radical pathways. The competing channels have different spin states: the singlet, the triplet, and the quintet. The singlet reaction pathway corresponds to a concerted mechanism and leads exclusively to a cis epoxide product. In contrast, the triplet and quintet reactions follow a stepwise mechanism and lead to a product mixture of cis and trans epoxides. We show that the experimentally observed dependence of isomer product ratios on electronic effects connected with the substitution of the catalyst ligands is due to changing the relative position and, hence, the relative activities of the channels with different cis-trans yields. Because the results and conclusions of the present work dramatically differ from the results and conclusion of the recent DFT theoretical investigation (Linde, C.; Akermark, B; Norrby, P.-O.; Svensson, M. J. Am. Chem. Soc. 1999, 121, 5083.), we studied possible sources for the deep contradictions between the two works. The choice of the DFT functional and a model has been shown to be crucial for accurate results. Using high level ab initio calculations (coupled cluster-CCSD(T)), we show that the computational procedure employed in this study generates significantly more reliable numerical results. It is also shown that the smaller cationic model without a chlorine ligand that was used by Linde et al. is too oversimplified with respect to our larger neutral model. For this reason, using the cationic model led to a qualitatively wrong quintet reaction profile that played a key role in theoretical postulates in the earlier work.     

Insights into the structure and reactivity of acylperoxo complexes in the Kochi-Jacobsen-Katsuki catalytic system. A density functional study

Structural properties of the acylperoxo complexes [(Salen)Mn(III)RCO(3)] (2) and [(Salen)Mn(IV)RCO(3)] (3), the critical intermediates in the Kochi-Jacobsen-Katsuki reaction utilizing organic peracids or O(2)/aldehydes as oxygen source, have been studied with the density functional theory. Four distinct isomers, cis(O,N), cis(N,O), cis(N,N), and trans, of these complexes have been located. The isomer 2-cis(O,N) in its quintet ground state, and nearly degenerate isomers 3-cis(O,N) and 3-cis(N,O) in their quartet ground states are found to be the lowest in energy among the other isomers. The O-O bond cleavage in the cis(O,N), cis(N,O), and trans isomers of 2 and 3 has been elucidated. In complex 3, the O-O bond is inert. On the contrary, in complex 2, the O-O bond cleaves via two distinct pathways. The first pathway occurs exclusively on the quintet potential energy surface (PES) and corresponds to heterolytic O-O bond scission coupled with insertion of an oxygen atom into an Mn-N(Salen) bond to form 2-N-oxo species; this pathway has the lowest barrier of 14.9 kcal/mol and is 15.6 kcal/mol exothermic. The second pathway is tentatively a spin crossover pathway. In particular, for 2-cis(O,N) and 2-cis(N,O) the second pathway proceeds through a crucial minimum on the seam of crossing (MSX) between the quintet and triplet PESs followed by heterolytic O-O cleavage on the triplet PES, and produces unusual triplet 2-cis(O,N)- and 2-cis(N,O)-oxo ([(Salen)Mn(V)(O)RCO(2)]) species; this pathway requires 12.8 kcal/mol and is 1.4 kcal/mol endothermic. In contrast, for the 2-trans isomer, spin crossing is less crucial and the O-O cleavage proceeds homolytically to generate 2-trans-oxo [(Salen)Mn(IV)(O)] species with RCO(2) radical; this pathway, however, cannot compete with that in 2-cis because it needs 21.9 kcal/mol for activation and is 15.3 kcal/mol endothermic. In summary, the O-O cleavage occurs predominantly in the 2-cis complexes, and may proceed either through pure high spin or spin crossover heterolytic pathway to produce 2-cis-oxo and 2-N-oxo species.     

The electronic structure and energetics of V(+)-benzene half-sandwiches of different multiplicities: Comparative multireference and single-reference theoretical study

Multireference [complete active space self-consistent field (CASSCF) and multiconfigurational quasidegenerate perturbation theory (MCQDPT)] and single-reference ab initio (Moller-Plesset second order perturbation theory (MP2) and coupled clusters with singles, doubles and noniterative triples [CCSD(T)]) and density functional theory (PBE and B3LYP) electronic structure calculations of V(C(6)H(6))(+) half-sandwich in the states of different multiplicities are described and compared. Detailed analyses of the geometries and electronic structures of the all found states are given; adiabatic and diabatic dissociation energies are estimated. The lowest electronic state of V(C(6)H(6))(+) half-sandwich was found to be the quintet (5)B(2) state with a slightly deformed upside-down-boat-shaped benzene ring and d(4) configuration of V atom, followed by a triplet (3)A(2) state lying about 4 kcal/mol above. The lowest singlet state (1)A(1)(d(4)) lies much ( approximately 28 kcal/mol) higher. MCQDPT calculated adiabatic dissociation energy (53.6 kcal/mol) for the lowest (5)B(2)(d(4)) state agrees well with the current 56.4 (54.4) kcal/mol experimental estimate, giving a preference to the lower one. Compared to MCQDPT, B3LYP hybrid exchange-correlation functional provides the best results, while CCSD(T) performs usually worse. Gradient-corrected PBE calculations tend to systematically overestimate metal-benzene binding in the row quintet相似文献   

Quantum chemical study of the mechanism of reaction between NH (X 3sigma-) and H2, H2O, and CO2 under combustion conditions

Reactions of ground-state NH (3sigma-) radicals with H2, H2O, and CO2 have been investigated quantum chemically, whereby the stationary points of the appropriate reaction potential energy surfaces, that is, reactants, products, intermediates, and transition states, have been identified at the G3//B3LYP level of theory. Reaction between NH and H2 takes place via a simple abstraction transition state, and the rate coefficient for this reaction as derived from the quantum chemical calculations, k(NH + H2) = (1.1 x 10(14)) exp(-20.9 kcal mol(-1)/RT) cm3 mol(-1) s(-1) between 1000 and 2000 K, is found to be in good agreement with experiment. For reaction between triplet NH and H2O, no stable intermediates were located on the triplet reaction surface although several stable species were found on the singlet surface. No intersystem crossing seam between triplet NH + H2O and singlet HNO + H2 (the products of lowest energy) was found; hence there is no evidence to support the existence of a low-energy pathway to these products. A rate coefficient of k(NH + H2O) = (6.1 x 10(13)) exp(-32.8 kcal mol(-1)/RT) cm3 mol(-1) s(-1) between 1000 and 2000 K for the reaction NH (3sigma-) + H2O --> NH2 (2B) + OH (2pi) was derived from the quantum chemical results. The reverse rate coefficient, calculated via the equilibrium constant, is in agreement with values used in modeling the thermal de-NO(x) process. For the reaction between triplet NH and CO2, several stable intermediates on both triplet and singlet reaction surfaces were located. Although a pathway from triplet NH + CO2 to singlet HNO + CO involving intersystem crossing in an HN-CO2 adduct was discovered, no pathway of sufficiently low activation energy was discovered to compare with that found in an earlier experiment [Rohrig, M.; Wagner, H. G. Proc. Combust. Inst. 1994, 25, 993.].     

Triplet ground state derivative of aza-m-xylylene diradical with large singlet-triplet energy gap

Organic molecules with a strong preference for triplet ground states, in which the triplet state is below the lowest singlet state by ≥10 kcal/mol, are typically short-lived and mostly detected as reactive intermediates. We now report a triplet ground state derivative of aza-m-xylylene diradical with a large singlet-triplet energy gap (ΔE(ST)) of ～10 kcal/mol, which is comparable to ΔE(ST) for the well-known reactive intermediate m-xylylene diradical. The aminyl diradical persists in solution at room temperature on the time scale of minutes.     

Reaction mechanism of the CCN radical with nitric oxide

To investigate the possibility of the carbyne radical CCN in removal of nitric oxide, a detailed computational study is performed at the Gaussian-3//B3LYP/6-31G(d) level on the CCN + NO reaction by constructing the singlet and triplet electronic state [C(2)N(2)O] potential energy surfaces (PESs). The barrierless formation of the chain-like isomers NCCNO (singlet at -106.5, triplet cis at -48.2 and triplet trans at -47.6 kcal/mol) is the most favorable entrance attack on both singlet and triplet PESs. Subsequently, the singlet NCCNO takes an O-transfer to form the branched intermediate singlet NCC(O)N (-85.6), which can lead to the fragments CN + NCO (-51.2) via the intermediate singlet NCOCN (-120.3). The simpler evolution of the triplet NCCNO is the direct N-O rupture to form the weakly bound complex triplet NCCN...O (-56.2) before the final fragmentation to NCCN + (3)O (-53.5). However, the lower lying products (3)NCN + CO (-105.6) and (3)CNN + CO (-74.6) are kinetically much less competitive. All the involved transition states for generation of CN + NCO and NCCN + (3)O lie much lower than the reactants. Thus, the novel reaction CCN + NO can proceed effectively even at low temperatures and is expected to play a role in both combustion and interstellar processes. Significant differences are found on the singlet PES between the CCN + NO and CH + NO reaction mechanisms.     

Exploring the lowest energy triplet potential energy surface for cyclic C4H4 isomers with the complete basis set ab initio method. Is the transformation of triafulvene into cyclobutadiene possible in their excited states?

High level computational studies were performed with the aim being to explore the possibility of converting the experimentally available triafulvene into the hard-to-detect cyclobutadiene. The method is based on the simple approach used to excite triafulvene into the triplet state, and then through various reaction channels, come to the aromatic triplet cyclobutadiene. Triplet cyclobutadiene is only a few kcal/mol higher in energy than singlet cyclobutadiene and should be easily relaxed into the latter. Several reaction pathways that include only a concerted mechanism, as well as reaction pathways that include the radical formation–recombination were also explored. Some possible approaches for experimentally obtaining the singlet or triplet cyclobutadiene were suggested.     

Mechanistic photodissociation of CO-ligated neuroglobin and subsequent rebinding processes: a theoretical study

In the present work, density functional theory (QM) and molecular mechanics (MM) method were used to study mechanistic photodissociation of CO-ligated neuroglobin (Ngb-CO). It was found that all the electronic states investigated here are bound with respect to the Fe-CO separation, except for a couple of near-degenerate states (1E) that are repulsive. Irradiation of Ngb-CO at 533 nm leads to the system in the lowest two excited singlet states (1Q), where non-adiabatic CO dissociation proceeds with high efficiency through the intersection between 1Q and 1E. Soret band (1B) is the strongest in the absorption spectra of Ngb-CO with the peak at 415 nm. The systems in the 1B states decay to the 1E states via fast internal conversion, which is followed by the CO dissociation. The CO dissociation induces a considerable change in the structure of the Ngb protein. The initial dissociation involves a rotation of CO, which is accompanied with movement of several residues. When the Fe-C distance is larger than a critical value of 3.0 A, the CO molecules transfer more freely into the cavity of the protein. The pentacoordinated heme was found to be a transient intermediate after CO dissociation.     

p-Phenylenecarbenonitrene and its halogen derivatives: how does resonance interaction between a nitrene and a carbene center affect the overall electronic configuration?

A series of para-conjugatively coupled phenylenecarbenonitrenes [(4-nitrenophenyl)methylene (3a), (4-nitrenophenyl)fluoromethylene (3b), (4-nitrenophenyl)chloromethylene (3c), and (4-nitrenophenyl)bromomethylene (3d)] were generated in argon matrix at low temperature (10 or 13 K) and characterized by IR and UV/vis spectroscopy. Density functional theory (B3LYP/6-31G(d)) and ab initio (MCSCF, CASPT2) methods were used to study the ground- and some low-lying excited states of 3a-d. The experimental and computational data suggest that 3a-d have singlet ground states (S0) and can be thought of as quinonoidal biradicals. In all cases, the lowest triplet (T1) and quintet (Q1) states lie about 2 kcal mol(-1) and 28-29 kcal mol(-1), respectively, higher in energy than S0. On the other hand the substituent is found to have a significant effect on the relative energy of the second excited triplet (T2) state. This state tends to become relatively more stable as the ability of the substituent to enforce a closed-shell configuration at the carbene subunit increases. Interestingly, the energy difference between the T2 and S0 states in 3a-d is found to depend linearly on the S-T gap of the corresponding phenylcarbenes 7a-d. This relationship is helpful in predicting when a substituted p-phenylenecarbenonitrene may have a triplet ground state instead of a singlet one.     

The almost bottleable triplet carbene: 2,6-dibromo-4-tert-butyl-2',6'-bis(trifluoromethyl)-4'-isopropyldiphenylcarbene.

Computations on 2,6-dibromo-4-tert-butyl-2',6'-bis(trifluoromethyl)-4'-isopropyldiphenylcarbene (1) using ab initio and density functional theory methods underscore the unusual stability of the triplet over the singlet state. At the B3LYP/6-311G(d,p) level, the triplet state had a slightly bent central C-C-C bond angle of 167 degrees, whereas this angle in the singlet was 134 degrees. The B3LYP singlet-triplet splitting (12.2 kcal/mol) was larger than that of the parent molecule (5.8 kcal/mol), diphenylcarbene (2), which also has a triplet ground state. The energy of a suitable isodesmic reaction showed the triplet and singlet states of (1) to be destabilized, by 6.3 and 12.5 kcal/mol, respectively, due to the combined effects of the CF3, Br, and alkyl substituents. The linear-coplanar form of (3)(1), which might facilitate dimerization or electrophilic attack at the more exposed diradical center, was prohibitively (35.9 kcal/mol) higher in energy. Our results confirm Tomioka's conclusion that the triplet diarylcarbene, ortho-substituted with bulky CF3 and Br substituents, is persistent due to steric protection of the diradical center. Dimerization and other possible reaction pathways are inhibited, not only by the bulky ortho substituents but also by the para alkyl groups. The increase in stability of the triplet ((3)(1)) state relative to the singlet ((1)(1)) state does not influence the reactivity directly.     

Relativistic DFT study on the reaction mechanism of second-row transition metal Ru with CO2

To estimate the importance of relativistic effects on the reaction mechanisms between Ru and CO2, the potential energy surfaces have been performed in the triplet and quintet electronic states using quasi-relativistic (Pauli), zero-order regularly approximated (ZORA), and nonrelativistic (NR) density functional theory (DFT) at the PW91/TZP level. The results demonstrate that there are two rival reaction mechanisms: one is an addition mechanism and the other is an insertion mechanism in the triplet state. The only mechanism in the quintet state is the insertion mechanism. The most favored reaction mechanism in Ru + CO2 is that the Ru atom in its ground state first attacks the CO bond of CO2, forming q-Ru(CO)O (5A') with the insertion mechanism, and then undergoes an intersystem crossing to t-Ru(CO)O (3A'). Then it crosses t-TS3 to produce t-ORuCO molecule. The relativistic effects are important for reactivity of the second-row transition metal to CO2. In the key step of t-Ru(CO)O via t-TS3 to t-ORuCO, relativistic effects reduce the barrier energy by 10.3 kcal/mol, which is nearly half the nonrelativistic barrier energy.     

Characterization of electronic structure and properties of a Bis(histidine) heme model complex

Ferric and ferrous hemes, such as those present in electron transfer proteins, often have low-lying spin states that are very close in energy. To explore the relationship between spin state, geometry, and cytochrome electron transfer, we investigate, using density functional theory, the relative energies, electronic structure, and optimized geometries for a high- and low-spin ferric and ferrous heme model complex. Our model consists of an iron-porphyrin axially ligated by two imidazoles, which model the interaction of a heme with histidine residues. Using the B3LYP hybrid functional, we found that, in the ferric model heme complex, the doublet is lower in energy than the sextet by 8.4 kcal/mol and the singlet ferrous heme is 6.7 kcal/mol more stable than the quintet. The difference between the high-spin ferric and ferrous model heme energies yields an adiabatic electron affinity (AEA) of 5.24 eV, and the low-spin AEA is 5.17 eV. Both values are large enough to ensure electron trapping, and electronic structure analysis indicates that the iron d(pi) orbital is involved in the electron transfer between hemes. M?ssbauer parameters calculated to verify the B3LYP electronic structure correlate very well with experimental values. Isotropic hyperfine coupling constants for the ligand nitrogen atoms were also evaluated. The optimized geometries of the ferric and ferrous hemes are consistent with structures from X-ray crystallography and reveal that the iron-imidazole distances are significantly longer in the high-spin hemes, which suggests that the protein environment, modeled here by the imidazoles, plays an important role in regulating the spin state. Iron-imidazole dissociation energies, force constants, and harmonic frequencies were calculated for the ferric and ferrous low-spin and high-spin hemes. In both the ferric and the ferrous cases, a single imidazole ligand is more easily dissociated from the high-spin hemes.     

Spin-forbidden ligand binding to the ferrous-heme group: ab initio and DFT studies

The potential energy surfaces (PESs) and associated energy barriers that characterize the spin-forbidden recombination reactions of the gas-phase ferrous deoxy-heme group with CO, NO, and H2O ligands have been calculated using density functional theory (DFT). The bond energy for binding of O2 has also been calculated. Extensive large basis set CCSD(T) calculations on two small models of the heme group have been used to calibrate the accuracy of different DFT functionals for treating these systems. Pure functionals are shown to overestimate the stability of the low-spin forms of the deoxy-heme model, and to overestimate the binding energy of H2O and CO, whereas hybrid functionals such as B3PW91 and B3LYP yield accurate results. Accordingly, the latter functionals have been used to explore the PESs for binding. CO binding is found to involve a significant barrier of ca. 3 kcal mol-1 due to the need to change from the deoxy-heme quintet ground state to the bound singlet state. Binding of water does not involve a barrier, but the resulting bond is weak and may be further weakened in the protein environment, which should explain why water binding is not usually observed in heme proteins such as myoglobin. NO binding involves a low barrier, which is consistent with observed rapid geminate recombination. The calculated bond energies are in good agreement with previous reported values and in fair agreement with experiment for CO and O2. The value for NO is significantly lower than the experimentally derived bond energy, suggesting that B3LYP is less accurate in this case.