Ab initio quantum chemical studies on the reactions of CF3O2 with OH

The potential energy surface for the CF₃O₂ + OH reaction has been theoretically investigated using the DFT (B3LYP/6-311G(d,p)) level of theory. Both singlet and triplet potential energy surfaces are investigated. The reaction mechanism on the triplet surface is simple. However, the reaction mechanism on the singlet surface is more complicated. It is revealed that the formation of CF₃O + HO₂ is the dominant channel on the triplet surface. The potential energy surface (PES) for this reaction has been given according to the relative energies calculated at the DFT/B3LYP/6-311G(d,p) level. Because this reaction involves both triplet and singlet states, triplet–singlet intersystem crossing (ISC) crossing also have been investigated in this paper.     

Theoretical Studies on Reaction Mechanism of CH3SO with NO

The potential energy surface （PES） of CH3SO radical with NO reaction has been studied at MP2/6-311G（2df, p） and QCISD/6-311G（2df, p） levels. Geometries of the reactants, transition states （TS） and products were optimized at B3LYP/6-311G （d,p） level. The geometries of the transition states were found for the first time. The calculated results show that the reaction can proceed via singlet-state or triplet-state PES. Because of the high energy barrier of triplet surface, the singlet surface reactions are dominant. The topological analysis of electron density shows that there are two kinds of structaral transition states （the bifurcation-type ring structure transition state and the T-shaped conflict structure transition state） in the titled reaction. The total electronic density of the reactants, TS and products and the spin electronic density on the triplet surface were also discussed in this paper.     

High‐level ab initio study of the N+(3P)+SH2 reactions in the gas phase: Role of spin‐forbidden pathways

The singlet and triplet potential energy surfaces involved in N⁺+SH₂ reactions have been explored using high‐level ab initio techniques. The geometries of the stationary points were optimized at the QCISD/6‐311G(df,p) level. The final energies were obtained in CCSD(T)/6‐311+G(3df,2p) single‐point calculations. The results obtained show that, although the N⁺(¹D)+SH₂ entrance channel is higher in energy than the N⁺(³P)+SH₂ one, most of the [H₂, S, N]⁺ singlet state cations are lower in energy than the corresponding triplets, due to their different bonding characteristics. Both singlet and triplet potential energy surfaces are quite close each other, and crossover between them can occur. The minimum energy crossing points were located by means of CASSCF(6,5) calculations. The spin‐orbit couplings show that the transition probability from the triplet to the singlet potential energy surface is significantly large. One of the most important consequences is that some of the products of the reaction, such as SH⁺, can be formed in typical spin‐forbidden processes. Since all the relevant structures along these pathways are much lower in energy than the reactants, this mechanism should be accessible even at low impact energies and therefore could be important in processes taking place in interstellar media. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Theoretical study of [Si,O,C,O] species: Prediction of new species on triplet potential energy surface

The intermediates [Si,O,C,O] of the Si + CO₂ reaction have been studied in detail using high level ab iniitio methods. Both singlet and triplet [Si,O,C,O] species are characterized structurally and energetically. On the singlet potential energy surface (PES), the vdw‐OSi–CO isomer and in the triplet PES, the bent‐SiOCO isomer is found to be thermodynamically as well as kinetically most stable species. All possible isomerization transition states (TS) are located on both singlet and triplet potential surfaces. On the triplet surface, the stability of the trans‐OSiCO isomer is comparable with that of the bent‐SiOCO isomer. A non‐planar cis‐SiOCO isomer is located on the triplet PES, which is predicted for the first time. Heats of formation at 0 K (Δ_fH°, 0 K) for all singlet and triplet species are computed using G3B3, G3MP2, and CBS‐Q theories. The discrepancy between G3B3 and the other two methods for the heat of formation value for triplet trans‐OSiCO is discussed. The PESs for singlet as well as triplet species with their dissociation asymptotes are explored at the CCSD(T)/6‐311G(d,p)//MP2/6‐311G(d,p) level of theory. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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     

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中脱去的.     

Structure and Stability of Interstellar Molecule C_3S

ＩｎｔｒｏｄｕｃｔｉｏｎＳｍａｌｌｃｌｕｓｔｅｒｓｃｏｎｔａｉｎｉｎｇｃａｒｂｏｎａｎｄｓｕｌｆｕｒ,ｓｕｃｈａｓＣＳ ,Ｃ2 ＳａｎｄＣ3 Ｓ ,ｗｈｉｃｈｐｏｓｓｅｓｓｌａｒｇｅｐｅｒｍａｎｅｎｔｄｉｐｏｌｅｍｏｍｅｎｔｓａｎｄｈａｖｅｂｅｅｎｉｄｅｎｔｉｆｉｅｄｉｎｔｈｅｃａｒｂｏｎｓｔａｒＩＲＣ+ 10°2 16ａｎｄｉｎｔｈｅＴａｕｒｕｓｃｏｌｄｍｏｌｅｃｕｌａｒｄｅｎｓｅｃｌｏｕｄＴＭＣ 1,1 7ｈａｖｅａｔｔｒａｃｔｅｄｍｕｃｈａｔｔｅｎｔｉｏｎｂｅｃａｕｓｅｏｆｔｈｅｉｒｉｍｐｏｒｔａｎｔｒｏｌｅｓ…     

Theoretical study on the potential energy surface of NC<Subscript>3</Subscript>P isomers

Density functional theory (DFT) calculations have been used to study the isomerization process in the NC₃P system. At the DFT/B3LYP/6-311G(d) level, 28 triplet and 28 singlet minima were obtained on their respective potential energy surfaces. The linear triplet ³NCCCP is the lowest-energy structure among the isomers. On the triplet PES, only linear isomers ³NCCCP, ³CNCCP, ³CCCNP, and ³CCNCP possess great kinetic and thermodynamic stabilities to exist under low-temperature conditions (such as in the dense interstellar clouds). At the same time, one chain-like and four three-membered-ring isomers on the singlet PES have been located with high kinetic and thermodynamic stabilities. Further CCSD(T)/6-311G(2df)//QCISD/6-311G(d), CCSD(T)/cc-pVTZ//DFT/B3LYP/cc-pVTZ, and CASPT2(14,12)/cc-pVQZ//CASSCF(14,12)/cc-p VQZ calculations are performed on the structures, frequencies, and energies of the relevant species. The bonding natures were analyzed and the results were compared with the analogous NC₃N and NC₂P molecules so as to aid their future experimental or astrophysical detection.     

CH_3SO+HO_2气相反应机理与主通道速率常数

在MPW1PW91/6-311G(d,p)水平上优化了标题反应各驻点物种的几何构型,并在相同水平上通过频率计算和内禀反应坐标(IRC)分析对过渡态结构及连接性进行了验证.采用QCISD(T)/6-311G(d,p)方法对所有驻点及反应路径的部分选择点进行单点能校正,分别构建了CH3SO+HO2反应体系的单、三重态反应势能剖面.研究结果表明,CH3SO+HO2反应体系存在6条反应通道7条路径,优势通道(1)R→3IM→P1(CH3SOH+3O2)发生在三重态势能面上,此通道包含两条路径,其表观活化能分别为12.01和-30.04kJ?mol-1,主路径(2)R→3IM→3TS2→P1(CH3SOH+3O2)是一个无势垒氢迁移过程.利用经典过渡态理论(TST)与变分过渡态理论(CVT)并结合小曲率隧道效应模型(SCT),分别计算了主路径(2)在200～2500K温度范围内的速率常数kTST,kCVT和kCVT/SCT,在此温度区间内的表观反应速率常数三参数表达式为kCVT/SCT=4.08×10-24T3.13exp(8012.2/T)cm3imolecule-1is-1,具有负温度系数效应.速率常数计算结果显示,变分效应在计算温度段内影响较小,而量子力学隧道效应在低温段有显著影响.     

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     

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势垒较低, 更容易发生. 最后一条涉及双自由基的反应活化能最大, 最难发生.     

Ab initio quantum chemical studies of the reactions of CF3CFHO2 with HO2

The potential energy surface (PES) for the CF₃CFHO₂+HO₂ reaction has been theoretically investigated using the DFT [B3LYP/6‐311G(d,p)] and B3LYP/6‐311++G(3df,3pd)//B3LYP/6‐311G(d,p) levels of theory. Both singlet and triplet PESs are investigated. The reaction mechanism on the triplet surface is simple. It is revealed that the formation of CF₃CFHOOH+³O₂ is the dominant channel on the triplet surface. On the basis of the ab initio data, the total rate constants for the reaction CF₃CFHO₂+HO₂ in the T = 210–500 K range have been computed using conventional transition state theory with Wigner's tunneling correction and have been fitted by a rate constant expression as k = 1.04 ×10^?12(cm³ molecule^?1 s^?1) exp (700.33/T). Calculated transition state rate constants with Wigner's tunneling correction for the reaction CF₃CFHO₂+HO₂ are in good agreement with the available experimental values. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Mechanism and rate constants for ammonia photochemical oxidation with O2 on the singlet and triplet potential energy surfaces

The mechanism and thermodynamic of NH₃ + O₂ reaction on the singlet and triplet potential energy surfaces (PES), were carried out using the RMP2 and CCSD (T)//RMP2 theoretical approaches in connection with the 6-311++G(d, p) basis set. Three pre-reactive complexes, ¹C1, ¹C2, and ³C1 on the singlet and triplet PES were formed between ammonia and molecular oxygen. With variety of pre-reactive complexes, six types of products are obtained, of which two types are found to be thermodynamically stable. The mechanistic properties of all products channels are discussed. Results show that production of HONO + H₂ and HN(OH)₂ are the main reaction channels in thermodynamic viewpoint with the Gibbs free energy of ?G° = ?34.681 and ?27.153 kcal/mol, respectively. Rate constants of the title reaction over the temperature range of (200–1000 K) show kinetic products are different from thermodynamic products.     

Computational mechanistic investigation of the gas phase C2H4 + CO reaction on the singlet and triplet potential energy surfaces

Reaction pathways of ethylene and carbon monoxide on the singlet and triplet potential energy surfaces (PESs) have been calculated at B3LYP/6-311++G (3df, 3dp), G3B3 and CCSD(T)//B3LYP levels. Reaction mechanisms have been investigated by analysis of various structures. Suggested reaction mechanisms reveal that ³P3(CH₂CHCHO) and ³P4(CH₃CCHO) are thermodynamically stable adducts with the negative value in Gibbs free energies on the triplet PES. In addition, results show that one intersystem crossing exists between triplet and singlet PESs, which are obtained by scanning of the C–C bond length in ¹IN3 and ³IN7 species.     

Exploring the Reaction Paths on the Potential Energy Surfaces of the S1 and T1 States in Methylenecyclopropane

The reaction paths of methylenecyclopropane 1 on the potential energy surfaces (PESs) of the lowest triplet (T₁) state and the lowest excited singlet (S₁) state, as well as that of the ground state (S₀), were explored by using the nudged elastic band method at the MRMP2//MCSCF/6‐31++G(d,p) and DFT(B3LYP)/6‐31++G(d,p) levels of theory. After vertical excitation of 1, three transition states on the PES of the lowest triplet state and one transition state on the S₁ PES were found along the reaction path to produce a carbene, cyclobutylidene 2. All of these transition states are lower in energy than the S₁ state produced by vertical excitation at the S₀ energy minimum in 1. Fast transition is predicted to occur from the T₁ state or from the S₁ state to the S₀ state due to strong spin‐orbit coupling or nonadiabatic coupling in the geometrical vicinity of 2. On the MRMP2 S₀ PES, the energy barriers of 5.0, 10.3 and 13.5 kcal mol^?1 were obtained for C migration reaction (backward reaction), 1,2‐H migration reaction to cyclobutene 3, and 1,3‐H migration reaction to bicyclopropane 4, respectively, started at 2. The introduction of phenyl groups makes the energy barriers smaller due to the π conjugation between the carbene center and phenyl groups.     

Ab Initio Study on the Potential Energy Surfaces of B4 Cluster

The singlet and triplet potential energy surfaces (PES) for the isomerization and dissociation reactions of B₄ isomers have been investigated using ab initio methods. Ten B4 isomers have been identified and of these 10 species, 4 have not been reported previously. The singlet rhombic structure ¹1 is found to be the most stable on the B₄ surface, in agreement with the results of previous reports. Several isomerization and dissociation pathways have been found. On the singlet PES, the linear ¹3b can rearrange to rhombus ¹1 directly, while ¹3c rearranges to ¹1 through two‐step reactions involving a cyclic intermediate. On the triplet PES, the capped triangle structure ³2 undergoes ring opening to the linear isomer ³3b with a barrier of 34.8 kcal/mol and 44.9 kcal/mol, and the latter undergoes ring closure to the square structure ³1 with a barrier of 30.4 kcal/mol and 33.0 kcal/mol at the MP4/6–311+G(3df)//MP2/6–311G(d) and CCSD/aug‐cc‐pVTZ//MP2/6–311G(d) levels of theory, respectively. The direct decomposition of singlet B₄ yielding to B₃+B is shown to have a large endothermicity of 87.3 kcal/mol (CCSD), and that producing 2B₂ to have activation energy of 133.4 kcal/mol (CCSD).     

Computational study on the reaction mechanism of the gas-phase atom-negative ion of S?+?NO2 ?: comparative study of mechanism with S?+?O3 reaction as isoelectronic and isostructure systems

The reaction mechanism of sulfur vapor (S) with nitrite ion (NO₂ ⁻) has been investigated theoretically on the triplet and singlet potential energy surfaces (PESs). All stationary points for the title reaction have been optimized at the B3LYP/6-311+G(3df) level. The energetic data have been obtained at the CCSD(T)//B3LYP level employing the 6-311+G(3df) basis set. Five stable collision complexes, ³IN1 (S–ONO⁻), ³IN2 (cyclic SONO⁻), ¹IN1 (cis S–ONO⁻), ¹IN2 (S–NO₂ ⁻), and ¹IN3 (trans S–ONO⁻), have been considered on the triplet and singlet PESs through barrier-less and exothermic processes. By starting from these complexes, a simple mechanism has been obtained on the triplet PES while a complex mechanism has been considered on the singlet PES. The calculated results show that there are no favorable paths for the reaction of S with NO₂ ⁻ on the singlet PES. Therefore, the S + NO₂ ⁻ reaction proceeds only on the triplet PES to produce ³SO + ³NO⁻ as main products. The results from the comparative study of S + NO₂ ⁻ reaction mechanism with S + O₃ (as isoelectronic and isostructure reactions) on the singlet PES show similarities in the overall trend of reaction mechanism and atom connectivity and differences in the stability of intermediates and the energy barriers of transition states.     

Theoretical study on structures and stability of HC2P isomers

 The structures and isomerization pathways of various HC₂P isomers in both singlet and triplet states are investigated at the B3LYP/6-311G(d,p), QCISD/6-311G(d,p) (for isomers only) and single-point CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p) levels. At the CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p) level, the lowest-lying isomer is a linear HCCP structure ³ 1 in the ³∑⁻ state. The second low-lying isomer has a CPC ring with exocyclic CH bonding ¹ 5 in a singlet state at 10.5 kcal/mol. The following third and fourth low-lying isomers are a singlet bent HCCP structure ¹ 1 at 20.9 kcal/mol and a bent singlet HPCC structure ¹ 3 at 35.8 kcal/mol, respectively. Investigation of the HC₂P potential-energy surface indicates that in addition to the experimentally known isomer ³ 1, the other isomers ¹ 1, ¹ 3 and ¹ 5 also have considerable kinetic stability and may thus be observable. However, the singlet and triplet bent isomers HCPC ¹ 2 and ³ 2 as well as the triplet bent isomer HPCC ³ 3 are not only high-lying but are also kinetically unstable, in sharp contrast to the situation of the analogous HCNC and HNCC species that are both kinetically stable and that have been observed experimentally. Furthermore, the reactivity of various HC₂P isomers towards oxygen atoms is briefly discussed. The results presented here may be useful for future identification of the completely unknown yet kinetically stable HC₂P isomers ¹ 1, ¹ 3 and ¹ 5 either in the laboratory or in interstellar space. Received: 5 November 2000 / Accepted: 25 November 2001 / Published online: 8 April 2002     

Quantum chemical study on insertion and abstraction reaction of dichlorocarbene with methyl alcohol and methyl mercaptan

The insertion and abstraction reaction mechanisms of singlet and triplet CCl₂ with CH₃MH (M=O, S) have been studied by using the DFT, NBO and AIM methods. The geometries of reactions, the transition state and products were completely optimized by B3LYP/6–311G(d, p). All the energy of the species was obtained at the CCSD(T)/6–311G(d, p) level. The calculated results indicated that the major pathways of the reaction were obtained on the singlet potential energy surface. The singlet CCl₂ can not only trigger the insertion reaction with C-H and M-H in four pathways, by which the products P1 [CH₃OCHCl₂, reaction I(1)], P3[Cl₂HCCH₂OH, reaction I(2)], P5[CH₃SCHCl₂, reaction II(1)] and P7[Cl₂HCCH₂SH, reaction II(2)] are produced respectively, but also abstract M-H, resulting P4 [CH₂O+CH₂Cl₂, reaction I(3)] and P8[CH₂S+CH₂Cl₂, reaction II(3)]. In addition, the important geometries in domain pathways have been studied by AIM and NBO theories. Supported by the National Natural Science Foundation of China (Grant No. 20335030) and Foundation of Education Committee of Gansu Province (Grant No. 0708-11)     

A theoretical study of H?H σ bond activation catalyzed by VO2+ in gas phase

The mechanism of H? H σ bond activation catalyzed by VO(¹A₁/³A′) has been investigated by using density functional theory at the B3LYP/6‐311G(2d, p) level and the single‐point energy calculations were done at the CCSD/6‐311G (2d, p)//B3LYP/6‐311G(2d, p) level of theory using the geometries along the minimum energy pathway. According to our calculation results, the different reaction mechanisms were found for the singlet and triplet potential energy surfaces (PESs). Specially, the crossing points (CPs) between the different PESs have been located by means of the intrinsic reaction coordinate approach used by Yoshizawa et al, and corresponding minimum energy CPs that we obtained by the mathematical algorithm proposed by Harvey et al. has also been employed. In addition, the orbital interaction for ion‐molecule complexes ¹IM1 and ³IM1 have been examined by fragment molecular orbital analysis. Finally, the frontier molecular orbital interaction analysis about ³TS1 and ³TS2 were used to gain useful information about the H? H σ bond activation by VO. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012