CH~2(^3B~1)+H~2→CH~3+H的反应途径及动态理论研究

用从头算的UHF/6-31G方法,反应途径哈密顿理论以及变分过渡态理论, 计算了反应CH~2(^3B~1)+H~2→CH~3+H的反应途径.沿反应途径的动态学性质和正则变分过渡态理论的速率常数,用变分过渡态方法处理效果明显;较低温度下考虑隧道效应更重要,而用小曲率近似的校正是有效的;H~2 分子的振动激发对反应速率常数有较大增进.     

NCO+H_2→HNCO+H的动态学和选模反应研究

用从头算方法计算了NCO+H_2→HNCO+H的反应途径。在此基础上,计算了沿反应途径的动态学性质,正则变分过渡态理论的速率常数和振动选模反应的速率常数。结果表明,电子相关能对反应活化位垒影响较大,用多体微扰方法进行校正是有效的;反应存在返回效应和隧道效应,用变分过渡态方法和小曲率近似方法分别进行校正也是有效的。反应中,两个能量较高的束缚振动模式与反应途径有较强的耦合和传能过程,激发这些振动,特别是H_2的伸缩振动,对反应有较大促进。     

CH~2(^3B~1)+H~2→CH~3+H的反应途径及动态理论研究

用从头算的UHF/6-31G方法,反应途径哈密顿理论以及变分过渡态理论, 计算了反应CH~2(^3B~1)+H~2→CH~3+H的反应途径.沿反应途径的动态学性质和正则变分过渡态理论的速率常数,用变分过渡态方法处理效果明显;较低温度下考虑隧道效应更重要,而用小曲率近似的校正是有效的;H~2 分子的振动激发对反应速率常数有较大增进.     

C2H2^＋—2→C2H3^＋＋H反应振动选态速率常数的理论研究

用从头计算法辅助以能量梯度法优化了反应Ｃ２Ｈ２＾＋＋Ｈ２→Ｃ２Ｈ３＾＋＋Ｈ的过渡态，用福井谦一的理论求出反应途径，用反应途径哈密顿理论及正则变分过渡态理论计算以应途径的动力学性质和反应速度数，在此基础上对涉及振动激发的选态反应速率常数进行计算，所得结果与现有的实验结果相符合。     

C_2H_2~＋＋H_2→C_2H_3~＋＋H反应振动选态速率常数的理论研究

用从头计算法辅以能一梯度法优化了反应Ｃ_２Ｈ＋Ｈ_２→Ｃ_２Ｈ＋Ｈ的过渡态，用福井谦一的理论求出反应途径，用反应途径哈密顿理论及正则变分过渡态理论计算沿反应途径的动力学性质和反应速率常数。在此基础上对涉及振动激发的选态反应速率常数进行计算，所得结果与现有的实验结果相符合。     

H2 + CN(n=0,1)→H + HCN振动选态反应

用从头算方法获得了Ｈ２＋ＣＮ反应的内禀反应坐标（ＩＲＣ），沿着ＩＲＣ，计算了各垂直于ＩＲＣ的简正模所对应的频率（Ｗ）以及沿ＩＲＣ运动与垂直ＩＲＣ运动的简正模之间的耦合常数（ＢＫＦ），根据传统过渡态，变分过渡态理论和选态公式，计算了ｎＣＮ＝０及ｎＣＮ＝１时反应的速率常数，并得到了实验相一致的结果，还计算了ｎＣＨ＝１及ｎＣＮ＝１的Ｈ＋ＨＣＮ→Ｈ２＋ＣＮ反应速率常数，可供实验工作者参考。     

H2O(HOD) + Cl→HCl(DCl) + OH(OD)反应动态学的理论研究

用从头算方法, 获得了H2O + Cl→HCl + OH(R1), HOD +Cl→DCl + OH(R2), HOD + Cl→HCl + OD(R3)反应的内禀反应坐标(IRC)。根据传统过渡态、变分过渡态理论及相应的隧道效应校正, 计算了反应的速率常数。对已有实验速率常数值的R1反应, 我们计算的结果和实验一致。根据Truhlar的振动选态公式, 分别讨论了激发HOD中OH, OD振动模式对反应速率的影响,得到激发HOD中的OH振动模式将有利于产物OD + HCl生成, 和实验的结论相一致。     

甲硅硫醛(H2SiS)脱氢反应的反应路径解析

本文用内禀反应坐标(IRC)方法讨论了甲硅硫醛分子的脱氢反应机理. 在RHF/4-31G基组上对此反应做了量子化学从头计算和反应路径解析, 得到H_2SiS的平衡几何、过渡态结构、反应势能曲线、活化能, 反应热以及沿反应坐标一些物理量的变化, 并对平衡几何和过渡态做了振动分析. 同时, 计算了反应的频率因子A和在500 K时的活化熵△S≠。     

CH3S自由基H迁移异构化及脱H2反应的直接动力学研究

采用密度泛函方法（MPW1PW91）在6．311G（d,p）基组水平上研究了CH3S自由基H迁移反应CH3S→CH2SH（R1）,脱H2反应CH3S→HCS＋H2（R2）以及脱H2产物HCS异构化反应HCS→CSH（R3）的微观动力学机理．在QCISD（t）／6．311＋＋G（d,p）／／MPW1PW91／6．311G（d,p）＋ZPE水平上进行了单点能校正．利用经典过渡态理论（TST）与变分过渡态理论（CVT）分别计算了各反应在200-2000K温度区间内的速率常数K^TST和k^CVT,同时获得了经小曲率隧道效应模型（SCT）校正后的速率常数萨k^CVT/SCT．结果表明,反应R1,R2和R3的势垒△E^≠分别为160．69,266．61和241．63kJ／mol。R1为反应的主通道．低温下CH3S比CH2SH稳定,高温时CH2SH比CH3S更稳定．另外,速率常数计算结果显示,量子力学隧道效应在低温段对速率常数的计算有显著影响,而变分效应在计算温度段内对速率常数的影响可以忽略．     

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

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

Seven dimensional quantum dynamics study of the H2+NH2-->H+NH3 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H2+NH2-->H+NH3 reaction using a seven dimensional model on an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values and nonreactive NH2 group keeps C2v symmetry and the rotation-vibration coupling in NH2 is neglected. The total reaction probabilities are calculated when the two reactants are initially at their ground states, when the NH2 bending mode is excited, and when H2 is on its first vibrational excited state, with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants and equilibrium constants are calculated for the temperature range of 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the reaction is dominated by ground-state reactivity and the main contribution to the thermal rate constants is thought to come from this state, (b) the excitation energy of H2 was used to enhance reactivity while the excitation of the NH2 bending mode hampers the reaction, (c) the calculated thermal rate constants are very close to the experimental data and transition state theory results at high and middle temperature, while they are ten times higher than that of transition state theory at low temperature (T=200 K), and (d) the equilibrium constants results indicate that the approximations applied may have different roles in the forward and reverse reactions.     

Full dimensional time-dependent quantum dynamics study of the H + NH3 --> H2 + NH2 reaction

A rigorous full dimensional time-dependent wave packet method has been developed for the reactive scattering between an atom and a tetra-atomic molecule. The method has been applied to the hydrogen abstraction reaction H+NH(3)-->H(2)+NH(2). Initial state-selected total reaction probabilities are investigated for the reactions from the ground vibrational state and from four excited vibrational states of ammonia. The total reaction probabilities from two lowest "tunneling doublets" due to the inversion barrier for the umbrella bending motion of NH(3) and from two pairs of doubly degenerate vibrational states of NH(3) are also inspected. Integral cross sections and rate constants are calculated for the reaction from the ground state with the centrifugal-sudden approximation. The calculated results are compared with those from the previous seven dimensional calculations [M. Yang and J. C. Corchado, J. Chem. Phys. 126, 214312 (2007)]. This work shows that the full dimensional rate constants are a factor of 3 larger than the corresponding seven dimensional calculated values at T=200 K and are overall smaller than those obtained from the variational transition state theory in the whole temperature region. The work also reveals that nonreactive NH bonds of NH(3) cannot be treated as spectators due to the fact that three NH bonds are coupled with each other during the reaction process.     

Thermal and state-selected rate constant calculations for O(3p) + H2 → OH + H and isotopic analogs

We present a new parametrization (based on ab initio calculations) of the bending potentials for the two lowest potential energy surfaces of the reaction O(³P) + H₂, and we use it for rate constant calculations by variational transition-state theory with multidimensional semiclassical tunneling corrections. We present results for the temperature range 250–2400 K for both the rate constants and the intermolecular kinetic isotope effects for the reactions of O(³P) with D₂ and HD. In general, the calculated rate constants for the thermal reactions are in excellent agreement with available experiments. We also calculate the enhancement effect for exciting H₂ to the first excited vibrational state. The calculations also provide information on which aspects of the potential energy surfaces are important for determining the predicted rate constants.     

A product branching ratio controlled by vibrational adiabaticity and variational effects: kinetics of the H + trans-N2H2 reactions

The abstraction and addition reactions of H with trans-N(2)H(2) are studied by high-level ab initio methods and density functional theory. Rate constants were calculated for these two reactions by multistructural variational transition state theory with multidimensional tunneling and including torsional anharmonicity by the multistructural torsion method. Rate constants of the abstraction reaction show large variational effects, that is, the variational transition state yields a smaller rate constant than the conventional transition state; this results from the fact that the variational transition state has a higher zero-point vibrational energy than the conventional transition state. The addition reaction has a classical barrier height that is about 1 kcal∕mol lower than that of the abstraction reaction, but the addition rates are lower than the abstraction rates due to vibrational adiabaticity. The calculated branching ratio of abstraction to addition is 3.5 at 200 K and decreases to 1.2 at 1000 K and 1.06 at 1500 K.     

Seven-dimensional quantum dynamics study of the H+NH3-->H2+NH2 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H+NH3-->H2+NH2 reaction using a seven-dimensional model and an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values. The total reaction probabilities are calculated for the initial ground and seven excited states of NH3 with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants are calculated for the temperature range 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the total reaction probabilities are overall very small, (b) the symmetric and asymmetric NH stretch excitations enhance the reaction significantly and almost all of the excited energy deposited was used to reduce the reaction threshold, (c) the excitation of the umbrella and bending motion have a smaller contribution to the enhancement of reactivity, (d) the main contribution to the thermal rate constants is thought to come from the ground state at low temperatures and from the stretch excited states at high temperatures, and (e) the calculated thermal rate constants are three to ten times smaller than the experimental data and transition state theory results.     

Vibrational energy relaxation of the bend fundamental of dilute water in liquid chloroform and d-chloroform

The population lifetimes of the bend fundamental of dilute water in liquid chloroform (8.5 ps) and d-chloroform (28.5 ps) display an interesting solvent isotope effect. As the lowest excited vibrational state of the molecule, the water bend fundamental relaxes directly to the ground state with about 1600 cm-1 of energy released to the other degrees of freedom. The strong solvent isotope effect along with the large energy gap indicates the participation of solvent vibrational modes in this vibrational energy relaxation process. We calculate the vibrational energy relaxation rates of the water bend in chloroform and d-chloroform using the Landau-Teller formula with a new potential model developed and parametrized self-consistently to describe the chloroform-water interaction. The computed values are in reasonable agreement with the experimental results, and the trend for the isotope effect is correct. It is found that energy transfer to the solvent vibrations does indeed play an important role. Nevertheless, no single dominant solvent accepting mode can be identified; the relaxation appears to involve both the bend and the C-Cl stretches, and frequency changes of all of these modes upon deuteration contribute to the observed solvent isotope effect.     

Reaction rates for O + HD → OH + D and O + HD → OD + H

We have calculated reaction rates for the reactions O + HD → OH + D and O + DH → OD + H using improved canonical variational transition state theory and least-action ground-state transmission coefficients with an ab initio potential energy surface. The kinetic isotope effects are in good agreement with experiment. The optimized tunneling paths and properties of the variational transition states and the rate enhancement for vibrationally excited reactants are also presented and compared with those for the isotopically unsubstituted reaction O + H₂ → OH + H. The thermal reactions at low and room temperature are predicted to occur by tunneling at extended configurations, i.e., to initiate early on the reaction path and to avoid the saddle point regions. Tunneling also dominates the low and room temperature reactions for excited vibrational states, but in these cases the results are not as sensitive to the nature of the tunneling path. Overbarrier mechanisms dominate for both thermal and excited-vibrational state reactions for T > 600 K. For the excited-state reaction (with initial vibrational quantum number n > 0) a transition state switch occurs for T > 1000 K for the O + HD(n = 1) → OD + H case and for T > 1500 K for the O + DH(n = 1) → OD + H reaction, and this may be a general phenomenon for excited-state reactions at higher temperature. In the present case the switch occurs from an early variational transition state where the vibrationally adiabatic approximation is expected to be valid to a tighter variational transition state where nonadiabatic effects are probably important and should be included.     

Theoretical rate constants and kinetic isotope effects in the reaction of methane with H, D, T, and Mu atoms

The rate constants and kinetic isotope effects of the reaction of methane with four isotopes of hydrogen, protium (H), deuterium (D), tritium (T), and muonium (Mu), were studied using variational transition state theory with multidimensional tunneling on an analytical potential energy surface, PES-2002, previously constructed by our group. For the four isotopes, our kinetics results agree reasonably with available experimental measurements, improving previous theoretical results that used different potential energy surfaces and/or theoretical approaches. In the comparison of the reactivity between protium and muonium, which is the most severe test of the surface and theoretical method due to the large mass difference between the two isotopes, some sources of discrepancy between theory and experiment were analyzed. These were the zero-point energy, tunneling effect, and the role of the reactivity from methane excited vibrational states.     

Including anharmonicity in the calculation of rate constants. 1. The HCN/HNC isomerization reaction

A method for calculating anharmonic vibrational energy levels in asymmetric top and linear systems that is based on second-order perturbation theory in curvilinear coordinates is extended to the bound generalized normal modes at nonstationary points along a reaction path. Explicit formulas for the anharmonicity coefficients, x(ij), and the constant term, E0, are presented, and the necessary modifications for resonance cases are considered. The method is combined with variational transition state theory with semiclassical multidimensional tunneling approximations to calculate thermal rate constants for the HCN/HNC isomerization reaction. Although the results for this system are not very sensitive to the choice of coordinates, we find that the inclusion of anharmonicity leads to a substantial improvement in the vibrational energy levels. We also present detailed comparisons of rate constants computed with and without anharmonicity, with various approximations for incorporating tunneling along the reaction path, and with a more practical approach to calculating the vibrational partition functions needed for larger systems.