An implementation of the variational quantum RRKM program is presented to utilize the direct ab initio dynamics approach for calculating k(tE, J), k(E) and k(T) within the framework of the microcanonical transition st...An implementation of the variational quantum RRKM program is presented to utilize the direct ab initio dynamics approach for calculating k(tE, J), k(E) and k(T) within the framework of the microcanonical transition state (μTST) and microcanonical variational TST (μVT) theories. An algorithm including tunneling contributions in Beyer-Swinehart method for calculating microcanonical rate constants is also proposed. An efficient piece-wise interpolation method is developed to evaluate the Boltzmann integral in calculation of thermal rate constants. Calculations on several test reactions, namely the H(D)2CO→ H(D)2 + CO, CH2CO→CH2 + CO and CH4 + H→CH3 + H2 reactions, show that the results are in good agreement with the previous rate constants calculations. This approach would require much less computational resource.展开更多
Variational transition state theory is used to study the kinetics of the reactions Li+HF and Li+HCl on the BO potential energy surfaces, fitted to ab initio values by Lagana et al. The saddle point, the location of ge...Variational transition state theory is used to study the kinetics of the reactions Li+HF and Li+HCl on the BO potential energy surfaces, fitted to ab initio values by Lagana et al. The saddle point, the location of generalized transition state and other kinetic features on the bottleneck of the po- tential energy surface are obtained. We point out that the location of generalized transition state tends to the direction in which the two atoms with larger reduced mass separate off. It can be determined whether a potential energy surface is accurate or not by a fine analysis of its bottleneck region and a comparison of kinetic calculation with experimental result. This is helpful for building a more accurate potential energy surface.展开更多
The reaction of C3H8＋O（^3p）→C3HT＋OH is investigated using ab initio calculation and dynamical methods. Electronic structure calculations for all stationary points are obtained using a dual-level strategy. The geo...The reaction of C3H8＋O（^3p）→C3HT＋OH is investigated using ab initio calculation and dynamical methods. Electronic structure calculations for all stationary points are obtained using a dual-level strategy. The geometry optimization is performed using the unrestricted second-order Moller-Plesset perturbation method and the single-point energy is computed us- ing the coupled-cluster singles and doubles augmented by a perturbative treatment of triple excitations method. Results indicate that the main reaction channel is C3Hs＋O（^3p）→i- C3HT＋OH. Based upon the ab initio data, thermal rate constants are calculated using the variational transition state theory method with the temperature ranging from 298 K to 1000 K. These calculated rate constants are in better agreement with experiments than those reported in previous theoretical studies, and the branching ratios of the reaction are also calculated in the present work. Furthermore, the isotope effects of the title reaction are calculated and discussed. The present work reveals the reaction mechanism of hydrogenabstraction from propane involving reaction channel competitions is helpful for the understanding of propane combustion.展开更多
基金supported by the Pre-research Special Project(Grant No.2001CCD03500)of the 973 Fundamental Research Project.
文摘An implementation of the variational quantum RRKM program is presented to utilize the direct ab initio dynamics approach for calculating k(tE, J), k(E) and k(T) within the framework of the microcanonical transition state (μTST) and microcanonical variational TST (μVT) theories. An algorithm including tunneling contributions in Beyer-Swinehart method for calculating microcanonical rate constants is also proposed. An efficient piece-wise interpolation method is developed to evaluate the Boltzmann integral in calculation of thermal rate constants. Calculations on several test reactions, namely the H(D)2CO→ H(D)2 + CO, CH2CO→CH2 + CO and CH4 + H→CH3 + H2 reactions, show that the results are in good agreement with the previous rate constants calculations. This approach would require much less computational resource.
基金Project supported by the National Natural Science Foundation of China.
文摘Variational transition state theory is used to study the kinetics of the reactions Li+HF and Li+HCl on the BO potential energy surfaces, fitted to ab initio values by Lagana et al. The saddle point, the location of generalized transition state and other kinetic features on the bottleneck of the po- tential energy surface are obtained. We point out that the location of generalized transition state tends to the direction in which the two atoms with larger reduced mass separate off. It can be determined whether a potential energy surface is accurate or not by a fine analysis of its bottleneck region and a comparison of kinetic calculation with experimental result. This is helpful for building a more accurate potential energy surface.
文摘The reaction of C3H8＋O（^3p）→C3HT＋OH is investigated using ab initio calculation and dynamical methods. Electronic structure calculations for all stationary points are obtained using a dual-level strategy. The geometry optimization is performed using the unrestricted second-order Moller-Plesset perturbation method and the single-point energy is computed us- ing the coupled-cluster singles and doubles augmented by a perturbative treatment of triple excitations method. Results indicate that the main reaction channel is C3Hs＋O（^3p）→i- C3HT＋OH. Based upon the ab initio data, thermal rate constants are calculated using the variational transition state theory method with the temperature ranging from 298 K to 1000 K. These calculated rate constants are in better agreement with experiments than those reported in previous theoretical studies, and the branching ratios of the reaction are also calculated in the present work. Furthermore, the isotope effects of the title reaction are calculated and discussed. The present work reveals the reaction mechanism of hydrogenabstraction from propane involving reaction channel competitions is helpful for the understanding of propane combustion.
