Capture and dissociation in the complex-forming CH + H2 → CH2 + H, CH + H2 reactions

The rate coefficients for the capture process CH + H(2)→ CH(3) and the reactions CH + H(2)→ CH(2) + H (abstraction), CH + H(2) (exchange) have been calculated in the 200-800 K temperature range, using the quasiclassical trajectory (QCT) method and the most recent global potential energy surface. The reactions, which are of interest in combustion and in astrochemistry, proceed via the formation of long-lived CH(3) collision complexes, and the three H atoms become equivalent. QCT rate coefficients for capture are in quite good agreement with experiments. However, an important zero point energy (ZPE) leakage problem occurs in the QCT calculations for the abstraction, exchange and inelastic exit channels. To account for this issue, a pragmatic but accurate approach has been applied, leading to a good agreement with experimental abstraction rate coefficients. Exchange rate coefficients have also been calculated using this approach. Finally, calculations employing QCT capture/phase space theory (PST) models have been carried out, leading to similar values for the abstraction rate coefficients as the QCT and previous quantum mechanical capture/PST methods. This suggests that QCT capture/PST models are a good alternative to the QCT method for this and similar systems.     

H+CH3NO2H2+CH2NO2反应途径和变分速率常数计算研究

采用MP2(FULL)/6-311G**从头算方法, 优化了H+CH3NO2H2+ CH2NO2反应的过渡态结构, 得出该反应的正逆反应的活化位垒分别是82.73和57.14 kJ*mol-1. 沿IRC分析指出该反应是一个H-H键生成和C-H键断裂的协同反应, 而且在反应途径上存在一个引导反应进行的振动模式, 这一反应模式引导反应进行的区间在-0.7～0.2( amu)1/2*a0之间; 在1 000～1 400 K温度范围内, 运用变分过渡态理论(CVT), 计算了该反应的速率常数, 计算结果与实验相一致.     

Capture and dissociation in the complex-forming CH(v = 0,1) + D2 → CHD + D, CD2 + H, CD + HD reactions and comparison with CH(v = 0,1) + H2

Rate coefficients for the CH(v = 0,1) + D(2) reaction have been determined for all possible channels (T: 200-1200 K), using the quasiclassical trajectory method and a suitable treatment of the zero point energy. Calculations have also been performed on the CH(v = 1) + H(2) reaction and the CH(v = 1) + D(2) → CH(v = 0) + D(2) process. Most of the results can be understood considering the key role played by the deep minimum of the potential energy surface (PES), the barrierless character of the PES, the energy of the reaction channels, and the kinematics. The good agreement found between theory and experiment for the rate coefficients of the capture process of CH(v = 0) + D(2), the total reactivity of CH(v = 1) + D(2), H(2), as well as the good agreement observed for the related CH(v = 0) + H(2) system (capture and abstraction), gives confidence on the theoretical rate coefficients obtained for the capture processes of CH(v = 1) + D(2), H(2), the individual reactive processes of CH(v = 1) + D(2), H(2), the abstraction and abstraction-exchange reactions for CH(v = 0) + D(2), and the inelastic process mentioned above, for which there are no experimental data available, and that can be useful in combustion chemistry and astrochemistry.     

用变分过渡态理论研究CH3SiH3+H→CH3SiH2+H2反应动力学

用变分过渡态理论对CH3SiH3与H的抽提反应进行了理论研究;利用从头算计算了反应体系的构型、振动频率和能量等信息;计算了温度在298 ～1700K内反应的速率常数和穿透系数。结果表明,在室温下,变分对于此反应影响较大,隧道效应特别明显,计算得到的速率常数和实验值符合得很好。     

Ab initio potential energy surface and quantum dynamics for the H + CH4 → H2 + CH3 reaction

A new full-dimensional potential energy surface for the title reaction has been constructed using the modified Shepard interpolation scheme. Energies and derivatives were calculated using the UCCSD(T) method with aug-cc-pVTZ and 6-311++G(3df,2pd) basis sets, respectively. A total number of 30,000 data points were selected from a huge number of molecular configurations sampled by trajectory method. Quantum dynamical calculations showed that the potential energy surface is well converged for the number of data points for collision energy up to 2.5 eV. Total reaction probabilities and integral cross sections were calculated on the present surface, as well as on the ZBB3 and EG-2008 surfaces for the title reaction. Satisfactory agreements were achieved between the present and the ZBB3 potential energy surfaces, indicating we are approaching the final stage to obtain a global potential energy surface of quantitative accuracy for this benchmark polyatomic system. Our calculations also showed that the EG-2008 surface is less accurate than the present and ZBB3 surfaces, particularly in high energy region.     

Bimolecular reaction rates from ring polymer molecular dynamics: application to H + CH4 → H2 + CH3

In a recent paper, we have developed an efficient implementation of the ring polymer molecular dynamics (RPMD) method for calculating bimolecular chemical reaction rates in the gas phase, and illustrated it with applications to some benchmark atom-diatom reactions. In this paper, we show that the same methodology can readily be used to treat more complex polyatomic reactions in their full dimensionality, such as the hydrogen abstraction reaction from methane, H + CH(4) → H(2) + CH(3). The present calculations were carried out using a modified and recalibrated version of the Jordan-Gilbert potential energy surface. The thermal rate coefficients obtained between 200 and 2000 K are presented and compared with previous results for the same potential energy surface. Throughout the temperature range that is available for comparison, the RPMD approximation gives better agreement with accurate quantum mechanical (multiconfigurational time-dependent Hartree) calculations than do either the centroid density version of quantum transition state theory (QTST) or the quantum instanton (QI) model. The RPMD rate coefficients are within a factor of 2 of the exact quantum mechanical rate coefficients at temperatures in the deep tunneling regime. These results indicate that our previous assessment of the accuracy of the RPMD approximation for atom-diatom reactions remains valid for more complex polyatomic reactions. They also suggest that the sensitivity of the QTST and QI rate coefficients to the choice of the transition state dividing surface becomes more of an issue as the dimensionality of the reaction increases.     

将SVRT模型应用到单原子多原子反应F+CH2 D2→CH2 D/CHD2 +DF/HF中

The semirigid vibrating rotor target（SVRT）model is applied for the reaction F+CH2D2→CH2D/CHD2+DF/HF. The time-dependent wave packet approach is also used in the calculation. Reaction probabilities,crosssections,and rate constants are calculated for the title reaction from the ground state of the reagent on the modified J1（MJ1）potential energy surface（PES）for both channels. Numerical calculation shows the oscillatory structures in the energy dependence of the calculated reaction probability. Those structures are generally associated with broad dynamical resonances. They are smooth in the energy dependence of integral cross-sections due to summation over partial waves. The calculated rate constants are in good agreement with the experimental measurement.     

Theoretical Studies on CH3SiH3+H→CH3SiH2+H2 Reaction with the Variational Transitional State Theory

The abstract reaction of CH3SiH3 with H has been studied by using the“direct dynamics”method of variational transition-state theory,which is based on the information on geometries,frequencies and energies calculated by ab initio along the minimum energy path.The rate constants and transmission coefficients were calculated for the temperature range 298～1700K.The result indicates that the variational effect on this reaction is great and the tunneling effect is very obvious at room temperature.The rate constants calculated match well with the experimental value.     

CH2O+H→CHO+H2反应途径和变分速率常数计算研究

采用QCISD/6-311G^** 从头算方法,优化了吸氢反应CH2O+H→CHO+H2的反应物、过渡态、产物几何结构,得出该反应的正、逆反应活化位垒分别是35.4kJ/mol和98.8kJ/mol。沿IRC分析指出该反应是一个C—H键断裂和H—H键生成协同进行的反应,而且在反应途径上存在一个引导反应进行的振动模式,这一反应模式引导反应进行的区间在—0.4～0.55(amu)^1/2之间。在300～3200K温度范围内,运用变分过渡态理论(CVT),计算了该反应的速率常数。     

Mean potential statistical theory of the H+ + D2 → HD + D+ reaction

The reactive collision process H(+) + D(2)(ν = 0, j = 0) → HD + D(+) is theoretically analyzed for collision energies ranging from threshold up to 1.3 eV. It is assumed that the reaction takes place via formation of a collision complex. In calculations, a statistical theory is used, based on a mean isotropic potential deduced from a full potential energy surface. Calculated integral cross sections, opacity functions, and rotational distributions of the HD products are compared with recent statistical and quantum mechanical calculations performed using a full potential energy surface. Satisfactory agreement between the results obtained using the two statistical methods is found, both of which however overestimate the existing quantum mechanical predictions. The effects due to the presence of identical particles are also discussed.     

CH2=CCIFhv→·CH=CCIF+H光解反应的理论研究

采用UHF,CIS和CASSCF方法,在aug-cc-pvdz基组水平上对CH2=CCIFhv→·CH=CClF+H的光解反应通道及其后续反应作了研究.计算表明:分子吸收一个光子后,在第一电子激发态(S1)经过一个过渡态解离与Cl原子同侧的C-H键,这与用CIS方法计算垂直激发得到的π→σ*C-H跃迁及其对Frank-Condon点的计算中分子的单占轨道和键电荷密度变化所预测的结果是一致的.光解产物·CH=CCIF(基态)还可再发生反应,经过渡态解离C-Cl键或是C-F键.     

OH+ C2H2N←C2H3 + NO→CH3 + NCO反应机理的密度泛函理论研究

应用密度泛函理论研究了反应通道(a)C2H3 NO→CH3 NCO和(b)C2H3 NO→OH C2H2N的反应机理．在B3LYP／6-31G(d)水平上优化了反应物、中间体、过滤态、产物的几何构型，通过频率分析确定了11个中间体和10个过渡态．所有的反应物、中间体、过渡态、产物都在CCSD／6-311 G(d，P)水平上进行了单点能较正．并讨论了反应的异构化过程．计算结果表明10是能量最低的中间体，比反应物的能量低308．479kJ／mol；过渡态1／3，2／5，3／4，4／8比反应物的能量高，其中3／4是能量最高的过渡态，比反应物的能量高91．894kJ／mol．通道(a)和(b)的理论放热值分别为111．059和96．619kJ／mol．     

The distonic ion ·CH2CH2CH+OH,keto ion CH3CH2CH=O +·, enol ion CH3CH=CHOH+·, and related C3H6O+· radical cations. Stabilities and isomerization proclivities studied by dissociation and neutralization-reionization

Metastable ion decompositions, collision-activated dissociation (CAD), and neutralization-reionization mass spectrometry are utilized to study the unimolecular chemistry of distonic ion ^·CH₂CH₂CH^?OH (2^+·) and its enol-keto tautomers CH₃CH=CHOH^?· (1 ^+·) and CH₃CH₂CH=O ^+· (3^+·). The major fragmentation of metastable 1^+·–3^+· is H^· loss to yield the propanoyl cation, CH₃CH₂C≡O⁺. This reaction remains dominant upon collisional activation, although now some isomeric CH₂=CH-CH⁺ OH is coproduced from all three precursors. The CAD and neutralization-reionization (⁺NR⁺) spectra of keto ion 3 ^+· are substantially different from those of tautomers 2^+· and 1^+·. Hence, 3^+· without sufficient energy for decomposition (i. e. , “stable” 3^+·) does not isomerize to the ther-modynamically more stable ions 2^+· or 1^+·, and the 1,4-H rearrangement H-CH₂CH₂CH=O^+·(3 ^+·) → CH₂CH₂CH⁺ O-H (2 ^+·) must require an appreciable critical energy. Although the fragment ion abundances in the ⁺ NR ⁺ (and CAD) spectra of 1 ^+· and 2 ^+· are similar, the relative and absolute intensities of the survivor ions (recovered C₃H₆O^+· ions in the ⁺NR⁺ spectra) are markedly distinct and independent of the internal energy of 1 ^+· and 2 ^+·. Furthermore, 1 ^+· and 2 ^+· show different MI spectra. Based on these data, distonic ion 2 ^+· does not spontaneously rearrange to enol ion 1 ^+· (which is the most stable C₃H₆O^+· of CCCO connectivity) and, therefore, is separated from it by an appreciable barrier. In contrast, the molecular ions of cyclopropanol (4 ^+·) and allyl alcohol (5 ^+·) isomerize readily to 2 ^+·, via ring opening and 1,2-H^? shift, respectively. The sample found to generate the purest 2 ^+· is α-hydroxy-γ-butyrolactone. Several other precursors that would yield 2 ^+· by a least-motion reaction cogenerate detectable quantities of enol ion 1 ^+·, or the enol ion of acetone (CH₂=C(CH₃)OH^+·, 6 ^+·), or methyl vinyl ether ion (CH₃OCH=CH ₂ ^+· , 7 ^+·). Ion 6 ^+· is coproduced from samples that contain the —CH₂—CH(OH)—CH₂— substructure, whereas 7 ^+· is coproduced from compounds with methoxy substituents. Compared to CAD, metastable ion characteristics combined with neutralization-reionization allow for a superior differentiation of the ions studied.     

Infinite order sudden approximation treatment of the H + D2 → HD + D reaction

Results of a study of the H + D₂ → HD + D reaction within the quantum infinite order sudden approximation (RIOS) are reported here for the total energy range 1.0317 ⩽ E_t ⩽ 2.1417 eV. We present the vibrationally selected integral and total integral cross sections, and the latter are compared with both classical and experimental results. The RIOS results are in excellent agreement with the experimental results at 1.0317 eV and in reasonable agreement (≈ 20% high) at 2.1417 eV.     

CH2+O2反应的反应机理

CH2与 O2的反应是不饱和碳氢化合物燃烧过程中的一个十分重要的反应 .CH2（基态,三重态）与 O2反应的实验研究工作已有不少报导 [1－ 5],其主要产物有两 组 :H2O＋ CO和 CO2＋ H2,这些产物表明反应在反应体系的单重态位能面上发生 .至今还未见到关于 CH2＋ O2反应机理的完整的理论研究报导 .本文用 CASSCF方法详细研究了 CH2＋ O2反应的机理,给出从反应物至两组不同最终产物的完整的反应途经的描写 .1计算方法 　　计算使用量子化学高斯 98 W软件 .用从头算 CASSCF(8,8)/6-31G(d,p)方法在 CH2＋ O2单重态位能面上以优化…     

Nonadiabatic quantum dynamics of C(1D)+H2→CH+H: coupled-channel calculations including Renner-Teller and Coriolis terms

The Renner-Teller (RT) coupled-channel dynamics for the C((1)D)+H(2)(X(1)Σ(g) (+))→CH(X(2)Π)+H((2)S) reaction has been investigated for the first time, considering the first two singlet states ??(1)A' and b(1)A' of CH(2) dissociating into the products and RT couplings, evaluated through the ab initio matrix elements of the electronic angular momentum. We have obtained initial-state-resolved probabilities, cross sections and thermal rate constants via the real wavepacket method for both coupled electronic states. In contrast to the N((2)D)+H(2)(X(1)Σ(g)(+)) system, RT effects tend to reduce probabilities, cross sections, and rate constants in the low energy range compared to Born-Oppenheimer (BO) ones, due to the presence of a repulsive RT barrier in the effective potentials and to long-lived resonances. Furthermore, contrary to BO results, the rate constants have a positive temperature dependence in the 100-400 K range. The two-state RT rate constant at 300 K, lower than the BO one, remains inside the error bars of the experimental value.     

CH2+O2反应的反应机理

The mechanisms of the CH2＋ O2→ H2O＋ CO and CH2＋ O2→ H2＋ CO2 reactions have been studied by performing ab initio CAS(8,8)/6-31G(d,p) calculations, and five intermediates(IMn) and eight transitions(TSn) have been located along the reaction paths. The predicted path for the CH2＋ O2→ H2O＋ CO is: CH2＋ O2→ TS1→ IM1→ TS2→ IM2→ TS3→ IM3→ TS4→ IM4a→ TS5→ H2O＋ CO. For the CH2＋ O2→ H2＋ CO2 reaction, there are two paths: (i) CH2＋ O2→ TS1→ IM1→ TS2→ IM2→ TS3→ IM3→ TS6→ H2＋ CO2 and (ii) CH2＋ O2→ TS1→ IM1→ TS2→ IM2→ TS3→ IM3→ TS4→ IM4a→ TS7→ IM4b→ TS8→ H2＋ CO2, with the latter path more favorable energetically.     

BaCl2-C6H12O6-H2O/(H2O+CH3CH2OH)体系在30 ℃时的相平衡研究

用等温溶解度法研究了BaCl2-C6H12O6-H2O/(H2O+CH3CH2OH)体系在30 ℃时的相平衡, 发现在体系中有一未见文献报道的不一致溶解的化合物BaCl2*2C6H12O6*2H2O生成. 研究结果为肌醇生产除杂质工艺提供了理论依据.