Theoretical study of the mechanism of CH2CO + CN reaction

The potential energy surface information of the CH₂CO + CN reaction is obtained at the B3LYP/6‐311+G(d,p) level. To gain further mechanistic knowledge, higher‐level single‐point calculations for the stationary points are performed at the QCISD(T)/6‐311++G(d,p) level. The CH₂CO + CN reaction proceeds through four possible mechanisms: direct hydrogen abstraction, olefinic carbon addition–elimination, carbonyl carbon addition–elimination, and side oxygen addition–elimination. Our calculations demonstrate that R→IM1→TS3→P3: CH₂CN + CO is the energetically favorable channel; however, channel R→IM2→TS4→P4: CH₂NC + CO is considerably competitive, especially as the temperature increases (R, IM, TS, and P represent reactant, intermediate, transition state, and product, respectively). The present study may be helpful in probing the mechanism of the CH₂CO + CN reaction. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

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     

A dual-level ab initio and hybrid density functional theory dynamics study on the unimolecular decomposition reaction C2H5O-->CH2O + CH3

We present a direct ab initio and hybrid density functional theory dynamics study of the thermal rate constants of the unimolecular decomposition reaction of C2H5O-->CH2O + CH3 at a high-pressure limit. MPW1K/6-31+G(d,p), MP2/6-31+G(d,p), and MP2(full)/6-31G(d) methods were employed to optimize the geometries of all stationary points and to calculate the minimum energy path (MEP). The energies of all the stationary points were refined at a series of multicoefficient and multilevel methods. Among all methods, the QCISD(T)/aug-cc-pVTZ energies are in good agreement with the available experimental data. The rate constants were evaluated based on the energetics from the QCISD(T)/aug-cc-pVTZ//MPW1K/6-31+G(d,p) level of theory using both microcanonical variational transition state theory (microVT) and RRKM theory with the Eckart tunneling correction in the temperature range of 300-2500 K. The calculated rate constants at the QCISD(T)/aug-cc-pVTZ/MPW1K/6-31+G(d,p) level of theory are in good consistent with experimental data. The fitted three-parameter Arrhenius expression from the microVT/Eckart rate constants in the temperature range 200-2500 K is k = 2.52 x 10(12)T(0.41)e(-8894.0/T) s(-1). The falloff curves of pressure-dependent rate constants are performed using master-equation method within the temperature range of 391-471 K. The calculated results are in good agreement with the available experimental data.     

DFT studies on the mechanism of the reaction of C2H5S with NO2

The mechanisms for the reaction of C₂H₅S with NO₂ are investigated at the QCISD(T)/6‐311++G(d, p)//B3LYP/6‐311++G(d, p) level on both single and triple potential energy surfaces. The geometries, vibrational frequencies and zero‐point energy (ZPE) corrections of all stationary points involved in the title reaction are calculated at the B3LYP/6‐311++G(d, p) level. The results show that the reaction is more predominant on the single potential energy surface, while it is negligible on the triple potential energy surface. Without barrier height in the whole process, the major channel is R → C₂H₅SONO (IM1 and IM2) → P1 (C₂H₅SO+NO). With much heat released in the formation of C₂H₅SNO₂ (IM3) and the transition state involved in the subsequent step more stable than reactants, P4 (CH₃CHS + t‐HONO) is subdominant product energetically. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Mechanistic and kinetic investigations of N2H4 + OH reaction

The reaction of N₂H₄ with OH has been investigated by quantum chemical methods. The results show that hydrogen abstraction mechanism is more feasible than substitution mechanism thermodynamically. The calculated rate constants agree with the available experimental data. The calculated results show that the variational effect is small at lower temperature region, while it becomes significant at higher temperature region. On the other hand, the small‐curvature tunneling effect may play an important role in the temperature range 220?3000 K. Moreover, the calculated rate constants show negative temperature dependence at the temperatures below 500 K, which is in accordance with Vaghjiani's report that slightly negative temperature dependence is found over the temperature range of 258?637 K. The mechanism of the major product (N₂H₃) with OH has also been investigated theoretically to understand the title reaction thoroughly. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Direct dynamic study on the hydrogen abstraction reaction CH3CN + OH-->CH2CN + H2O

The reaction of acetonitrile with hydroxyl has been studied using the direct ab initio dynamics methods. The geometries, vibrational frequencies of the stationary points, as well as the minimum energy paths were computed at the BHandHLYP and MP2 levels of theory with the 6-311G(d, p) basis set. The energies were further refined at the PMP4/6-311+G(2df, 2pd) and QCISD(T)/6-311+G(2df, 2pd) levels of theory based on the structures optimized at BHandHLYP/6-311G(d, p) and MP2/6-311G(d, p) levels of theory. The Polyrate 8.2 program was employed to predict the thermal rate constants using the canonical variational transition state theory incorporating a small-curvature tunneling correction. The computed rate constants are in good agreement with the available experimental data.     

Theoretical study on the rate constants for the C2H5 + HBr --> C2H6 + Br reaction

The reaction C(2)H(5) + HBr --> C(2)H(6) + Br has been theoretically studied over the temperature range from 200 to 1400 K. The electronic structure information is calculated at the BHLYP/6-311+G(d,p) and QCISD/6-31+G(d) levels. With the aid of intrinsic reaction coordinate theory, the minimum energy paths (MEPs) are obtained at the both levels, and the energies along the MEP are further refined by performing the single-point calculations at the PMP4(SDTQ)/6-311+G(3df,2p)//BHLYP and QCISD(T)/6-311++G(2df,2pd)//QCISD levels. The calculated ICVT/SCT rate constants are in good agreement with available experimental values, and the calculate results further indicate that the C(2)H(5) + HBr reaction has negative temperature dependence at T < 850 K, but clearly shows positive temperature dependence at T > 850 K. The current work predicts that the kinetic isotope effect for the title reaction is inverse in the temperature range from 200 to 482 K, i.e., k(HBr)/k(DBr) < 1.     

Theoretical study on the OH + CH3NHCOOCH3 reaction

The multiple-channel reactions OH + CH3NHC(O)OCH3 --> products are investigated by direct dynamics method. The optimized geometries, frequencies, and minimum energy path are all obtained at the MP2/6-311+G(d,p) level, and energetic information is further refined by the BMC-CCSD (single-point) method. The rate constants for every reaction channels, R1, R2, R3, and R4, are calculated by canonical variational transition state theory with small-curvature tunneling correction over the temperature range 200-1000 K. The total rate constants are in good agreement with the available experimental data and the two-parameter expression k(T) = 3.95 x 10(-12) exp(15.41/T) cm3 molecule(-1) s(-1) over the temperature range 200-1000 K is given. Our calculations indicate that hydrogen abstraction channels R1 and R2 are the major channels due to the smaller barrier height among four channels considered, and the other two channels to yield CH3NC(O)OCH3 + H2O and CH3NHC(O)(OH)OCH3 + H2O are minor channels over the whole temperature range.     

Theoretical study of the reaction of CN with C2H2+

 A theoretical study of the reaction of CN with C₂H₂ ⁺ has been carried out at three levels of theory, namely G2, B3LYP and CCSD(T). The main conclusion is that this is a feasible process under interstellar conditions, but only linear species may be produced. The most favourable product is HCCCN⁺, followed by CCCNH⁺. Production of HCCNC⁺ is predicted to be slightly endothermic; therefore, the reaction of CN + C₂H₂ ⁺ may produce precursors of HC₃N and C₃N in space. Furthermore, the B3LYP level is found to perform rather well compared with G2 and even better than CCSD(T). Received: 14 September 1999 / Accepted: 3 February 2000 / Published online: 12 May 2000     

CH_3自由基与C_2H_5CN高温反应机理的理论研究

用G3(MP2)//B3PW91/6-311G(d,p)双级别方法研究了CH_3自由基与C_2H5_CN的反应机理和动力学性质.计算表明反应存在抽氢、加成-消除和取代3种机理7条反应通道.用CVT方法计算了所有反应通道在1 000K~3 000K温度范围内的速率常数,结果表明计算值与实验值符合得很好.     

CH_3SH+CN微观反应机理的理论研究

采用BMC-CCSD//B3LYP/6-311G(d,p)方法对CH3SH+CN反应机理进行了详细的理论研究.反应中涉及的各稳定点的构型、振动频率和零点能在B3LYP/6-311G(d,p)水平下计算得到,计算结果表明,该反应存在两种反应机理,5条可能的反应通道.SN2机理由于能垒太高,与直接氢抽提机理相比可以忽略.该反应的最可行通道为CN中的C原子进攻SH中的H原子经由一个前期和一个后期分子络合物生成产物CH3S和HCN.计算得到的反应焓变与已有实验值非常吻合.     

Ab initio MO study of potential energy surface of NH2 with CN reaction

An extensive quantum chemical study of the potential energy surfaces (PES) for the association reaction of NH₂ with CN and the subsequent isomerization and dissociation reactions has been carried out using density functional theory (DFT)/B3LYP/6‐311++G(3df,2p) level of theory on both singlet and triplet states. The reaction mechanism on the triplet surface is more complicated than that on the singlet surface. A total of 19 isomers and 46 transition states have been identified and characterized on the triplet PES. Among them, IM2 (IM2a), IM3 (IM3a, IM3b), and IM10 are the lowest‐lying isomers with thermodynamic stability. Twenty available dissociation channels, depending on the different initial isomers, have been identified. On the singlet surface, only 12 isomers and 16 transition states have been found, and among them IM1(S) and IM2(S) are the lowest‐lying isomers. The higher isomerization and dissociation barriers on the singlet surface indicate that the addition and the subsequent reactions of NH₂+CN are most likely to occur on the triplet PES because of the lower barriers. A prediction can be made for the possible mechanism explaining the production of H+HNCN. Besides HNCN, other major products are NH+HCN and NH+HNC, which are produced by direct dissociation reactions from triplet IM2 and IM3, respectively. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

自由基HO2与C2H2反应势能面的理论研究

应用量子化学从头计算和密度泛函理论(DFT)对HO2+C2H2反应体系的反应机理进行了研究.在B3LYP/6-311G**和CCSD(T)/6-311G**水平上计算了HO2+ C2H2反应的二重态反应势能面.计算结果表明,主要反应方式为自由基HO2的H原子和C2H2分子中的C原子结合,经过一系列异构化,最后分解得到主要产物P1 (CH2O+ HCO).此反应是放热反应,化学反应热为-321.99 kJ·mol-1.次要产物为P2 (CO2 +CH3),也是放热反应.     

Theoretical study on the Br + CH3SCH3 reaction

The multiple-channel reactions Br + CH(3)SCH(3) --> products are investigated by direct dynamics method. The optimized geometries, frequencies, and minimum energy path are all obtained at the MP2/6-31+G(d,p) level, and energetic information is further refined by the G3(MP2) (single-point) theory. The rate constants for every reaction channels, Br + CH(3)SCH(3) --> CH(3)SCH(2) + HBr (R1), Br + CH(3)SCH(3) --> CH(3)SBr + CH(3) (R2), and Br + CH(3)SCH(3) -->CH(3)S + CH(3)Br (R3), are calculated by canonical variational transition state theory with small-curvature tunneling correction over the temperature range 200-3000 K. The total rate constants are in good agreement with the available experimental data, and the two-parameter expression k(T) = 2.68 x 10(-12) exp(-1235.24/T) cm(3)/(molecule s) over the temperature range 200-3000 K is given. Our calculations indicate that hydrogen abstraction channel is the major channel due to the smallest barrier height among three channels considered, and the other two channels to yield CH(3)SBr + CH(3) and CH(3)S + CH(3)Br are minor channels over the whole temperature range.     

星际媒介H2NCH2CN与H2O反应中的H2O分子偶合质子转移机理和氢隧道效应

H2NCH2CN+H2O→H2NCH2C(OH)NH是一个重要的反应, 涉及到星际媒介中甘氨酸的形成, 与早期地球上的氨基酸起源有关. 如果没有考虑氢隧道效应, 在MP2/6-311+G(d,p)级别上计算反应能垒是254.7 kJ·mol-1, 在星际媒介中该气相反应很难进行. 在星际媒介冰颗粒表面上, 水分子催化反应增强了该化学反应的活性. H2NCH2CN与(H2O)3反应中的两个水分子作为催化剂降低活化能77.5 kJ·mol-1和活化自由能70.9 kJ·mol-1, 并且通过氢键桥协同传递质子. 量子氢隧道对于该反应进行至关紧要,采用小弯曲隧道(SCT)近似和正则变分过渡态理论(CVT)方法研究. 温度50 K时, 速率常数kSCT/CVT为1.86×10-23 cm3·molecule-1·s-1, 表明在星际媒介中通过质子隧道机理该反应容易进行. 研究结果与地球上的氨基酸起源于地球本身物质的观点相一致.     

Ab initio Study of the Potential Energy Surface and Product Branching Ratios for Reaction of O（^1D） with C2H5CI

The potential energy surface of O（^1D）＋C2H5Cl reaction was studied using QCISD（T）/6- 311＋＋G（d,p）//MP2/6-31G（d,p） method. The calculations reveal an insertion-elimination mechanism. The insertion reaction of O（^1D） and C2H5Cl produces two energy-rich intermediates, IM1 and IM2, which subsequently decompose into various products. The calculations of the branching ratios of various products formed through the two intermediates were carried out using RRKM （Rice-Ramsperger-Kassel-Marcus） theory at the collision energies of 0, 20.9, 41.8, 62.7, 83.6, 104.5, and 125.4 kJ/mol. HCl is the main decomposition product for IM1; CH2OH is the main decomposition product for IM2. Since IM1 is more stable than IM2, HCl is probably the main product of the O（^1D）＋C2H5Cl reaction.     

Assignment on the X,A, B,C, and D states of the C6H5Br+ cation based on high‐level calculations

Geometries, frequencies, and energies of the 1²B₁, 1²A₂, 1²B₂, 2²B₁, 2²B₂, and 1²A₁, of the C₆H₅Br⁺ ion were calculated by using CASSCF and CASPT2 methods in conjunction with an ANO‐RCC basis. The CASPT2//CASSCF adiabatic excitation energies and CASPT2 relative energies for the six states are in good agreement with experiment. The X, A, B, C, and D electronic states of the C₆H₅Br⁺ ion were assigned to be X²B₁, A²A₂, B²B₂, C²B₁, and D²B₂ based on the CASSCF and CASPT2 calculations. The assignment on the D state of the C₆H₅Br⁺ ion is different from the previously published works. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Theoretical study and rate constant calculation for reaction of CF(3)CH(2)OH with OH

The reaction mechanism of CF(3)CH(2)OH with OH is investigated theoretically and the rate constants are calculated by direct dynamics method. The potential energy surface (PES) information, which is necessary for dynamics calculation, is obtained at the B3LYP/6-311G (d, p) level. The single-point energy calculations are performed at the MC-QCISD level using the B3LYP geometries. Complexes, with the energies being less than corresponding reactants and products, are found at the entrance and exit channels for methylene-H-abstraction channel, while for the hydroxyl-H-abstraction channel only entrance complex is located. By means of isodesmic reactions, the enthalpies of the formation for the species CF(3)CH(2)OH, CF(3)CHOH, and CF(3)CH(2)O are estimated at the MC-QCISD//B3LYP/6-311G (d, p) level of theory. The rate constants for two kinds of H-abstraction channels are evaluated by canonical variational transition state theory with the small-curvature tunneling correction (CVT/SCT) over a wide range of temperature 200-2000 K. The calculated results are in good agreement with the experimental values in the temperature region 250-430 K. The present results indicate that the two channels are competitive. Below 289 K, hydroxyl-H-abstraction channel has more contribution to the total rate constants than methylene-H-abstraction channel, while above 289 K, methylene-H-abstraction channel becomes more important and then becomes the major reaction channel.     

Theoretical Studies on Structures and Stabilities of C4H2+ Isomers

The structures, energies, stabilities and spectroscopies of doublet C₄H₂⁺ cations were explored at the DFT/B3LYP/6-311G(d,p), CCSD(T)/6-311+G(2df,2pd)(single-point), and G3B3 levels. Ten minimum isomers including the chainlike, three-member-ring, and four-member-ring structures are interconverted by means of 15 interconversion transition states. The potential energy surface was investigated. At the CCSD(T)/6-311+G(2df,2pd) and G3B3 levels, the global minimum isomer was found to be a linear HCCCCH. The structures of the stable isomer and its relevant transition state are further optimized at the QCISD/6-311G(d,p) level. The bonding nature and structure of isomer HCCCCH were analyzed.     

·C2H3+O2→HC·O+H2CO 的密度泛函理论研究

应用密度泛函理论研究了@C2H3+O2→HC@O+H2CO的反应机理.在DFT(B3LYP/6-31G*)水平上对反应过程中所有反应物、中间体、过渡态和产物的几何构型进行优化,通过频率振动分析确认中间体和过渡态.计算IRC反应路径的能量,分析了中间体的异构化过程和各主要原子的自旋密度.