Ab initio calculation of the lowest singlet and triplet states in CH2, CHF,CF2, and CHCH3

The lowest singlet and triplet states of the radicals CH₂, CHF, CF₂, and CHCH₃ have been investigated both in SCF and IEPA approximation (“independent electron pair approach” to account for electron correlation). The SCF calculations yield triplet ground states for CH₂, CHF, and CHCH₃, and a singlet ground state for CF₂. Electron correlation stabilizes the singlet state by about 14 kcal/mole with respect to the triplet for all four radicals leading to a singlet ground state also for CHF. The final triplet-singlet energy separations are 10, 6, ?11, ?47 kcal/mole for CH₂, CHCH₃, CHF, CF₂, respectively. Values for equilibrium bond angles, ionization potentials and bond energies are also given.     

Direct ab initio MD study on the hydrogen abstraction reaction of triplet state acetone from methanol molecule

Solvent re-orientation process of triplet acetone/methanol complex and intermolecular hydrogen atom abstraction reaction on the triplet state energy surface, (CH₃)₂C=O (T₁) + CH₃OH → (CH₃)₂C–OH + CH₂OH in gas phase, have been investigated by means of density functional theory (DFT) and direct ab initio molecular dynamics (MD) methods. The static DFT calculation of hydrogen abstraction reaction at the T₁ state showed that the transition state is 16.4 and 30.9 kcal/mol lower than the energy levels of S₁ and S₂ states, respectively, and 9.2 kcal/mol higher than the bottom of T₁ state. The product state, (CH₃)₂C–OH⋯CH₂OH, is 8.4 kcal/mol lower in energy than the level of T₁ state. The direct ab initio MD calculation showed that the product is rapidly formed within 150 fs and the separated products (CH₃)₂C–OH + CH₂OH were formed. The mechanism of reaction dynamics of the triplet acetone/methanol complex was discussed on the basis of theoretical results.     

RN (R＝CH3, CH3CH2)异构化及脱氢反应的量子拓扑研究

利用量子化学从头算CASSCF方法在6-311＋G (d, p)基组水平上对单线态和三线态RN (R＝CH₃, CH₃CH₂)异构化反应及RN脱氢反应的微观机理进行了理论研究. 在MP2/6-311＋G (d, p)和CCSD/6-311＋G (d, p)水平上进行了单点能校正. 单态和三态势能面的交叉点(ISC)的存在清楚地说明了基态反应物³RN异构化为基态产物¹R'NH (R'＝CH₂, CH₃CH)的过程. 电子密度拓扑分析显示在整个异构化过程中有两种类型的结构过渡态: 单态反应通道为T型过渡态, 三态反应通道为环状过渡态. 单线态RN脱氢反应通道中“原子-分子键”的存在说明两个H原子是以H₂的形式从RN中脱去的.     

Mechanism and rate constants of the CH2 + CH2CO reactions in triplet and singlet states: A theoretical study

Ab initio and density functional CCSD(T)-F12/cc-pVQZ-f12//B2PLYPD3/6-311G** calculations have been performed to unravel the reaction mechanism of triplet and singlet methylene CH₂ with ketene CH₂CO. The computed potential energy diagrams and molecular properties have been then utilized in Rice–Ramsperger–Kassel–Marcus-Master Equation (RRKM-ME) calculations of the reaction rate constants and product branching ratios combined with the use of nonadiabatic transition state theory for spin-forbidden triplet-singlet isomerization. The results indicate that the most important channels of the reaction of ketene with triplet methylene lead to the formation of the HCCO + CH₃ and C₂H₄ + CO products, where the former channel is preferable at higher temperatures from 1000 K and above. In the C₂H₄ + CO product pair, the ethylene molecule can be formed either adiabatically in the triplet electronic state or via triplet-singlet intersystem crossing in the singlet electronic state occurring in the vicinity of the CH₂COCH₂ intermediate or along the pathway of CO elimination from the initial CH₂CH₂CO complex. The predominant products of the reaction of ketene with singlet methylene have been shown to be C₂H₄ + CO. The formation of these products mostly proceeds via a well-skipping mechanism but at high pressures may to some extent involve collisional stabilization of the CH₃CHCO and cyclic CH₂COCH₂ intermediates followed by their thermal unimolecular decomposition. The calculated rate constants at different pressures from 0.01 to 100 atm have been fitted by the modified Arrhenius expressions in the temperature range of 300–3000 K, which are proposed for kinetic modeling of ketene reactions in combustion. © 2018 Wiley Periodicals, Inc.     

Geometry optimization of the ring-opened oxirane diradical: mechanism of the addition reaction of the triplet oxygen atom to olefins

Geometry optimizations of several low-lying diradical states of the ring-opened oxirane (·CH₂CH₂O·) were performed by using the energy gradients of the UHF MINDO/3, STO-3G and 4-31G solutions. Both the STO-3G and 4-3 IG methods predict that the most stable form is the triplet state of the non-twisted σπ conformation in which the unpaired spins localized on the terminal carbon and oxygen atoms are oriented perpendicularly to each other. The singlet σσ diradical state in which both the radical-site p orbitals are coplanar with the molecular framework is only 2.3 (STO-3G) and 1.2 (4-31G) kcal/mol less stable than the triplet σπ diradical state. It is found that the geometry of the singlet σσ diradical is unique in that the C-C-O angle is noticeably small as compared with various other diradical states. Implications of these results to the mechanism of the oxirane-forming O(³P) + C₂H₄ reaction are discussed.     

Ligand Effects on the Mechanisms of Thermal Bond Activation in the Gas‐Phase Reactions NiX+/CH4→Ni(CH3)+/HX (X=H,CH3, OH,F). Short Communication

The thermal ion‐molecule reactions NiX⁺+CH₄→Ni(CH₃)⁺+HX (X=H, CH₃, OH, F) have been studied by mass spectrometric methods, and the experimental data are complemented by density functional theory (DFT)‐based computations. With regard to mechanistic aspects, a rather coherent picture emerges such that, for none of the systems studied, oxidative addition/reductive elimination pathways are involved. Rather, the energetically most favored variant corresponds to a σ‐complex‐assisted metathesis (σ‐CAM). For X=H and CH₃, the ligand exchange follows a ‘two‐state reactivity (TSR)’ scenario such that, in the course of the thermal reaction, a twofold spin inversion, i.e., triplet→singlet→triplet, is involved. This TSR feature bypasses the energetically high‐lying transition state of the adiabatic ground‐state triplet surface. In contrast, for X=F, the exothermic ligand exchange proceeds adiabatically on the triplet ground state, and some arguments are proposed to account for the different behavior of NiX⁺/Ni(CH₃)⁺ (X=H, CH₃) vs. NiF⁺. While the couple Ni(OH)⁺/CH₄ does not undergo a thermal ligand switch, the DFT computations suggest a potential‐energy surface that is mechanistically comparable to the NiF⁺/CH₄ system. Obviously, the ligands X act as a mechanistic distributor to switch between single vs. two‐state reactivity patterns.     

Singlet-triplet separation in borylene. comparison with the methylidyne cation,methylene and the amino cation

Borylene (BH), which is known to have a ¹Σ⁺ ground state, is predicted to have a low-lying ³Π triplet state with an energy separation of 31.9 kcal/mol. The ab initio molecular-orbital method used is shown to give results which are consistent with theoretical and experimental data on CH⁺, CH₂ and NH⁺₂.     

Theoretical study on the gas‐phase reaction mechanism between nickel monoxide and methane for syngas production

The comprehensive mechanism survey on the gas‐phase reaction between nickel monoxide and methane for the formation of syngas, formaldehyde, methanol, water, and methyl radical has been investigated on the triplet and singlet state potential energy surfaces at the B3LYP/6‐311++G(3df, 3pd)//B3LYP/6‐311+G(2d, 2p) levels. The computation reveals that the singlet intermediate HNiOCH₃ is crucial for the syngas formation, whereas two kinds of important reaction intermediates, CH₃NiOH and HNiOCH₃, locate on the deep well, while CH₃NiOH is more energetically favorable than HNiOCH₃ on both the triplet and singlet states. The main products shall be syngas once HNiOCH₃ is created on the singlet state, whereas the main products shall be methyl radical if CH₃NiOH is formed on both singlet and triplet states. For the formation of syngas, the minimal energy reaction pathway (MERP) is more energetically preferable to start on the lowest excited singlet state other than on the ground triplet state. Among the MERP for the formation of syngas, the rate‐determining step (RDS) is the reaction step for the singlet intermediate HNiOCH₃ formation involving an oxidative addition of NiO molecule into the C? H bond of methane, with an energy barrier of 120.3 kJ mol^?1. The syngas formation would be more effective under higher temperature and photolysis reaction condition. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Theoretical study on the reaction mechanisms of CH2SH + NO reaction

The mechanisms for the CH₂SH + NO reaction were investigated on both of the singlet and triplet PES at the BMC-CCSD//B3LYP/6-311+G(d,p) level. The results indicate that the singlet PES is much lower than the triplet PES energetically; therefore, the reaction occurs on the singlet PES dominantly. The most favorable channel on the singlet PES takes place by a barrierless addition of N atom to CH₂SH radical to form HSCH₂NO. Subsequently, the rearrangement of the initial adduct HSCH₂NO (IM1) to form another intermediate IM3 via a four-center transition state, followed by the C–O bond fission in IM3 leading to the major product CH₂S + HNO. Due to high barriers, other product including HC(N)SH + HO, HON + CH₂S, and HNO + CHSH could be negligible. The direct abstraction channel was also determined to yield CH₂S + HON. With high barrier (33.3 kcal/mol), it is not competitive with the addition channel, in which all stationary points are lower than reactant energetically. While on the triplet PES, with the lowest barrier height (18.8 kcal/mol), the direct N-abstracted channel to form CH₂S + HNO is dominant. However, it is not competitive with the channels on the singlet PES. Our results are in good accordance with experimental conclusions that the reaction proceeds via addition mechanism.     

Theoretical study on the reaction mechanism of (CH<Subscript>3</Subscript>)<Subscript>3</Subscript>CO<Superscript>.</Superscript> radical with NO

The reaction mechanism of (CH₃)₃CO^. radical with NO is theoretically investigated at the B3LYP/6-31G* level. The results show that the reaction is multi-channel in the single state and triplet state. The potential energy surfaces of reaction paths in the single state are lower than that in the triple state. The balance reaction: (CH₃)₃CONO⇔(CH₃)₃CO^.+NO, whose potential energy surface is the lowest in all the reaction paths, makes the probability of measuring (CH₃)₃CO^. radical increase. So NO may be considered as a stabilizing reagent for the (CH₃)₃CO^. radical.     

Insight into the reaction mechanism of CH2SH with HO2: A density functional theory study

Density functional theory was adopted in this work to reveal the reaction mechanism of CH₂SH with HO₂. Reaction rate constants were computed from 200 to 2000 K using the transition state theory combined with Wigner and Eckart tunneling correction. Moreover, localized orbital locator, atoms in molecules and Mayer bond order analyses were used to study the essence of chemical bonding evolution. Eleven singlet paths and three triplet ones are located on the potential surface (PES). The results show that the main products on the singlet PES are ¹CH₂S and H₂O₂, whereas on the triplet PES they are CH₃SH + ³O₂, which are coincident with the similar reaction of CH₃S and HO₂. This conclusion is also supported by rate constant calculation results. Interestingly, all the possible paths are involved in the hydrogen transfer. The results have provided underlying insights to the analogous reactions and further experimental studies.     

A density functional theory study on the atmospheric reaction of CH3O2 with HS: Mechanism and kinetics

The reaction mechanism of CH₃O₂ and HS was systematically investigated by density functional theory (DFT). Six singlet pathways and seven triplet ones are located on the potential surface (PES). The result indicates that the main products are CH₃O and HSO both on the singlet and triplet PES, different from the CH₃O₂ + OH reaction. Moreover, deformation density (ρ_def) and atoms in molecules (AIM) analyses were carried out to further uncover the nature of chemical bonding evolution in the primary pathways. Furthermore, reaction rate constants were calculated in the temperature range from 200 to 1000 K using the transition state theory with the Wigner and Eckart tunneling corrections. Our results can shed light on the title reaction and offer instructions for analogous atmospheric reactions, as well as experimental research in the future.     

The DFT study on mechanisms for NCO with C2H5 reaction

A detailed investigation has been performed at the QCISD(T)/6‐311++G(d,p)//B3LYP/6‐311+G(d,p) level for the reaction of NCO with C₂H₅ by constructing singlet and triplet potential energy surfaces (PES). The results show that the title reaction is more favorable on the singlet PES than on the triplet PES. On the singlet PES, the initial addition processes are barrierless and release lots of energy. The dominant channel occurs via the fragmentations of the initial adduct C₂H₅NCO and C₂H₅OCN to form C₂H₄ + HNCO and HOCN, respectively. With higher barrier heights, other products such as CH₄ + HNC + CO, CH₃CHNH + CO, CH₃CH + HNCO, and CH₃CN + H₂ + CO are less competitive. On the triplet PES, the entrance reactions surpass significant barriers; therefore, it could be negligible at the normal atmospheric condition. However, the most feasible channel on the triplet PES is the direct hydrogen abstraction channel to form CH₂CH₂ + HNCO. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Theoretical study on the atmospheric reaction of CH3O2 with OH

A quantum chemical investigation on the reaction mechanism of CH₃O₂ with OH has been performed. Based on B3LYP and QCISD(T) calculations, seven possible singlet pathways and seven possible triplet pathways have been found. On the singlet potential energy surface (PES), the most favorable channel starts with a barrierless addition of O atom to CH₃O₂ leading to CH₃OOOH and then the O? O bond dissociates to give out CH₃O + HO₂. On the triplet PES, the calculations indicate that the dominant products should be ³CH₂O₂ + H₂O with an energy barrier of 29.95 kJ/mol. The results obtained in this work enrich the theoretical information of the title reaction and provide guidance for analogous atmospheric chemistry reactions. © 2015 Wiley Periodicals, Inc.     

Ion-Velocity Map Imaging Study of Photodissociation Dynamics of Acetaldehyde

The photodissociation dynamics of acetaldehyde in the radical channel CH3＋HCO has been reinvestigated using time-sliced velocity map imaging technique in the photolysis wavelength range of 275-321 nm. The CH3 fragments have been probed via （2＋1） resonance-enhanced multiphoton ionization. Images are measured for CH3 formed in the ground and excited states （v2=0 and 1） of the umbrella vibrational mode. For acetaldehyde dissociation on T1 state after intersystem crossing from S1 state, the products are formed with high translational energy release and low internal excitation. The rotational and vibrational energy of both fragments increases with increasing photodissociation energy. The triplet barrier height is estimated at 3.8814-0.006 eV above the ground state of acetaldehyde.     

Mechanistic and kinetic study the reaction of O(<Superscript>3</Superscript>P) + CH<Subscript>3</Subscript>CFCH<Subscript>2</Subscript>

The triplet potential energy surface of the O(³P) + CH₃CFCH₂ reaction has been investigated at the BMC-CCSD//MP2/6-311++G(d,p) level. Multichannel RRKM theory and transition state theory are employed to calculate the rate constants over a wide range of temperatures and pressures. The total rate constants show positive temperature dependence and pressure independence. At pressure of 10 Torr with He as bath gas, the addition/elimination on triplet potential energy surface is a dominant pathway. Major predicted end products are CH₃CFCHO (I) and H at the temperatures between 200 and 3,000 K; the direct hydrogen abstraction leading to OH + CH₂CFCH₂(I) plays an important role at higher temperatures. The calculated overall rate constants are in good agreement with the available experimental data.     

Theoretical study of model elementary reactions of dissociative addition of light hydrocarbons to the Ti-doped aluminide cluster Al<Subscript>12</Subscript>Ti

The potential energy surfaces (PES), energies E, and activation barriers h of elementary reactions of dissociative addition of CH₄ and C₂H₆ molecules to the Al₁₂Ti cluster with a marquee structure in the singlet and triplet states were calculated within the B3LYP approximation of the density functional theory using the 6-31G* basis set. The first stage of the reaction Al₁₂Ti + CH₄ leads to the adsorption complex CH₄ · Al₁₂Ti with the R(TiC) distance of ～2.4 Å. The methane molecule is coordinated as a tridentate ligand the singlet state and as a bidentate ligand in the triplet state, although both coordination modes are close in energy. In the transition state, the CH₄ molecule is coordinated through its active C-H bond to an inclined Ti-Al edge of the cluster, and the C-H bond is significantly elongated and weakened. The activation barrier height h referenced to the CH₄ complex is ～9 and ～19 kcal/mol for the singlet and triplet, respectively, and that referenced to the primary products Al₁₂Ti(CH₃)(H) is ～21 kcal/mol. The barrier to migration of the CH₃ group around the metal cluster is estimated at ～10 kcal/mol. At the initial stage of the reaction Al₁₂Ti + C₂H₆, two types of C₂H₆ · Al₁₂Ti adsorption complexes are formed. In one of them, the ethane molecule is coordinated through a methyl group (as the methane molecule); and in the other type, the coordination is through the C-C bond. This reaction can proceed through two paths by means of insertion into C-H or C-C bonds to give Al₁₂Ti(C₂H₅)(H) or Al₁₂Ti(CH₃)₂, respectively. The second path is impeded by a high barrier (～30 kcal/mol) and is possible, if at all, only at high temperatures. Conversely, the insertion into a C-H bond in ethane is somewhat more favorable than in methane. Analogously, the PES of addition of the second methane molecule to Al₁₂Ti(CH₃)(H) was calculated. The second molecule is adsorbed and dissociates by the same mechanism as the first CH₄ molecule, but with somewhat lower barriers and energy effect of formation of Al₁₂Ti(CH₃)₂(H)₂. The addition of propane and longer hydrocarbons is briefly considered. The results are compared with the results of previous analogous calculations of the PES of related reactions of dissociative adsorption of dihydrogen on the Al₁₂Ti cluster, which are more exothermic, have lower barriers, and can occur under milder conditions.     

Theoretical study on the atmospheric reaction of CH3SH with O2

A detailed study on the reaction mechanism of CH₃SH with O₂ was carried out using quantum chemical methods. Eleven singlet pathways and four triplet pathways were found based on CCSD(T)//M06-2x calculations. The nature of chemical bonding evolution was also studied using electron localization function and atoms in molecules analysis. Moreover, reaction rate constants were calculated between 200 and 800 K at the level of the transition state theory by Wigner tunneling correction. The results suggest that the main products should be CH₂SO, H₂O, CH₃OH, SO, CH₄, and SO₂, respectively, basically coinciding with the experimental results. The corresponding feasible pathways are channels R7, R8, and R9, respectively, with an effective energy barrier of 56.21 kJ/mol. Obviously, given the low energy barrier similar to the main paths mentioned above, the products CH₂SH and HO₂ should assume a definite proportion in all possible products, although such species were not yet detected in experiment.     

CH2CO＋CH2的多通道反应的理论研究

利用密度泛函(DFT)和自然键轨道理论(NBO)及高级电子耦合簇[CCSD(T)]和电子密度拓扑(AIM)方法, 对单重态和三重态CH₂与CH₂CO反应的微观机理进行了研究. 在B3LYP/6-311＋G(d,p)水平上优化了反应通道各驻点的几何构型. 在CCSD(T)/6-311＋G(d,p)水平上计算了各物种的单点能量, 并对总能量进行了校正. 计算表明, 单重态CH₂与CH₂CO的C—H键可发生插入反应, 与C＝C、C＝O可发生加成反应, 存在三条反应通道, 产物为CO和C₂H₄, 从能量变化和反应速控步骤能垒两方面考虑, 反应II更容易发生. 对反应通道中的关键点进行了自然键轨道及电子密度拓扑分析. 三重态CH₂与CH₂CO的反应存在三条反应通道, 一条是与C-H键的插入反应, 另一条是三重态CH₂与C＝C发生加成反应, 产物为CO和三重态C₂H₄, 通道II势垒较低, 更容易发生. 最后一条涉及双自由基的反应活化能最大, 最难发生.